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that reflects reality. For an underwater virtual world that simulates acoustic detection, a physically based 
sonar propagation model is required if ranges in excess of tens of meters are expected. 

This thesis creates an application programming interface (API) for real-time 3D computation 
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The simulation is programmed in Java and runs either as a stand-alone program or as a script in a web 
browser. This program generates Virtual Reality Modeling Language (VRML 97) compliant code that 
can be viewed from any VRML-capable browser. This approach allows the characteristics of the energy 
propagation to be calculated with high precision and observed in 3D. 

As sonar system information bandwidth becomes larger, more intuitive ways of presenting 
information to a user will be required. Higher information density in a more intuitive format can free the 
user from integrating the data himself and allow quicker reaction times. This thesis and the API provide 


the foundation for fundamental advances in sonar modeling and visualization. 
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I. INTRODUCTION 


A. OVERVIEW 


This thesis provides an interactive three-dimensional (3D) sonar propagatıon and rendering 
application programming interface (API) for use in virtual worlds as well as stand-alone simulation and 
visualization programs. The interactive program allows a simulation to send virtual sonar energy pulses 
into a virtual world for detection of target objects. This allows the simulation to have real-time and 
realistic sonar information from the virtual world. The stand-alone programs allow the compilation of 
data on how the sonar pulse energy interacts with the ocean environment. These programs allow the user 
to view the data as 3D virtual worlds. These worlds allow the researcher the ability to immerse himself in 
the data and develop a deeper understanding of the data set. The accuracy of this physically based sonar 
model is very comparable to other widely مم‎ sonar propagation models. Its accuracy is achieved by 
using the recursive Tay acoustics (RRA) algorithm [Ziomek, 1996]. 

The Phoenix autonomous underwater vehicle (AUV) is a research robot at the Naval 
Postgraduate School (NPS). It is designed to navigate autonomously and perform a variety of tactical 
functions in the ocean. In order to perform complex testing without getting the robot ^wet," a virtual 
world for the robot to interact with was created. This world responds to the robot's control signals, 
allowing it to navigate in the virtual world, and it provides external sensor stimulation information to the 
robot for feedback and information collection. In this type of world two models are highly important: the 
hydrodynamics and the sonar. The hydrodynamics model is a full-featured six degree-of-freedom model 
which accounts for cross body-flow. The current sonar model is a straight-line geometric model. In 
essence, the model depicts acoustic energy as traveling in a conical beam through the ocean. This sonar 
model is satisfactory for distances on the order of a few hundred meters, but for ranges on the order of 


kilometers will not provide realistic results. 


B. MOTIVATION 


Real time simulation of sonar information and visualization of that information are important for 
several reasons. The use of a sonar simulator for training of crews and AUV's is enhanced by physically 
based, realistic sonar information. In this instance not only must the data be realistic but it also must be 
real time. Tactical planning and visualization is another field that benefits from realistic 3D sonar 
information. If an area of ocean can be accurately simulated then more detailed plans can be made for the 
searching and controlling of that area. The ability to generate an API that can be used in these capacities 


is the motivation behind this thesis. 


e: PURPOSE 


The purpose of this thesis is to provide an application programming interface (API) for use in 
investigating the propagation of acoustic energy in shallow water regions of the ocean. This API is to 
provide a convenient and intuitive interface to the user to allow him to easily integrate physics based 
ocean acoustics into a simulation. The API will provide methods for the user to specify a virtual world of 
arbitrary complexity, propagate acoustic energy in the world, collect information on any objects 
encountered and visualize the energy as it propagates in the world. The user is able to vary all aspects of 


the virtual world and of the acoustic energy pulse. 


D. ORGANIZATION 


This thesis provides all the information necessary to understand the motivation for accurate sonar 
simulations, the construction of the sonar model, the implementation of the model in Java, 3D 
visualization and the accuracy of the model compared to other generally accepted models. 

Chapter II, Related Work, provides information about the motivation for this thesis. It describes 
previous work in the fields of tactical, virtual world and sonar modeling and simulation. Chapter HI, 
Problem Statement, relays in detail the problem to be solved, possible solutions and the chosen solution. 
Chapter IV, Recursive Ray Acoustics (RRA) Derivation, describes the development of the ray acoustic 


2 


approximation to the linear acoustic wave equation. It develops the equations that govern the propagatıon 
of acoustic energy in the sonar model. Chapter V, Modeling of Recursive Ray Acoustics (RRA), describes 
each sub-model in the RRA, recursive ray acoustics, model and the methods and data required in each of 
these sub-models. Chapter VI, Sonar Visualization, provides the same information on visualization that 
Chapter V provided on RRA. Chapter VII, Model Implementation and Integration, discusses the method 
and vehicle of implementation of the RRA and visualization models. Chapter VIII, Simulation Results, 
compares time, position and energy transport results from the model to other generally accepted models. 
Chapter IX, Conclusions and Recommendations, summarizes the thesis and provides information about 


improvements and extensions that can be made to the simulation. 


aa 
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Il. RELATED WORK 


A. INTRODUCTION 


Many past works play an important role in the development of thıs thesis. They are summarized 
in this chapter in order to provide the motivation and background needed by the reader. The Phoenix 
AUV virtual world is the original project motivating this thesis. The linearly geometric sonar model of 
the virtual world needs to be made more robust so that real-world results can be obtained. In the Manta 
UUV project, tactical visualization simulations show that multiple mine clearing tactics can be evaluated 
both analytically and visually to find an optimal tactic. This project also lacked a realistic sonar model 
and again provided motivation to pursue a real-time sonar simulation. The Manta UUV concept vehicle 
shows that implementation of an autonomous unmanned vehicle is seriously being considered from a 
military stand point and the requirements placed on this type of vehicle helped to refine the features of the 
shallow water sound propagation model. An essential overview of acoustic modeling [Etter, 1996] 
presents the advantages and disadvantages of all types of sonar models. The model evaluation criteria 
provided in Etter [1996] show that the choice of RRA is equal to or better than all other choices, based on 
the criteria of environment complexity, range of propagation and calculation speed. Seminal work by 
Stewart demonstrates the importance of 3D visualization and the derivation of target solutions from 3D 


data via the stochastic back-projection algorithm [Stewart, 1988]. 
B. AUV VIRTUAL WORLD 


1. What It Is 


The NPS AUV virtual world, [Brutzman, 1994] was created to provide a realistic “out of the 


water” simulation environment for the Phoenix AUV. The Phoenix AUV was created at the Naval 





Postgraduate School in Monterey, California and a description of the vehicle is given in the dissertation. 


Figure 2.1 is a drawing of the Phoenix AUV. 
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Figure 2.1 External drawing of the Phoenix AUV 


The basic idea behind virtual simulation is that every time a new change is made to the hardware 
or software of the AUV, testing of the change can be made without the time, expense and inherent danger 
of putting the vehicle in the water. The simulation models all aspects of the effects of the vehicles control 
surfaces on its position, orientation, speed and rate of change of orientation. In addition it models several 
aspects of the ocean environments effects on the AUV. For example, the effects of cross vehicle fluid flow 
are simulated. By attaching the computer in the Phoenix AUV to the computers simulating the ocean 
environment, simulated operation of the AUV in water is achieved. Previous papers [Brutzman, February 
1997] and [Brutzman, 1994], have been devoted to validating the virtual environment against actual 


operation of the AUV in the water. More information on the Phoenix AUV is found in [Brutzman, 1998]. 


2: What Are Its Acoustic Limitations 


The original modeling of sonar detection in the virtual world simulation was rather simplistic. 
Sonar geometry of the test tank was trigonometrically derived for any posture inside the tank. This 
approach is insufficient for most applications, but for the initial abilities of the Phoenix AUV the sonar 
was acceptable. The more general model derived in [Davis, 1996] can handle arbitrarily large and 
cluttered environments. In essence the sonar propagation was just a sonar beam comprised of several 
straight-line ray segments. These ray segments were “shot” into the 3D graphics scene database which 
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comprised the virtual world and any object that they encountered was detected. Figure 2.2 shows the 
AUV in a virtual test tank, with sonar rays coming from its virtual sonar and lines of fluid flow coming 


from its cross-body thrusters and propellers. This is an accurate solution within a few hundred meters. 





Figure 2.2 Virtual Phoenix AUV in a virtual test tank. The red lines 
are sonar rays and the green lines are water flow from the 
propellers and cross-body thrusters. (Brutzman, 1994] 


For several reasons this scheme will not work for long-range propagation. As is described in 
later chapters, sound speed in the water varies tremendously and because of this variability the rays 
undergoing long-range sound propagation do not travel in straight lines. This effect typically is noticeable 
(i.e. > 1 meter in variation) at ranges on the order of hundreds of meters. Also when distances from the 
source of the sound get large enough, the area front of the beam tube gets large. If this area starts to get 
larger than the object the sonar beam is trying to detect, missed detections occur in the numerical 
modeling. The reason for missing returns in a ray-trace model is that the rays may surround but not 
intersect the target even though the beam front might intersect. This poses a serious problem if simulated 
detections are required at ranges on the order of a kilometer or so, depending on the size of the target and 


the beam width of the sonar beam. Improved sonar modeling to cope with such scenarios is the next 
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logical step for the Phoenix AUV sonar model, if Phoenix is to realize its role as an autonomous mine- 


hunting platform. 


C. AUV TACTICAL MINEFIELD SEARCH 


1. What It Is 


The AUV tactical minefield search program is another simulation developed by students 
working with Dr. Brutzman. The program is detailed in [Brutzman, February 1997], where the efficiency 
of nine different minefield search maneuvers is presented. In the case of clearing a mine field with a 
presupposed mine distribution, having a good idea of what search tactics produce the best results is 
essential. If an exhaustive search of the area is desired then all tactic types will produce essentially the 
same results. However, if there is a time or spatial constraint placed on the situation, not all tactics will 
produce the same results. The tactical minefield search simulation shows and analytically compares search 
effectiveness to determine optimal tactics. Though simple in form, it goes a long way toward showing that 


computer simulations can give meaningful structure to the amorphous question of which tactic is best for 


a variety of situations. Figure 2.3 shows a tactical search in progress. 





Figure 2. 3 Simulation - tactics evaluation for mine hunting. This 3D rendering allows the tactician to see the 
simulated search in progress. The search vessel is in the foreground with an oval sonar beam front 
approaching mine targets 








2 What Are Its limitations 


Again the limitations of this simulation apply mostly to the area of sonar detection. In this 
simulation the area to be searched is split into a two-dimensional grid of smaller areas. The areas are 
chosen so that an AUV traversing the area has a high probability of detecting a mine in the area. No 
account for sound speed profile (SSP), search depth or bottom composition was made. To get a truly 
accurate indication of efficiency of detection, these are parameters for which account must be made. For 
example, if there existed a strong surface layer and the AUV is in that surface layer, almost all of the 
energy projected by the AUV stays in the layer. Any mines that are below this layer may remain 


undetected, giving the particular search tactic a lower efficiency than indicated. 


D. MANTA CONCEPT 


1. What It Is 


Manta is a concept unmanned underwater vehicle (UUV) that is being examined by engineers at 
the Naval Undersea Warfare Center (NUWC) Newport, Rhode Island. The concept of operations has four 
Mantas attach to a submarine hull for transit and then detach to augment the submarine capabilities when 
in the area of operations. Figure 2.4 shows a VRML rendering of the Manta UUV. The role of this 
vehicle is still in development but may include mine detection, surveillance and و"‎ eme of hostile 
targets. In any of these situations sonar is important for target acquisition, target tracking and navigation. 
While there exist many sonar models that are used for sonar performance evaluation and off-line 
interaction between vehicle and environment, there are few (if any) models that provide accurate real- 
time active acoustic simulation of the ocean environment. Such a capability is very important as money to 
develop and test new concepts remains scarce. Research time and resources can be utilized more 
efficiently if the initial mechanical, acoustical and tactical testing of the autonomous vehicle are 


performed with high fidelity in the laboratory environment. 
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Figure 2. 4 Manta Concept Vehicle, à) is starboard beam and b) is above the forward port quarter 


2- Why Is It Important 


Since the break up of the Soviet Union, the U.S. Navy's strategic role has shifted from deep 
water operations to a shallow water littoral environment. As several conflicts in recent years have shown. 
minefields can have a devastating cffect on the projection of naval power ashore as well as a devastating 
effect on the ships that encounter them. While the traditional means of mine sweeping work to a limited 
degree. in most cases a covert mapping and sweeping ability might enable rapid neutralization of a 
minefield. The Manta vehicle is the first look at an unmanned underwater vehicle (UUV) capable of 
operating in the littoral environment. While submarines can and have worked in the littoral 
environment, the potential cost both in loss of human life and monetary loss has made this option very 
undesirable. Although a fullv functioning Manta is not cheap. it is less expensive than a. submarine 
because there are no operators on board. Developing such a capability is invaluable and inevitable. More 


information on Manta minefield searching is found in [Brutzman, April 1998]. 
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E. SONAR MODELING AND VISUALIZATION 


sonar modeling attempts to recreate the highly variable ocean acoustic environment in an 
abstract form. Many different sonar models exist primarily because the real world is complex and no one 
model simplification is adequate to describe the real world in a reduced form. Sonar visualization applies 
scientific visualization techniques to the results of the sonar simulation, with the goal of enhancing the 
understanding of the interested researcher. This section identifies related work in sonar modeling and 


visualization. 


1. Sonar Modeling 


Deciding on the proper sonar model is very important and highly dependent on how the model is to be 
used. [Etter, 1996] states that no one sonar model is adequate to describe all situations and that a model 
(or combination of models) must be chosen to fit the particular set of environmental and operational 
conditions. Etter thus identifies a model hierarchy that can be used to systematically define the model to 


be used in a specific situation. Figure 2.5 shows Etter's modeling hierarchy. 
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Figure 2. 5 Sonar Modeling Hierarchy [Etter, 1996] 
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The figure shows how the model proceeds from general to specific as higher and higher layers of 
abstraction are applied. The lowest level of the model deals with the specific geographic and physical 
constraints in the problem to be studied. The next layer identifies the particular propagation routine used 
to transport acoustic energy through the environmental model defined on the lower level. On top of this 
layer are built models of noise and reverberation that directly integrate into the specific acoustic 
propagation model. The last layer of the overall sonar solution is the performance model, which receives 
accurate simulated sonar information from the basic acoustic model and performs signal processing on the 
data. Etter then examines many example models for each layer in the hierarchy. The number of 
permutations of this hierarchy is staggeringly large, which accounts for the relatively few attempts made 
to form a cohesive model. Some comprehensive model operating systems exist that address these issues of 
complexity and interoperability by switching among models, but they lack scalability and generality. Few 
(if any) models are computationally tractable in real time. The overall message in [Etter, 1996] is that the 
ocean is a complex acoustic environment and that models must be carefully considered to ensure realistic 
results. The implication for real-time 3D sonar modeling and visualization is that model selection is 


challenging and critically important. 


2: Sonar Visualization 


{Stewart, 1988] presents a largely statistical approach to modeling, analyzing and visualizing 
sonar return data. His stochastic back-projection technique adaptively forms a visual model of the world 
from inherently noisy sonar information. This sonar data is characterized as high bandwidth, high noise 
and highly redundant. It is redundant in the fact that many propagation paths are available for acoustic 
energy in the ocean environment. The particular advantage of this model is that immediate results are 
always available and that the accuracy of the visualization increases as time progresses. 

Extensive research in the area of sonar visualization shows that a lack of other resources are 
available. As evidence of this statement, a five year literature search of the Acoustical Society of America 


Journal (JASA), IEEE Transactions on Visualization and Computer Graphics, and most published books 
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on scientific visualization revealed only a handful of diagrams rendered in 3D. It is proposed that real- 


time 3D sonar visualization is an important new field of study lacking computational foundation. 


F. SUMMARY 


This chapter reviewed the work related to the design of a shallow water sonar propagation and 
visualization model. Several current virtual world simulations and tactical simulations were presented, as 
well as two systematic approaches to sonar and visualization modeling. For those interested in tactical 
search evaluation, Appendix D and Appendix E contain an annotated slide set and HTML pages on 
tactical visualization using the Manta UUV. Specific faults of various models were addressed and a 
strategy for addressing these faults was presented. In particular [Etter, 1996] describes a systematic 
technique for the development of an accurate and realistic sonar model. Finally, real-time 3D sonar 


visualization is proposed as an important new field of study. 
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HI. PROBLEM STATEMENT 


A. INTRODUCTION 


Chapter II, Related Work, showed that there are several programs and simulations in existence 
that have limited or no realistic sonar modeling. When a virtual world is created to mimic the real world, 
an account of all of the important physical parameters must be made. Performing an underwater acoustic 


simulation with an unrealistic acoustic model ensures unrealistic results. 


B. PROPOSED FRAMEWORK 


A sonar model that has its origins in the physics of underwater acoustics is necessary for any 
high-resolution virtual world simulation of the real world. All of the parameters vital to underwater 
sound propagation must be included: bottom type, water column sound speed profile (SSP) and surface 
interactions. Due to the complexity of acoustic energy propagation in real-world underwater 
environments, the sonar model solution must possess robust visualization capabilities as well as the ability 
to deliver the vast amount of information in the simulation to an outside client. The sonar model needs to 
be realistic in energy transfer through the water and must also provide real-time response for user 


visualization and robot interaction. 


C. FACTORS AFFECTING THE SOLUTION 


The problem at hand is to generate a real-time, three-dimensional (3D) acoustic energy 
propagation simulation. This is in response to the need for autonomous underwater vehicles that can 
perform tactical operations in the underwater environment. This environment is highly complicated, and 
therefore acoustic simulation in other than a physically rigorous manner is of little use. Additionally, due 
to the complexity of the real-world environment, it 1s desirable to build realistic physics-based virtual 


worlds. Why is a virtual world desirable? The development of a virtual world provides many features. 
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First, visualization of the acoustic energy propagation in a spatially oriented manner becomes possible. 
Second, since the scene can be animated, visualization of the time-dependent nature of the acoustic wave 
front can be seen as it develops in both time and space. Lastly, when the sonar model is allowed to run in 
an interactive mode as a server to some client program, then more complex real-world simulations can be 
created. In this last form, the ability to enhance the Phoenix AUV virtual world and the AUV search 
tactics evaluation simulation programs results in data that is more consistent with the real world and its 
processes. However, physics-based real-time solutions to any wave-energy propagation model requires 
sufficient computing power. Surprisingly, adequate computational power can be provided by currently 
available personal computers to give the amount of resolution required for sufficiently accurate results. To 
increase the computational power further, either throughput in an individual processor needs to be 
increased or the number of processors working on the job needs to be increased. Both of these options are 
viable and need to be explored. Increasing throughput is a job being handled quite well by microprocessor 
manufacturers, so no further discussion is warranted here. Increasing the number of processors falls into 
one of two categories: building a massively parallel computer or interconnecting many individual 
computers via a high speed local-area network (LAN). For most researchers , the second option may 
already exist or can be created relatively easily. With the real-time 3D sonar visualization problem well 


defined, how is the solution to be obtained? 


D. POSSIBLE SOLUTIONS 


This section addresses many possible ways of solving the problem. Several different sonar 
propagation models exist, such as ray tracing, normal mode, parabolic equation and others. In addition 
many different ways of simulation exist; direct analog simulation and digital computer simulation are but 
two of the many ways to simulate a physical process. Once the specifics of the implementation are worked 
out, what to produce as the results of the simulation must be decided. In other words, what useful function 
does the simulation provide that was not available before implementation. These aspects and others are 


explored in this section. 
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l Available Sonar Models 


Deciding on the proper sonar model is very important and very dependent on how the model is to 
be used. [Etter, 1996] Chapter 4 is devoted to explaining the advantages and disadvantages of the 
common sonar-propagation models. The chapter presents ray models, normal mode models, multipath 
expansion models and parabolic equation models. A short description of each of the model types is given 
in the following paragraphs. 

Ray theory models start with the Helmholtz equation and make the assumption that a 
propagating wave can be thought of as a surface of constant phase emerging from a source. This surface, 
called the eikonal, can be represented by rays that travel perpendicular to the surface of the wave front. 
The direction of travel of these rays depends completely on the initial direction of travel and the sound 
speed profile (SSP) in the water column. Ray theory requires that the frequency of the acoustic signal be 
large compared to the rate of change of sound speed in the water column and that the wavelength be small 
compared to the characteristic length of structures on the ocean surface and bottom. [Etter, 1996] 
develops Equation (3.1), 


y تن‎ HE (3.1) 


where f is the minimum accurate frequency, c 1s the speed of sound in the ocean and H is the depth of the 
ocean. This gives a good rule of thumb for the low limit of the accuracy of the ray-tracing algorithm. 
Another important characteristic of a sonar modeling theory is that of range dependence. Ray theory 
models are generally range-dependent models. This means that the acoustic structure of the ocean can 
vary with distance from the source in a 2D or 3D manner. 

Normal mode methods are derived from an integral representation of the linear acoustic wave 
equation. To arrive at practical solutions, cylindrical symmetry is assumed which limits range 
dependence to two dimensions. This پا‎ is then solved using separation of variables. This technique 
develops a solution that depicts acoustic waves as traveling waves in the horizontal direction and standing 


waves in the vertical. Many times in order to develop a closed form solution to the normal mode problem, 
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no dependence is placed on range, this is termed range independence. For many shallow water regions, 
especially near shelf break zones, range independence is an invalid assumption. 

Multipath expansion methods expand the acoustic field integral representation into an infinite 
orthogonal set of integrals. This sounds like a bad idea; to go from one integral to an infinite set of 
integrals. However, each integral is simpler than the original and with the proper approximations only a 
certain number are used in the final solution. By integrating only a limited number of integrals, an angle- 
limited source model is produced. This solution technique has many similarities to ray theory, but is 
better suited to deep water. One advantage over ray theory is that caustic and shadow zones are properly 
computed, instead of incorrectly producing infinite pressure fields as ray theory can. 

Parabolic equation methods find their roots in the 1940’s when it was first applied to long-range 
radio propagation. The first application to acoustic waves is in [Hardin and Tappert, 1973]. The basic 
theory is to make an approximation that there is little back scattering of acoustic energy in the water 
column. This approximation transforms the elliptic wave equation into the parabolic wave equation. 
From the form of the parabolic wave equation, many other assumptions about symmetry and sound speed 
profile can be made to simplify the mathematics. Several different numerical integration techniques are 
then used to evaluate the parabolic equation form of the wave equation. Two of the most popular are the 
split-step Fourier algorithm and the finite-element method. The split-step Fourier method has the distinct 
advantage of speed of computation, while the finite-element technique is typically more accurate at higher 
propagation angles. 

[Etter, 1996] provides a table that summarizes the applicability and practicability of each of the 
models in relation to water depth, source frequency and environmental dependence on range. A summary 
of the table from the book is presented in Figure 3.1. As was described earlier, the focus of this model and 
simulation is for the shallow water environment and so the deep-water portion of the figure was left off. 
Therefore, the two independent variables in this table are frequency and range dependence. A range- 
dependent environment has fundamental characteristics that change as the distance from the source 
varies, whereas a range-independent environment does not. Given that the goal is to develop a shallow- 


18 


water sonar model that is robust in a range-dependent environment and which works well with a 
frequency in the 1000 Hz range and higher, then Figure 3.1 indicates that a ray-theory model is best 
suited. If however the frequencies of interest are below 500 Hz and range dependence is required, the 


parabolic equation model is best. 


Applications Low Freq. - « 500 Hz 
Model Type Shallow Water High Freq - > 500 Hz 


RI - Range Independent 
RD- Range Dependent 


RD [RI [RD | 
Normal Mode ee | | 9 | 
Multipath | O | O | © j O | ©@ Applicable and Practical 
Fast Field 8۵ | ا ۵ | 9 | ن‎ (D) Limitations in Accuracy or Speed 
۰.801 0 ۱ 01 


Parabolic Eq. C) Neither Applicable or Practical 





Figure 3. 1 Applicability and practicability of model types, adapted from [Etter, 1996] 


Of particular note is a paper published in the Journal of the Acoustical Society of America, 
{Smith et al., 1991], which states that due to sensitive dependence on initial conditions, ray-theory 
solutions tend to exhibit chaotic behavior. The paper goes further to show that the chaotic behavior 
begins to become apparent somewhere around 100 km from the source in deep ocean environments with 
no boundary interactions. Chaotic behavior for most shallow-water environments will occur at shorter 
ranges due to extensive ocean bottom interaction. For this reason, care must be taken in evaluating the 


results of complex tactical shallow-water environments. 


D Type Of Simulation 


Choosing the type of simulation is important to any modeling problem. To simulate a shallow 
water environment one might use a physical simulation or an abstract simulation. An example of a 


physical simulation of a shallow-water environment is to construct a tank of water with a scaled version of 
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the environment to be studied. An example of an abstract simulation is a computer simulation of a sonar 
computational model, where the actual simulation is occurring as a bunch of electrons flowing through 
computer chips. Each method of simulation has its own advantages and disadvantages. The physical 
simulation has the advantage of being closely linked to reality, providing that proper scaling and 
terminating techniques are used. However physical simulations require massive amounts of material and 
time reconfiguring for a new environment. Abstract simulations are just the opposite in that they are 
further removed from reality even when the initial and boundary conditions are specified properly, but 


they are easıly reconfigured for different environments by changing the initial and boundary conditions. 


3. Networking Considerations 


The choice here is to opt for or against network interaction. Networked simulations have a 
distinct advantage over single-computer simulations. Networked simulations can be made more widely 
available and can as a result have more computers working on a given problem. In large-scale 
simulations this increased computing power and ease of access allow the formation of multiple participant 
simulations. The advantage here is that the simulation is closer to reality and any simulation that is closer 
to the real world progression of a system is a better simulation. However, networked models must also 
pay performance penalties associated with bandwidth capacity, message delivery latency and | 
synchronization overhead costs. The decision to network or not is primarily a function of whether or not 


computing power needs to be increased or if interaction with distant parties 1s desired. 


4. Visualization 


The options for visualization range far and wide. Output options include text-based tables, two- 
dimensional graphs and three-dimensional virtual worlds. The first two options are modest and well 
within the capabilities of most computer software that is available. The third option requires specialized 
rendering software. This software comes in many forms, such as stand-alone programs that read in 


textual tables and generate information, interactive programs that are directly called from the simulation, 
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or a combination of the two as can be found in common virtual reality (VR) browsers. Each form has its 
distinct advantages and disadvantages. One recent advantage, that stands out, is provided by VR 
browsers: widespread availability of a common standard for 3D graphics representation, using the Virtual 
Reality Modeling Language (VRML 97). This compatibility supports the needs of large-scale simulations 


presented in the previous section. 


E. SUMMARY 


Implementation considerations for a shallow water sonar propagation and visualization 
simulation are many. Proper choices for propagation model, simulation type, networking and 
visualization must be made. The remainder of this thesis discusses and develops these important choices 


in detail. 
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IV. RECURSIVE RAY ACOUSTICS (RRA) DERIVATION 


A. INTRODUCTION 


This chapter is included as a benefit to the reader interested in the mathematical derivation of the 
Recursive Ray Acoustics (RRA) equations. The nomenclature used is similar to [Ziomek, 1996] and each 
equation in the stated derivation comes from Chapter 5 of that text book. The flow of the mathematical 
development presented here is straight to the point. Some of the more difficult (and long) tangential 
derivations are left for the reader to explore in Ziomek's text. Although this derivation is streamlined, it 
is a complete and rigorous treatment of the transformation of the linear acoustic wave equation into a 


linear ray acoustic equation. 


B. THE DIFFERENTIAL EQUATION 


[Ziomek, 1996] Chapter 5 is devoted to solving the linear acoustic wave equation using ray 
theory simplifications. This section derives the pertinent equations of ray acoustics. The first assumption 
is that the acoustic source is a time-harmonic oscillator. When this is the case, the starting point of the 


derivation is the Helmholtz equation: 
V^o(r)* kon*(r)g(r) 2 O. (4.1) 
The acoustic pressure is then assumed to have the simple geometrical form which locally appears as a 


plane wave front of amplitude a(r), as in: 

olr) = a(r)e 7 (4.2) 
The function W(r) defines a surface of constant phase and typically is called the eikonal, which means 
image in Greek. Upon substituting Equation (4.2) into Equation (4.1) and simplifying, the following 


equation is obtained: 


۳ > ke In? (r) - IVW(r) IZO - jk, |a(r)V?W(r) +2Valr) e VW (r)| U (43) 
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Now since the left hand side of the equation is identically zero, both the real and imaginary parts of 
Equation (4.3) must be zero. This provides two equations: 
va(r) + be) - Iva la(r) = o (4.4) 
and 
a(vr)V^W(r) 4 2Va(r) e VW(r) 2 0. (4.5) 
Equation (4.4) is the equation that tells how the wave front of the acoustic signal travels through the water 
and Equation (4.5) describes the flow of energy of the acoustic signal. Equation (4.4) is now the focus of 


the following development. 


Rearranging this equation yields an alternate form: 


سی 1 


7 222) a(r) 


From this form if it is assumed that the amplitude function of Equation (4.2) vanes much more 


1 — ki vre =0. (4.6) 


slowly than the phase function, then 





y? 
ar) «« kl vw) (4.7) 


a(r) 





is true and Equation (4.6) reduces to 
¡VW(r)]* - n'(r) , (4.8) 


the eikonal equation. As a direct result of the assumption in Equation (4.7) Ziomek shows in his book 


that 


u >> k’(r) = En (4.9) 


and also indicates that the approximation amounts to a high-frequency approximation. 





From vector mathematics it is known that the gradient of the eikonal is 
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O i 0 ۵ O 2 
= — —— es ‚10 
VW (r) ay W Ar) xt a ٦ت‎ SDE (4.10) 


Since it is a vector it can be expressed as a magnitude and a direction, 


VW(r) - |VW(r)|s(r), (4.11) 
where sr) defines the local direction of the ray path. This can be rewritten, using Equation (4.8), as 


VW(r) - n(r) sx). (4.12) 
where n(r) is the index of refraction. To further clarify Equation (4.12) a ray is defined to be a path 


perpendicular to the phase fronts at a given position, as in Figure 4.1. 





Phase 


Front 
Figure 4. 1 Ray Definition 


From this graphic it is easy to see that the following equation is true: 





)4.13( نات DE‏ رھت = ir‏ ای 
ds Os Os o‏ 


since |dr| = ds. Equation 4.13 is the unit vector in the direction of Equation (4.10). Computing the 


directional derivative of W(r) and using the chain rule yields the change in phase along the ray path, 


I ZEIT ENT + Zu). am 


From inspection of Equations (4.10) and (4.13), it can be seen that 
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Q ^ 
A W(r)- VW(r)e s(r) (4.15) 
is equivalent to (4.14). Substitution of (4.12) into this result produces 
a W(r)=n(r). (4.16) 
e 


This equation can be integrated quite easily to yield the solution of the eikonal equation. In 
essence the equation says that a small change in the wave front can be generated from the product of the 


index of refraction and a small change in the path length of a ray traveling normal to the front. 


C: THE DIFFERENCE EQUATION 


To develop this equation into a form better suited for calculation on a digital computer and to 


solve for a more meaningful physical quantity, the gradient operator 1s applied to Equation (4.16) to yield, 


C V n(r). (4.17)‏ کاب 
OS‏ 


Since the two operators preceding W(r) commute, Equation (4.17) can be rewritten as, 


fo DS V n(r ). (4.18) 
Os 


Equation (4.12) is then substituted into Equation (4.18) yielding, 


2 


—( (r)s(r)) = Vn(r), (4.19) 


which can be rewritten as a difference equation, 
A(n(r)s(r)) = Vn(e)As, (4.20) 


instead of a differential equation. If the path length step As is kept small, then the path can be considered 


a straight line, and Equation (4.20) can be rewritten again as, 


n(r,)s(r,)- n(r, )s(r, ) = Vn(rXr,-r,). (4.21) 
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Substitution of, 


Co 


c(r)' 


and solving for the normal vector to the wave front at the final position yields, 





n(r) = 


(4.22) 





s(r,) = or E ہہت‎ c(r”), (4.23) 


where c(r ) is the average of the sound speed at the final and initial points. 

This is the final form of the solution and is referred to as the recursive solution, since to get the 
final solution the previous step must be solved recursively all the way back to the starting position. The 
recursive technique is relatively simple in that, if one starts with a normal to the wave front and takes 
small-path-length steps, calculating a new normal at each step, the solution can be achieved. It is this 


difference equation that is implemented in the RRA code of this thesis. 


D. SOLUTION TO THE TRANSPORT EQUATION 


Equation (4.5) is now addressed to discover how the acoustic energy is transferred (i.e. 


transported) with the ray. The first step is to multiply the equation by a(r) which yields, 
a^ (r)V?W(r) 4 2a(r)Va(r) e VW(r) Ve [a^ (1) VW(r)| — (4.24) 
Calculating the time averaged intensity vector of, 


olt,r)= p(rje?*, (4.25) 


where p(r) is Equation (4.2), the following average intensity is arrived at; 
1 
I, (r) E 5 ko Polea) Wir) (4.26) 


[Ziomek, 1996], Chapter 5, page 336 shows the development of the time-averaged intensity vector in 


detail. When Equation (4.26) is substituted into Equation (4.24), the result is, 
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Ve 0 É | = (427) 


Now suppose that instead of a considering single ray proceeding along a path, a bundle of rays 
such as Figure 4.2 exists. When Equation (4.27) is integrated over the volume defined in the figure and 


the divergence theorem is applied, 





Figure 4. 2. (a)Rays forming the bundle, (b) Surfaces enclosing the bundle — [Ziomek, 1996] 


۲ s Em = ار‎ edS < 0 , (4.28) 


is the result. The integral over the surface S;, in Figure 4.2, is assumed to be zero since the rays 
themselves form the surface and by definition cannot leave the surface. Therefore, it 1s easily seen that the 
integral over the three surfaces result in showing that all of the energy stays in the bundle and as such the 


energy in the bundle crossing any plane is a constant. In equation form this concept is written as, 


D (r 1 )s 1 == (r, )s 2 (4.29) 
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where the S’s are the surface areas perpendicular to the ray bundle and the p’s are the density at the two 
different positions. This equation allows the calculation of the intensity of the bundle at any position 


along its path, at any time. It is important to note that no summations of incremental intensity variations 
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along the path are necessary, only the knowledge of the values of three of the four quantities allow 


calculation of the other at any two positions in the water column. 


E. SUMMARY 


As can be seen from the preceding derivation, the Recursive Ray Acoustic (RRA) approximation 
is straight forward. The particularly fascinating outcome of this entire chapter is that when any acoustic 
wave meets the assumptions made in the mathematical development, the wave can be completely 
described by just two simple Equations, (4.23) and (4.29). These correspond directly to Equation (3) in 
[Ziomek, 1993] and Equation (5.2-170) in [Ziomek,1996]. Thus the RRA formulations appear to be both 
very general and computationally efficient. These attributes make the RRA model an excellent (and 


perhaps optimal) choice as a computational engine for real-time 3D sonar beam generation 
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V. MODELING OF RECURSIVE RAY ACOUSTICS (RRA) 


A. INTRODUCTION 


This chapter extends the recursive ray acoustics (RRA) algorithm into a complete model suitable 
for real-time simulation and 3D sonar visualization. The full set of equations that define the RRA model 
are stated and elaborated upon. Each sub-model of the RRA (sonar server) model is explained from a 
modeling perspective. The composition of the state vectors and operations for each sub-model are defined 
and the mathematical developments for curvature determination, reflection detection and target detection 
Bie given. This chapter lays out the formal structure of the RRA model, which is further described in 


Chapter VII, Model Implementation and Integration. 


B. MODEL ENGINEERING 


At the outset, the difference between the term ‘model’ and the term ‘simulation’ needs to be 
defined. A model is the abstract representation of the dynamic system in question. Simulation simply 
refers to a the performance of a model over time. Implementation of the model can predict the dynamic 
characteristics of the system, from which we gain insight about the original model. It is quite clear from 
the definitions that the formulation of a model must precede the development of a simulation. Therefore, 
before one can jump into the mechanics of programming the entire model into a simulation, model 
engineering must be carefully considered before a single line of code is written. Fishwick states, 

How do we engineer models? While there are many modeling techniques for simulation, we are 

often in a quandary as to which model technique to use and under what conditions we should use 

it. First start with a concept model of whatever dynamic system 1s being investigated. Break the 

system into a hierarchy of abstractions and then choose models to represent those abstraction 

levels. The final solution to a model will be a multi-model of the system since no one model type 


will be sufficient to describe a system except in only the most elementary circumstances. 
[Fishwick, 1995] (page 5) 
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Figure 5.1 shows the top-level architecture of the shallow-water sonar model. The sonar server 
model is considered ın this chapter, the sonar visualization model is considered in Chapter VI and the 


virtual world display model is considered in Chapter VII. 










Recursive Ray 
Acoustics 
Shallow Water 
Sonar Model 










Virtual World 
Display Model 


Sonar Server Visualization 
Model Model 


Figure 5. 1 Shallow-water RRA sonar model architecture 


The sonar server model is further decomposed into the hierarchy of models depicted in Figure 5.2. Asan 
overview, the base models in this hierarchy are the target model, the ray model, the ocean bottom model, 
the ocean surface model and the sound speed profile model. From the ray model the beam model is 
derived and from the beam model the lobe model is derived. Each sub-model and its implementation in 


the simulation is discussed in detail in the following sections. Figure 5.3 shows the relationship among 


Sonar Server 
Model 
Target Model Lobe Model Bottom Model | | Surface Model| | Sound Speed 
: Profile Model 


Beam Model 
Ray Model 


Figure 5. 2 Sonar Server Class Hierarchy 


rays, beams and lobes. 
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a) b) C) 


Figure 5. 3. a) A single ray, b) Four rays forming a beam and c) Twenty beams forming a lobe 


e. THE COMPLETE RRA EQUATIONS 


Before further elaboration on sub-model design and implementation, it is necessary to revisit the 
ray theory solution to the wave equation. Equation (4.23) is the central equation in the algorithm 
intended for eventual implementation. [Ziomek, 1996] shows the final form of the equations in a finite 
difference format vice the differential format. The equations are restated with minor variations from the 
originals for ease of understanding. 


r,=r_,+Ar_, (5.1) 


Ár, ,— Ar , Sii, (5.2) 


^ ^ i 1 o 
Si = Si-1 A Are(r Jaga Vela). (5.3) 
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where the unit normal vector S is the same as Equation (4.13) and i = 1,2,3,... Since most real-world rays 
are represented as an elevation and an azimuth angle, a coordinate system conversion must be made. 
These conversions are summarized in the following equations, 


Ox 


ERE sin(8 )cos($), 65) 
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- cos(£). (5.6) 


Y 
وج‎ 


and 


2 
5 = sin(B)sin (6), (5.7) 


where f is the elevation angle and ¢ is the azimuth angle relative to the initial ray position, as shown in 


Figure 5.4. 


Z - cross range 





X - down range 
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Y - depth 


Figure 5. 4 Ziomek coordinate system definition 0۶ ۵ and d. ß 1s measured from the y-axis and ¢ 
is in the x-z plane 


Equations (5.5)-(5.6) assume the model uses the coordinate system of Figure 5.4, referred to as the 
Ziomek coordinate system to distinguish it from systems used later. Figure 5.5 shows the Ziomek 
coordinate system and the traditional naval coordinate system. Using the preceding equations and 


conventions specified it is now possible to unambiguously specify each sub-model of the RRA model. 
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Figure 5. 5 Ziomek and Naval coordinate systems 


D. RAY MODEL 


With the construction of any model, one must first work out what data is important and limit the 
model to produce just this set of information. While the ability to calculate every foreseeable parameter 
hypothetically might be desirable, the reality of having to implement a practical model to such a level of 
detail proves difficult. In general, with an excessive level of detail, the time taken to run the simulation 
and collect the data increases dramatically. Therefore, what parameters must be modeled in the case of a 
sound wave traversing an ocean environment? To adequately describe the state of the model, and how it 
transitions from one state to the next, the designer must clearly state what the desired output of the model 
is. As seen in the previous chapter, a wave front traveling through the ocean can be thought of as a 
collection of sound rays traveling through the water. It is just these rays that are considered in this 
model. 

From the ray model the information that is available at each time step includes position, duration, 
direction of propagation, phase, time, total distance traveled, and total attenuation of the ray. In an ideal 
case, all information at every step of the calculation might be recorded, but for a single ray this can add up 
to megabytes of data. Obviously, this is not practical due to the excessive disk access and computational 
overhead involved in keeping that amount of data. Instead the model only keeps the amount of data 
necessary to visualize accurately what is developing in the model. Since a ray traces out a curved line in 


the ocean with occasional abrupt changes in direction at reflections, the model needs to retain enough 
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points to faithfully represent the ray curvature, especially for neighboring points in the vicinity of an 
abrupt change in direction of the ray. Therefore, in addition to the required output information, local 
curvature and reflection information must be kept. Thus the state vector of the ray model must include the 
position of the leading edge of the ray segment, position of the trailing edge of the ray segment, direction 
that the ray segment is traveling, simulation time, phase of the ray at the wave front, total distance 
traveled, and total attenuation. With the state vector now defined, the operations on this vector must be 
defined to enable the state vector to change in time. The desired operations needed are to initialize, 
propagate, reflect, calculate attenuation, calculate total curvature, save the current state and print any 
required state of the ray. Some of these operations are self-explanatory while others (propagation, 
reflection, attenuation calculation and curvature calculations) require more discussion. The reflection 
operation is discussed as part of the surface and bottom models later in this chapter. State vector 


parameters and operations for the ray model are enumerated in Figure 5.6 


Ray Model 
Operation 


Position Set Any Parameter 
Trailing Edge Position Get Any Parameter . 


Propagation Direction Propagate 


Simulation Time Reflect 
Phase Calculate Attenuation 


Attenuation Calculate Curvature 


Curvature Save State 
Print State 





Figure 5. 6 Ray model state vector and required operations summary 


The propagation of a ray first requires implementation of Equations (5.1)-(5.4). Equations (5.5)- 
(5.7) are used to initialize the unit normal vector to the wave front, and find no further use during the 
propagation of the ray. Starting with a given state vector, propagation involves applying Equations (5.1) - 


(5.4) (in reverse order) for every given change in path length. In order to keep the simulation time 
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between rays the same, the incremental path length, Ar, in any given time step, At, is calculated as 
follows, 

Ar = At-c(r’). (5.8) 
This substitution allows for stepping in regular discrete time intervals to calculate changes in the state 
vector. 

Attenuation of sound in sea water has been extensively studied and is due to many physical 
phenomena occurring simultaneously in various parts of the ocean. Relaxation processes and scattering of 
energy from resonant micro-bubbles [Clay and Medwin, 1977] are but a few of the loss mechanisms in the 
ocean volume. This model is concerned solely with the relaxation process, which occurs when a passing 
sound wave interacts with a media in a non-adiabatic fashion. In other words a net transfer of energy 
from the sound wave to molecules in the water has occurred, which results in an attenuation of energy in 
the wave. In Chapter 7 of [Kinsler et al., 1982] sound attenuation due to relaxation processes in sea 
water are shown empirically to follow Figure 5.7. Since the value of the attenuation changes with 
temperature' and pressure, those effects must be taken into account during every time step of the 
propagation. Therefore a simple summation of the product of the change in path length and the local 
value of the attenuation coefficient must be made. The formula for calculating the attenuation coefficient 
is complex and requires excessive processing time so an approximation is made. At the frequency of most 
active sonar systems (less than 100KHz), the attenuation coefficient is fairly negligible for all temperature 
and pressure conditions, thus a representative value is chosen from the range of possible values. This 
value is the attenuation coefficient as read from Figure 5.7. A more detailed discussion of attenuation and 


the relaxation attenuation coefficient is found in [Kinsler et al., 1982]. 
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Figure 5. 7 Attenuation coefficient in sea water [Kinsler et al., 1982] 


Curvature calculation for the ray is based on storage efficiency and the aesthetic appearance of 
visualizing the ray. As stated before, it is important to minimize the amount of information one carries 
along in a simulation. In the case of visualizing sound waves in the ocean, only the number of points 
necessary to make the ray appear smooth when seen from a distance must be saved for later use. The trick 
then is to calculate the cumulative radius of curvature of the ray and when it reaches a predefined small 
value, a detailed vector of state data is saved and the cumulative curvature is reset to zero. The 
development of this simple curvature relationship is presented in the rest of the section. In Figure 5.8, the 
left-hand graphic shows two ray segments laid head to tail, as they are when the ray propagates. The 


right-hand graphic shows the angular difference between the two ray segments. 
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Figure 5. 8 Successive segments in the ray path 


These vectors represent two successive steps in the tracing of a ray. Ifthe vectors are rearranged as in the 


right-hand graphic in the figure, it is easily seen that 
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where |r| = |r1| = |r2| when the constant time step is small and there are no abrupt changes in direction of 
propagation. In Equation (5.3), the values of c(ri), c(r;) and c(r*) are assumed to be equal, provided that 
At is small, therefore they disappear from Equation (5.3). Then Equations (5.3) and (5.8) are substituted 


into Equation 5.92, yielding 


A@ = isi ArVe - sı - |AIVC . (5.9b) 








It is also apparent that if the second normal does not change much from the first normal as the ray 
propagates, and the path length of each ray is small and nearly equal, then the distance from the mid- 
point of the first vector to the second is equal to the arc length of one of the rays, 

As = At-c(r). (5.10) 
Curvature then can be defined as the amount of angular change per amount of change in distance traveled. 


Thus Equations (5.9b) and (5.10) can be combined to arrive at 


)5.11( کو رت [curvature]‏ 
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which shows that for small changes in ray position with small changes in direction, the increment in 
curvature can be defined completely in terms of the sound speed profile of the medium. Since it is already 
necessary to take small step sizes and there are no abrupt changes in sound speed in the water column, 
this approximation is satisfactory. The next task is to then find the limit to which the cumulative 
curvature is to be taken before the ray state vector 1s recorded. This tends to be a matter of preference as 
to how smooth one wants the curve and how much overhead is desired. Evaluating visual renderings 
through a series of empirical measurements, a value of .004 radians/meter gave good results of curvature 
versus overhead. Now that the state vector and the processes to change the state vector in time are known, 


the ray tracing model is ready for simulation. This simulation is described in Chapter VII. 


E. BEAM MODEL 


The beam model is similar to the ray model in that it is just one level of abstraction away from 
the specifics of how rays are calculated. In other words, the beam model composes multiple ray models. 
The state vector for a beam is that of a composition of four matched rays which form a quadrilateral ray 
bundle. The only fields added to the state vector that have not been previously discussed are transmission 
loss, detect time and echo level. The transmission loss part of the state vector is used to define the 
logarithmic change in the intensity of the energy contained within the ray bundle, while the detect time 
and echo level parameters are used to document the time and intensity of a target detection. These new 
fields and operations are listed in Figure 5.9. 

The methods needed to change the beam state vector over time and to produce usable output are 
thus extensions of some of the methods developed for the individual rays. These methods are 
initialization, propagation, calculation of transmission loss and echo level, and production of usable 


output. 
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Beam Model 


Ray1 Set Any Parameter 
Ray2 Get Any Parameter 


Ray3 Calculate a Beam 
Ray4 Calculate TL 
Calculate Echo Level 
Detect Time Save State 

Echo Level Print State 


Transmission Loss (TL) 





Figure 5. 9 Beam model state vector and required operations summary 


Initialization and propagation are just extensions of the ray model initialization and propagation 
methods; these methods perform the initialization and propagation methods on each of the four rays in 
the beam state vector. The visualization methods are used to generate abstract visualizations of the data 
stored in the beam's state vector. The explanation of the visualization model is provided in Chapter VI. 
The only methods remaining are those of calculating the transmission loss of the bundle of rays and the 
echo level of the target. Since the echo level is determined from the transmission loss, the transmission 
loss is discussed first. 

In order to calculate the transmission loss of the ray bundle, it is essential that the ray segments 
are synchronized in time. As shown in Equation (4.29), given knowledge of the initial intensity, initial 
area and final area, the value of the intensity of the acoustic signal in the ray bundle can be calculated. 


This allows direct calculation of the transmission loss which is, 


IL- | | = lolog S] ۱ )5.12( 
avg \*1 


where the areas, S; and S;, are the same as in Equation (4.29). The values of the initial intensity and the 
initial area can be known from the sound source, while the value of the final area must be calculated for 


each spatial position of interest. 
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The method of calculating an approximation to the area of the ray bundle is to divide the 
quadrilateral beam front into two triangles. The area of these two triangles are then added together as the 
total area. The question is, how is an area calculated when all that is known is the positions of the four 
corners of the area? The straightforward solution is to use vector algebra to find the area. Figure 5.10 a) 
shows a quadrilateral formed by four points in space. These points are not necessarily coplanar, so they 
are divided into two triangles that are guaranteed to be coplanar. Figure 5.10 b) and c) show the two 
triangles as being transformed into parallelograms. Vector V1 is crossed into vector V2, which gives the 
vector whose length is equal in magnitude to the area of a parallelogram defined by V1 and V2. This 
value is then halved to form the area of triangle T1. The area of triangle T2 is found in the same manner 
as that of T1 using vectors V3 and V4. It is noted that the area of a parallelogram being equal to the 
magnitude of the cross product of the vectors defined by two adjacent sides is a fundamental property of 
the parallelogram. 


V4 
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Figure 5. 10 Division of beam front into triangles for area calculation. a) beam front b) triangle 
T1 transformed into a parallelogram c) triangle T2 transformed into a 
parallelogram 


The echo level, EL, is defined as 


EL = SL -2TL + TS, (5.13) 
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where SL is the source level and TS is the target strength. For this model, the users calling program 
keeps track of the SL. Thus, the model calculates EL-SL. The target strength as found in Figure 15.51 


from [Kinsler et al., 1982] is, 
= olo) | (5.14) 
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Where the symbol ‘o’ is defined to be the acoustic cross-section of the target. 


Again it is shown that the beam model propagation of sound in the water column is defined by 
the RRA model plus the addition of three simple equations, Equation (5.12), (5.13) and (5.14). Here the 
power of ray tracing begins to show, through the simplicity and ease of calculation of pertinent 
parameters. Simple calculations provide both reduced simulation time and the ability to propagate more 


rays in real time. 


F. LOBE MODEL 


This model is the simplest by far. The lobe model merely composes several beams. A lobe is 
defined as a horizontal and vertical beam width and a horizontal and vertical count of the beams in each 
width. Thus the state vector of a lobe consists of a matrix of beams. This matrix of beams determines the 
accuracy of the simulation. Since ray tracing is a space-time accurate algorithm, the number of beams 
calculated has no bearing on position and time accuracy. When the number of beams is increased, the 
size of the area of the tube Dif dies is kept small and proper interference patterns can be maintained. 
Minimization of beam front area is also important for reflection purposes. As the beam gets larger the 
individual beams that comprise it get farther apart. The beams can, if large enough, strike a surface at 
significantly different times or can even strike completely different surfaces. The current reflection and 
detection models do not account for the size of the beam fronts and thus drastic errors can be introduced 
into the simulation if beam front sizes are not maintained small. The operations that can be performed on 


a lobe are initialization, resetting, calculating and printing. All of these operations are self evident in that 
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they will rely upon the same type of operations in the beam sub-model to perform their operations. Each 
operation will call the sub-model operation for each beam in the lobe. The initialization and resetting 
operations also perform an energy distribution function. Each beam receives a fraction of the energy of 
the lobe in accordance to the position of the beam in the lobe. In this model all beams receive the exact 
same energy, but this can be changed to any energy distribution pattern desired. Figure 5.11 shows the 


fields and operations for the lobe model. 


Lobe Model 


m by n matrix of beams Set Any Parameter 


Beam widths Get Any Parameter 
Calculate a Lobe 
Print State 





Figure 5. 11 Lobe model state vector and required operations summary 


G. SURFACE MODEL 


The surface model defines the ocean surface with which the rays interact. The current model specifies the 
ocean surface as a smooth pressure-release surface. For most cases this 1s sufficiently accurate. However, 
if surface effects are being studied, the surface model can be easily changed and re-implemented, since it 
has been constructed as a separate entity from the beginning. The state vector of the surface model 1s that 
of a variable that indicates whether or not the ray in question has intersected the surface. There are two 
methods associated with this state vector and one utility method that outputs a representation of the 
surface itself. This final method will not be considered in this discussion since it is self explanatory. The 
first method is deciding when an intersection of the ray with the surface has occurred. The second method 
is the actual performance of the reflection. The reason for thinking of this as two methods instead of one 
is simply semantics. The first performs a query into the current state of the system and the second 
actually makes changes to the state. This separation was provided so that if the user model referencing 
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the surface model needed to make this distinction, the option was available. The reflection query 
operation is simply performed using vector algebra. Figure 5.12 provides a graphical representation of 
how a surface interaction is determined to have occurred and how a reflection is made, with arrows 


indicating the reflected ray path. 


P 


previous 







Plane of Reflection 


em 


reflected 
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Figure 5. 12 Reflection of a ray from a plane 


From the ray model it is known that the current state of a ray can be partially defined by a 
position and a direction normal that indicates what direction the point is moving. It is also known from 
vector algebra that a plane surface can be represented as a point and a normal to the surface. From this 
information a detection can be determined. First the vector from the defining point of the plane to the 
current ray position is constructed as, 


S-P. -P 


ray plane * 


(5.15) 


Then the perpendicular distance of the ray from the plane is calculated using, 
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N 


0 < -5 )5.16( 
This distance then indicates whether or not the ray has penetrated the bottom. If the value of d is negative 
an interaction has not occurred and if d is positive the ray has penetrated into the bottom. 
The reflection problem is founded in Snell's law, so it is there that the derivation must start. 


Snell's law is stated here as, 


sin(4)  sin(.) (5.17) 


— . 


C, C 


1 r 





Notice, however, Snell's law is stated in an inherently two-dimensional scalar form and that the ray model 
is stated in a three-dimensional vector form. The question is then applying Snell's law to ray model 
vector form. It is important to note that the law was derived for three-dimensional space and this two 
dimensional-simplification is valid for the reflection since it occurs in the plane formed by the incoming 
ray and the normal to the surface. Figure 5.12 shows a two-dimensional representation of the three- 
dimensional problem. Since Ci equals Cr, both the incident and reflection angles in Equation (5.17) are 
equal as well. From the figure it can be seen that the only correction to the ray needed is in the depicted 
plane in the direction of the normal to the reflecting surface. This reflection position can be found by 
adding twice the perpendicular distance of the ray position from the surface, to the ray position, which 


results in the reflected ray position. Mathematically this is shown by, 


A 


P enectea = Pray + 2d N piane . (5.18) 


r 


This forces the reflected normal to change in a like manner as shown, 
N reflected = Dray — a(n ray ® 9 Il plane - (5.19) 


These few equations govern all reflections in the surface model and enable precise calculation of reflected 


rays at the ocean surface. 
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H. BOTTOM MODEL 


The ocean bottom model is similar to the ocean surface model. The major difference is that the 
bottom is modeled as a rigid surface instead of a pressure-release surface. Thus the only difference is in 
the phase of the reflection from the bottom surface. Instead of a reflection phase change of x, the 
reflection phase change is zero. Also the bottom has two contour configurations. The "no slope" contour 
is much the same as the surface, but a depth needs to be specified, unlike the surface model. The “sloped” 
contour is a little more complicated; it has a flat deep bottom, a steeply sloping wall and a gently sloping 
shelf that tapers off to a depth of zero. All of the parameters of the “sloped” contour are changeable 
which makes for a highly versatile model. Figure 5.13 shows a three-dimensional representation of the 
surface and the “sloped” bottom models. The bottom model uses the exact same state vector and 
operations as the surface model. Therefore, it also uses the exact same equations as the surface model and 
as such a discussion of the specifics of a bottom model is not warranted. 

A final statement as to the generality of the bottom model needs to be made. Any surface can be 
broken down into a near plane surface by subdividing the surface into infinitely small pieces. The bottom 
model can be broken into a piece-wise planar (gridded) surface, with each plane specified by a point on 
the plane and a normal to the plane as well. With careful planning, almost any bottom interface can be 
modeled with ease and low computational overhead. Care must be taken however to ensure that the 
characteristic size of the bottom structure does not become so small that it approaches the scale of the 
sonar wavelength. This is because errors due to phase interactions become large, in this model, when 
bottom structure and wavelength are of the same order. An important area of future work is the addition 


of methods which can import standard bathymetric databases and compute corresponding normal vectors. 
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Figure 5. 13 Example sloped bottom contour with surface model 


I. SOUND SPEED PROFILE (SSP) MODEL 


This model is a simple one. It has no state vector and thus supplies information about the current 
time and position only. Its purpose is to return the sound speed and the rate of change of the sound speed 
with position. A detailed discussion of SSP is available in [Urick, 1975]. It has five basic profiles that 
come from [Ziomek, 1996] as the standard profiles. Again due to the highly modular breakdown of the 
RRA model, changes in the sound speed profile model are made easily. Chapter VIII Simulation Results, 
details in specific the SSP models used in the current simulation. An important area of future work is the 
development of an SSP database linked to the bottom model, import of real-time SSP data and SSP data 
visualization utility operations. It is important to note that the broad generality and accuracy of RRA is 


directly dependent on the quality and accuracy of local bottom, surface and SSP models. 


J. TARGET MODEL 


The target model is designed to keep track of all targets in the virtual world and to provide 
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operations that determine if a target has been detected by a sonar pulse. A target is defined as any object 
that is not environmental, geological or biological. Therefore the target model gives us information about 
the man-made objects in the virtual world. It does so by keeping information about each target in the 
virtual world and comparing that to position information provided by the beam model. The state vector of 
the targets contains all the information about the targets. 

In the current implementation the targets are randomly placed moored mines and submarines 
that have no movement. Thus in the current target model the target state vector does not change. This 
state vector is composed of the position and radial size of the mine. As for the submarine, it is considered 
to be composed of several mines and as such just generates a chain of mines the size of the submarine. 
This model is intentionally simplistic and much further work is warranted in target modeling. 
Additionally, information pertaining to the cumulative area of targets detected in the current time step in 
the beam is saved. The operations affecting the state vector are initializing, resetting, detecting a 
collision, getting the area of detection and outputting the target information. Again initializing, resetting, 
getting detection area and outputting are self explanatory. 

The only operation left is computationally efficient collision detection. This model uses the same 
equation for detection of a possible interaction between the sonar beam and the targets as the surface 
model uses to detect interactions between the ray model and the surface. Equations (5.15) and (5.16) are 
used to indicate when the plane of the beam front passes through the target. Therefore the perpendicular 
distance from the plane of the target, d, is used to tell if the target has the possibility of being detected. 
When the sign of d changes from positive to negative this means that the wave front has encountered the 
target in the current time step. Since the plane passing approximating the wave front of the beam is 
infinite in extent further operations are necessary to check ıf the target (or any portion of it) is in the beam 
tube. This detection is performed by checking to see if the target lies within the X, Y and Z coordinates of 
each pair of points on the box that make up a diagonal in the box. Figure 5.14 visually depicts the 


bounding box concept. The presence or absence of a target 1s reported at the end of every check. 
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Determining the distance to every object is also useful in another way. The distance of the closest target 
during each check is recorded and this distance is divided by a nominal ocean sound speed of 1500 meters 
per second. This time is then divided by the time per propagation step. The result is the number of 
iterations that can occur without checking to see if a collision occurs. Essentially, this number of time 
steps that can be skipped in between checks for target interactions is vet another member of the state 


vector. This has the potential to be a dramatic speed enhancement over checking all targets every time 
step. 
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Figure 5. 14 Bounding box target detection. A target in the sonar ping volume is considered detected 


K. PINGSERVER AND PINGER MODELS 


The pingserver and pinger models are included as an example to show how the RRA class library 
can be used. The pinger model is a stand-alone sonar ping model that produces text files that are either 
raw data or three-dimensional renderings of the sonar lobe. The pingserver model is where the true power 
of the RRA model comes into play. The pingserver uses a client-server communications approach, where 
the client asks for a sonar beam to be projected into the virtual world and then the server calculates the 
necessary data and returns this data back to the client. A major benefit of this approach is that a proxy 
server can be placed in between the client and server to improve overall system scalability. This proxy 
server appears to the server as a client and to the client as the server (hence the name). Since the RRA 
algorithm is a space-time based algorithm (no transforms to other domains), the propagation calculations 
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can easily be partitioned and shared among multiple pingserver implementations. Thus, the proxy server 
has the job of getting the request for the sonar ping and then distributing the calculation of the many 
beams required to the many available servers. With the advent of large-scale network clustering, the 
possibility of distributing real-time RRA calculations is highly attractive. 

At this point, a step back from modeling to discuss computer implementation seems appropriate. 
Network clusters, [Hill, et al., January 1998], are collections of large numbers of moderately fast general 
purpose computers connected via a fast local-area network (LAN). The motivating idea is that when a 
highly parallel calculation needs to be done, the task can be split up among the many computers. Such 
clustering can provide gigaflop performance from machines that provide only megaflop performance 
individually. The most amazing aspect of this clustering is that since the individual computers are mass 
produced and general purpose, each one has a very low price. The cost for the creation of one of the 
clusters, Loki, was around $50,000 for a peak performance of 1.4 gigaflops. This cluster was initially 
made around 1996 and production of this same system as of 1998 1s estimated to cost 50% less. It is seen 
that the RRA model lends itself to this type of massively parallel computing and it is the pingserver model 
that provides the example of how easy this parallel programming can be. Thus, the prospect of high- 
resolution sonar beam calculation in real time appears quite feasible. Distributing RRA sonar beam 


calculations over many computers is an important area for future work. 


L. SUMMARY 


This chapter has presented the basic model for the implementation of a RRA simulation. The use 
of the models presented has provided the ability to simulate a virtual sonar environment. The ray model 
provides the basic structures to propagate a sonar ray in real-time through a virtual world. The beam 
model combines four rays to form a beam tube and then adds operations for attenuation and collision 
detection. The lobe model combines a variable number of beams and propagates all of them as one unit. 
The ocean bottom and surface models provide reflection information to the ray model. The SSP model 
provides information on sound speed to the ray and target models. The target model maintains a target 
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database and provides efficient operations to the beam model to determine ıf a detection has occurred. 
The pingserver and pinger models combine the various models into a cohesive simulation and provide 
example implementations. 

The next step in the process is to convert these model ideas into a working simulation. Chapter 
VII addresses the topic of implementation and integration of the RRA model with the visualization model. 


Appendix A contains the simulation code for the RRA models. 
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VI. SONAR VISUALIZATION 


A. INTRODUCTION 


Simulations typically require some form of visualization. Whether the visualization is in the 
form of text tables or graphic plots, the simulation must produce data in a logical and orderly grouping of 
similar information. In the early days of scientific discovery, information was collected in tables and 
occasionally plotted on paper, usually by hand. With the coming of the digital computer, the mainstay of 
data visualization shifted from text tables to graphic plots, usually in two dimensions (2D). With the 
increase ın computational power of digital computers and the widespread availability of three-dimensional 
(3D) rendering software, the time has come to shift once again to a higher form of visualization. That 
higher form of visualization is 3D rendering. 

In the fleet today, sonar tactical planning is still performed from text tables and 2D plots. This 
approach was fine for twenty years ago when less data was processed and evaluated, but today the sheer 
volume of data and the finer resolution of information obtained demand a more robust visualization of the 
incoming data. Since we live in a 3D world, the next step ought to be 3D plots. [Karahalios, 1991] 
discusses the history of data visualization. In that thesis she quotes from [Defanti, 1987] stating that half 
of the human neo-cortex is devoted to visual information processing. Thus, since our brain is wired for 
visual data input and that typically is 3D in nature, the natural conclusion is again 3D representations. 
This chapter explores the possible 3D visualization model for presentation of the results of the RRA 


simulation. 
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B. VISUALIZATION CONSIDERATIONS 


1. Tactical Visualization 


In the U.S. Navy fleet today, visualization of four dimensional space-time is performed using two 
dimensional plots. These plots have various combinations of positional and time information as the axes 
of the plots. Time-range, time-bearing, and geographic plots are but a few of the examples in widespread 
use. In addition to these time-space plots there are also time-frequency and space-frequency plots, again 
in two dimensions. The usual practice for evaluating these plots in real time is to show several varieties of 
them at the same time on separate computer displays. There are typically several console operators and 
one supervisor interpreting such plots. The job of the supervisor is to integrate the plots mentally and 
develop a four-dimensional (space-time) mental model of the current sonar-based tactical picture in the 
real world. 

In years past when shipping was not nearly as heavy and detection ranges were much shorter, the 
incoming volume of sonar information was not overwhelming. Today, with increased shipping and 
increased detection ranges, the volume of information processed by a ship's sonar system has dramatically 
increased. This has often led to an overload of information on the operators themselves. To complicate 
matters, the quieting of modern submarines has led to decreased detection ranges for the primary threat of 
interest. With the increase in noise and decrease in signal, these shortened detection ranges result in 
shortened reaction time. Though the human brain is wonderfully created, expecting an information- 
overloaded supervisor to make snap decisions while mentally integrating several plots into one multiple- 
dimension plot is extremely difficult. 

While Gary Kasparov proved that a man using intuition could frequently outthink Deep Blue, 
IBM's chess-playing supercomputer, one must remember that a chess board is infinitely less complex than 
the real world. Deep Blue on the other hand proved that raw processing power and proper programming 


eventually can defeat one of the best chess-playing minds. This digression points to the fact that the 
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human brain has limitations in its ability to process large volumes of information and infer likely 
outcomes. This fact leads one to believe that more innovative means of visualization are necessary to 
allow the operator and supervisor to step one more level of abstraction away from the raw data, so that 
timely and accurate decisions can be made. 

In the scientific research community, sonar visualization has really not progressed much past the 
state of the art in the Naval Service. Typically, the extent of visualization is colorful 2D plots, which in 
essence are pseudo 3D plots. Nevertheless, when complex space-time processes are being investigated, 
ingenious use of 3D plots, color, intensity, and transparency can provide the researcher with pseudo four 
and five-dimensional plots. All of these properties can be displayed using existing rendering programs. 
As with most physics solutions, judicious simplification and abstraction are a must if a solution is to be 
obtained. Scientific visualization of sonar data can benefit from careful simplification and abstraction as 


well. 


D. High-Dimensional Space 


Acoustic data is typically a volumetric data set rather than a surface data set. In addition sonar 
has various parameters of interest that one can consider as additional dimensions. Sound intensity level, 
time of travel, and frequency are but a few examples of these extra dimensions. In the visualization | 
discipline this is termed a “high-dimensional data space.” For the dimensions just mentioned a six- 
dimensional plot needs to be generated to simultaneously render these parameters. How can such a plot 
be conceived? In the world of 2D plotting, special techniques of coloring and use of contours transform 
these plots into pseudo three-dimensional plots. Similar techniques can be applied in 3D renderings. By 
mapping the additional dimensions to other graphics quantities besides spatial, pseudo high-dimensional 


space plots can be generated. 


3. Available Graphics Parameters 


The question then is how to map these extra-dimensional parameters. In most 3D rendering 
languages, the objects that are created have rendering properties themselves that can be used as 
dimensional extensions, much the same as contour lines are used on geographic mappings. The 
properties that objects possess include color, intensity, transparency, smoothness, reflectivity and 
emissivity. The natural tendency at this point is to imagine that with all of these parameters a pseudo 
nine-dimensional space is possible. This is just not so. Many of the parameters rendered at the same time 
produce unpredictable results or do not elucidate the information intended. For example, no matter what 
color is given to an object, if the object is transparent, it cannot be seen. So indeed, when one decides to 
implement high dimensional space data rendering, careful attention must be given to ensuring that the 
parameters employed do not overlap in such a way as to result in incomprehensible data. These 


fundamental challenges form the basis of scientific visualization research. 


4. Mapping Data to Graphics Parameters 


The task of mapping, for the reasons stated in the last section, must be undertaken with great 
care. Suppose that time of travel of a sonar pulse was linked to the transparency parameter and that the 
intensity of the sonar pulse was mapped to the color. If care is not taken, long before the sonar pulse 
becomes insignificant, as far as intensity is concerned, the pulse may not be visible in the rendering. This 
premature invisibility produces a loss of information in the scene that cannot be recovered. Now suppose 
that the two rendering dimensions and parameters are swapped. Thus, when the sonar pulse intensity 
becomes negligible the time (mapped to color) will no longer be visible. This situation is of little concern 
since the information on intensity is of greater importance that the information on time. In other words, if 
there is no perceptible sonar pulse, the time that it arrived does not matter. This example indicates that a 
hierarchy of informational importance and a hierarchy of object parameters, likely exists, and is variably 


dependent on the goals of the user. The hierarchy of parameters must be chosen so that the parameters on 


56 


the lower end of the hierarchy have little influence on the parameters that are on the higher end. Careful 
use of such a hierarchy may prevent the rendering engine from eliminating important data from the 


visualization. 


S: No Single Right Answer 


Numerous scientific visualization research projects are presently being conducted that are 
investigating high-dimensional data space and the techniques used to make sense of the vast quantities of 
data contained therein. A particularly interesting finding is that there is no single right answer as to how 
the data must be viewed. The type of data and how it is to be used directly influence the way in which the 
data is presented. The format for presenting information from which snap decisions are made will likely 
prove to be different from the format from which keen insights into physical processes are made. Indeed, 
it is often user interaction with the visualization data that provides insight rather than any single 
presentation method. Thus, the most important goal and contribution of this project is support for user 
exploration of high-dimensional sonar information using real-time 3D graphics and scientific 


visualization techniques. 


6. Initial Visualization Recommendations 


The initial recommendation for visualizing sonar data is that the three spatial components are 
mapped to the three spatial coordinates of the visualization system. The sound intensity is mapped to 
transparency in a dynamic visualization, while color is mapped to detectability. In a static beam, 
intensity, distance from source and detectability are mapped to color. All three can be mapped to color by 
providing a means of selecting the desired parameter. A full investigation of sonar parameters and the 
means of best visualizing them definitely needs to be continued as future work. These initial visualization 
recommendations are provided as a starting point for evaluation of the most important sonar parameters 


provided by the RRA model. 
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C. RRA SONAR VISUALIZATION MODEL 


1. Interactive Server Model 


Figure 6.1 depicts the standard visualization model. This hierarchy provides the general high- 
level concept of what modules comprise a visualization model. Interactive visualization also suggests the 


need for an interactive user interface model, which is depicted in Figure 6.2. 








Visualization 
Model 





User Interface | | Data Source Rendering 


Engine 


Figure 6.1 Visualization conceptual hierarchy 


a. Data Flow 


The data flow in the interactive model is quite simple. After all communications 
between the sections have been established, which happens transparently to the user, the simulation 1s | 
waiting for user interaction. Figure 6.2 shows the hierarchy of the user interface model. The user enters 
the specifics of the characteristics of the sonar pulse desired and the desired parameters of the platform 
containing the virtual sonar system. When all initialization data is entered, the user can execute the 
ordered platform and sonar pulse commands by pressing the appropriate input buttons. This returns 
control from the control panel to the battle ee manager which relays the sonar ping requests to the 
sonar server and executes the platform commands. The battle scene manager then waits for return of the 
requested sonar information from the sonar server. When the requested data arrives, it is processed and 
forwarded to the rendering system via the rendering system interface. For responsiveness threads are used 


extensively in the processing of the interactive model. 
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User Interface 
Model 


Control Request Data Receiver Rendering 
Panel Sender Interface 


Figure 6.2 General user interface model hierarchy 





Threads are processes that can be run in parallel and are used when multiple concurrent 
processes are required. Thus the platform controller has a timer thread that updates the vehicle position 
every tenth of a second. Each of the interfaces to the battle scene manager, request sender, data receiver, 
rendering interface and the control panel, are also threads. This allows the battle scene manager to 
execute each phase of the management process in parallel and only control the timing of when the results 


of each interface are added to the current state of the system. 


b. Example Execution 


This section describes how the individual components in Figure 6.1 are executed. The 
interactive sonar server is executed on a dedicated remote computer, which allows the server to be located 
on the most computationally capable computer. The proxy server, a bridge between the user interface 
client and the sonar server, is typically executed on the client computer where the user determines the 
visualization is to be run. This is done because most web browsers have security features that prevent 
them from communicating directly with other networked computers. Finally, the local web browser is 
started and the virtual world is loaded. The loading of the virtual world performs two functions; the 
rendering system is loaded and the user interface, which is linked to the virtual world, is started in the 
browser. At this point the entire interactive simulation is initialized and the control panel is ready to 
accept user input. 

As an example of simulated sonar operation suppose a single sonar ping with a three 


degree horizontal and vertical lobe width is desired. Each lobe is divided into nine sub-lobes called 
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beams. Thus each beam has a one degree horizontal and vertical beam width. Figure 6.3 shows a control 
panel with all the choices filled in and thus defines the requested ping. When the ‘ping’ button is pressed 
the control panel thread passes the pulse specifications to methods in the battle scene manager. The battle 
scene manager invokes the request sender thread and the data receiver thread. The ping sender thread 
communicates with the proxy server which in turn communicates with the sonar server that is running on 
a separate machine. The sonar server calculates the positions of the sonar beams over time and detects 
any collisions with objects 1n the virtual world. While this computation is occurring, the data receiver 
thread is waiting for information and the ping sender thread finishes execution and terminates. When 
RRA calculations are finished, information on the sonar beam positions and target detections is delivered 
over the local area network to the proxy server, which forwards the information to the data receiver 
thread. The data receiver thread returns the information to the battle scene manager and then terminates 
execution. The battle scene server then executes a VRML generation thread that converts the raw data 
into a valid VRML node structure. This node is then sent to the VRML scene. When the node is 
processed the target information appears on the screen as well as a representation of the beam propagation 
through the water. This process can be repeated over and over again to explore the entire virtual world. 
Additionally, if the battle scene manager is modified to eliminate user input and an artificial intelligence 


module is added to control the sonar and platform, an autonomous virtual vehicle can be created. 
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Figure 6.3 Interactive server control panel 


2: Stand-alone Server 


Figure 6.1 also serves as a depiction of the stand-alone server model. As with any useful 


program it must accept input, access data and return information to the user. 


a. Data Flow 

Data flow in this model is extremely simple. Appendix B shows source code for a 
typical example. In this example the virtual world is constructed and filled with targets. The platform 
position and sonar characteristics are established in the lobe class. The lobe class then constructs the 
beams that comprise it. The beam class then constructs the rays that form it and the simulation is ready to 
run. The lobe position (and thus the beam and ray positions) is calculated over time and the subsequent 
positional data 1s returned as valid VRML 97 code. Then the VRML code representing the surface, 


bottom and lobe is printed to the console. 
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b. Example Execution 


The only action that the user has to perform once the stand-alone model is coded and 
compiled ıs to execute the program. When the program is executed, the surface, bottom, sound speed 
profile, target and lobe model implementations are initialized. The lobe parameters (and thus the beam 
and ray parameters) are initialized and the lobe is calculated by calculating the time behavior of each 
beam in the lobe. Each beam is then calculated by calculating how each ray interacts with the surface, 
bottom, water column and the targets in the virtual world. Upon completion of the calculations, methods 
in each model class are called that produce VRML 97 compliant code that allow visualization of the most 
recent lobe calculation. The simple examples provided in this thesis produce compelling results. Much 
more sophisticated applications and exciting research in real-time 3D sonar visualization now appears 


possible. 


D. CONCLUSIONS 


With current sonar data collection and processing abilities continuing to accelerate, something 
needs to be done with visualization to ease the information burden on the person evaluating the incoming 
data. Creative and ingenious ways of displaying data in three dimensions must be devised if even higher 
information density is to be processed by a human. This new level of abstraction will provide enhanced 
understanding of the problem under consideration and enhanced ability to make hoe rapid decisions. A 
series of example programs which use the RRA API are presented. Initial results provide useful insight 
and are extremely encouraging. Real-time 3D sonar visualization is an important new field which now 
appears computationally feasible. A great deal of research and experimentation is needed to capitalize on 


this new tool. 
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Vil. MODEL IMPLEMENTATION AND INTEGRATION 


A. INTRODUCTION 


This chapter discusses the practical factors involved in implementing a simulation. Choosing the 
computer language with which the model 1s implemented into a simulation is a critical task. The speed of 
execution of the code, the ease of integration with some form of visualization system and the network 
interface methods are discussed as they pertain to the RRA simulation. Implementation of the sonar 
server model (RRA model) and the visualization model is summarized and the integration of these two 


models 1s discussed in the last section. 


B. CHOICE OF PROGRAMMING LANGUAGE 


Using the models specified in the preceding chapters, it is time now to turn to implementation of 
the simulation. As with most simulations in the recent past, this simulation is implemented on a general 
purpose digital computer. The primary question is choosing the best computer language to use to 
accomplish the goal. One goal is that the simulation model components might be distributed as widely as 
possible. This is a daunting task, since any two computers are rarely configured compatibly. The design 
criteria for this thesis are provision for cross-platform compatibility, computational speed, real-time 3D 


visualization and network programming. 


1. Cross-platform Compatibility 


With the creation of the Internet, cross-platform compatibility has become a rallying cry. In 
response many private sector companies have developed tools with the network and inter-connectivity in 
mind. Due to the nature of the RRA model and its multiple potential uses, choosing a programming, 


language designed for network programming is a major requirement. In the next three sections, two 
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complimentary programming languages and two networking protocols are discussed that have recently 


come into existence that allow all goals to be met. 


2: Number-Crunching Speed 


As with any inherently physical simulation of a natural process, computational speed is essential 
if real-time interaction is required. Even if real-time interaction is not required, computational speed is 
still an important consideration. However, the source of speed is no longer restricted to highly optimized, 
near machine language, already-compiled code running on supercomputers. The source of speed is the 
networked interconnection of many computers. In a sense, the network is the computer. With this in 
mind, Java, developed by Sun Microsystems, is a sensible choice for a programming language. Even 
though compiled Java code is 5-15% slower on a given single computer compared to optimized 
FORTRAN or C/C++, Java from the start has been tightly integrated with the Internet and network 
computing. This tight coupling leverages the power of the Java language to be faster in completing a 
distributed task than either of the two standard scientific languages. Although C/C++ has seen good 
success in ا‎ computing, networking in C/C++ is not standardized completely and not designed into 
the language from the very beginning. Another factor weighing heavily is that the Java development 
environment is completely free as compared to most full-featured FORTRAN and C/C++ compilers. The 
free aspect of Java allows a person developing an application to know that money will not stop a potential 
user from either not using it or from trying to use a different version of a given type of compiler that is 
already owned. When a potential user wants to modify a RRA application, he must merely retrieve the 
Java runtime environment via the Internet, install it and run the application. Since an unlimited speed 
potential and widespread availability of the common operating environment exist in the form of Java, it is 


the logical choice for implementation. 
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3. Visualization 


Though Java does provide some facilities for data visualization, the goal for this simulation was 
three-dimensional (3D) visualization. Java does not yet supply all the needed constructs to form a truly 
3D virtual world and is not used for scientific visualization (although the Java 3D libraries have recently 
become available). Once again another freely available software language came to the forefront. This 
product is the Virtual Reality Modeling Language (VRML), which typically runs and is rendered inside a 
web browser. It is usually a plug-in program to the standard web browsers provided by Microsoft and 
Netscape. Again, this combination of VRML rendering engine and web browser comes at the right price: 
no cost. VRML gives all the constructs to form an immersive 3D environment as well as constructs to 
make the environment change over time with user input. As an added bonus, due to the wide acceptance 
of Java as the language of the web, direct interface from Java to the VRML web browser is provided. 
Other graphics-rendering programs exist and are free to distribute as well, such as openGL, but none are 
so tightly interfaced, as widely available or as popular as VRML. This browser interface allows the 
interconnection of standalone Java programs to VRML scenes through the web browser. Clearly the 
choice of visualization language is VRML since it integrates so tightly with the Java compiler already 
chosen and can immediately render RRA results on any personal computer. More information on 


integration of Java and VRML is found in [Brutzman, 1998]. 


4. Networking 


Java already has many fine networking features. Typically the data that can be passed through 
network links is formatted in a way that only reader/writer programs written together can decode the 
information. In highly structured simulations, entire state vectors need to be passed. Typically state 
vector changes render all previous versions of à program useless. Two protocols are available which 
address standardization of data passing. These protocols are Common Object Request Broker API 


(CORBA) and Distributed Information Simulation (DIS) protocol. These protocols are robustly designed 
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and freely available, but also are complex and difficult to learn. Another candidate is the dial-a-behavior 
protocol work 1۳ progress at NPS which simplifies formal specification of over-the-wire data formats 
[Brutzman, August 1997]. Thus for the purposes of this thesis the networking facilities provided by Java 
proved to be sufficient. In the future, as this simulation is integrated into a larger distributed simulation, 
integrating one of these data-oriented protocols will be useful and necessary. More information on 


graphics internetworking is found in [Brutzman, August 1997]. 


C. MODEL IMPLEMENTATION 


Now that the necessary software tools have been chosen to transform the many models into a 
simulation, it is time to present some aspects of the implementation that shaped the final Java class 
library for the RRA simulation. The following two sections will discuss some of the particulars of the 


implementation of the RRA and visualization models into a full-fledged API. 


1. RRA Implementation 


Implementation of the RRA model was straightforward since Java is an object-oriented 
programming language. Each sub-model of the RRA model was developed as a separate class, since 
classes in Java are an encapsulation of variables and methods in one common object. Member variables 
contained in the object keep track of the state of the object and methods contained in the object allow 
changing of or accessing of the data in the class. This correlates well to the model concept of state 
vectors and operations. À completed class then becomes a building block for other higher-level classes. 
As an example, suppose the ray model has been implemented as the ray class. This class itself becomes 
the definition of a ray object template. As with any object that can be produced from a template, many 
exact duplicates of the object can be made from the template. This object creation from the class 
definition is termed instantiation. Thus, with the ray class, many individual rays can be created and as 
with all objects the state vector describing each object can have no influence on the state vector of any 


other object. It is seen that the formation of the beam class requires the creation of four individual ray 
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objects. In addition to the ray objects, methods must be made that manipulate the ray objects, ın 
accordance with the allowed operations on ray objects, and that produce the desired effect on the beam 
object. This object-oriented layering continues through the lobe model to the pinger and pingserver 
models. 

The only methods yet to add to all of these models are those that allow each object to sense 
interactions and perform actions based on the perceived interaction. For example, when a sound wave 
strikes a rigid surface, the wave bounces off the surface in accordance with Snell's law. Thus a method is 
added to the bottom model that allows a ray object to communicate with the bottom object. This 
communication allows the bottom object to retrieve data from the ray in regards to its position and 
direction of motion. With that retrieved information the bottom model is then able to tell if an interaction 
has occurred and can report its findings back to the ray object. The ray object is then able to alter its state 
with the reflection information provided from the bottom object. The full RRA model is then simulated 
following these rules of layering and object interaction. The full RRA code is found in Appendix A. 
Several stand-alone visualization example programs are found in Appendix B. The source code for the 
pinger and pingserver programs are found in Appendix C. Excerpts from the very useful auto-generated 


Javadoc software documentation are found in Appendix F. 


2 Visualization Implementation 


Implementation of the visualization model requires implementation of each of the visualization 
sub-models just as with the RRA implementation. Whereas implementation of each of the RRA sub- 
models is quite similar, implementation of each of the visualization sub-models is not. 

The user interface sub-model is ORE in two different ways, depending on whether or not 
the model is stand alone or interactive. The stand alone model uses the web browser user interface. Each 
desired data set to be viewed is retrieved via a web page or the web browser's file selection routines. At 
this point the interactive model requires that user input be supplied, however the standard web browser 


does not provide this kind of input. The VRML world can be programmed to perform user input, but it is 
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difficult compared to functionality already built into Java. Thus, a Java class called the battle scene server 
provides a control panel for the user to enter detailed information about the desired sonar parameters. 
This class is much the same as the classes in the RRA class package, but differs in the ability to be 
executed in a VRML web browser as a program script file. 

The data source sub-model also has its implementation determined on the basis of the form of 
user interaction. In standalone visualization, the source of data is from a file provided by the local 
operating system as requested via the web browser user interface. In interactive mode the source of data is 
read via sockets across the Internet by the battle scene server object. 

The data from the two different sources is provided to the graphics rendering model. This model 
is implemented via the VRML browser operating inside the web browser. In standalone mode 
information from the data source is passed to the VRML browser by the web browser. This data 1s parsed 
and rendered by the VRML browser and appears on the computer screen as a virtual world. In the 
interactive mode the data is parsed by the battle scene server via a method call to the VRML browser. 
This parsed data is directly inserted into the VRML browser scene graph via another method call to the 


VRML browser. The data then appears in the currently open virtual world. 


D. INTEGRATION 


1. Standalone 


Integration of the RRA and visualization model implementations into a complete simulation is 
the final task. The point of interaction and thus the point of integration in the RRA implementation is the 
pinger class. The pinger class provides for the text being written to a system file. This text is in the form 
of VRML compliant code that specifies the requested 3D visualization of the desired data, which is stored 
to a file on the computer. This file can then be viewed with the VRML browser in a semi-interactive 
manner; the user can explore the data from any perspective, but cannot interactively modify the sonar 


beam parameters. 
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2: Interactive 


In interactive mode, data transfer occurs across the Internet with the pingserver of the RRA class 
package acting as the server and the battle scene server of the visualization class package acting as the 
client. A VRML world is created, and an object in it has sonar capability and thus has a connection out of 
the world to the web browser. The battle scene server interfaces with the user, the object in the VRML 
world and the pingserver. When the user requests a sonar pulse, it tells the pingserver the specifics of the 
request and the pingserver starts calculating. When the calculations are complete the information is 
forwarded from the pingserver to the battle scene server in the form of VRML compliant text strings. The 
battle scene server then converts the text to compiled VRML code. This compiled code is sent to the 


object in the VRML world where it is automatically added to the existing scene and then rendered. 


E. CONCLUSIONS 


The choice of the programming language depended on many requirements. Java from Sun 
Microsystems met each of the requirements and was selected. It was selected for its cross-platform 
compatibility, quick execution, visualization capabilities and network capabilities. The choice of Java 
allowed the RRA model to be converted to an API for use in simulation. The highly structured object- 
oriented features of Java allowed the API to be built in a modular fashion. The ability of VRML web 
browsers to execute Java code created a tight integration between the VRML browser and the Java RRA 
API. The ability to integrate E with Java is the main reason that VRML was chosen as the rendering 
language. The integration of the API and the visualization model implementation developed into two 
programs. One that allowed interactive control of the RRA API parameters and VRML parameters and 


the other only the VRML parameters. 
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VIH. SIMULATION RESULTS 


A. INTRODUCTION 


This chapter presents the simulation results. The accuracy of the RRA API is shown to be nearly 
identical to the original RRA program as presented in [Ziomek, 1993]. In addition, the accuracy of 
energy propagation is compared to the normal mode model. This comparison proves that the RRA API is 
reasonably accurate in energy propagation. Detection of virtual world objects provided expected results 
with a two percent error rate for randomly placed mines. Visualization results were very good, but they 


barely scratch the surface of the abilities of real-time 3D sonar visualization. 


B. RRA TEMPORAL AND POSITIONAL ACCURACY 


In [Ziomek, 1993] a comparison of the output of the RRA code to that of a set of four coupled 
ordinary differential equations (ODEs) was shown. In summary it showed that the outputs down to 
several decimal places were identical for many different bottom depths and sound speed profiles. Thus it 
was shown that the RRA solution was not erroneous in its development. Also the RRA algorithm is a full 
3D propagation algorithm that allows for range-dependent environments. These features of RRA allow it 


to be a very general high resolution and high accuracy solution for a wide variety of ocean conditions. 


£. SIMULATION ACCURACY: JAVA VERSUS FORTRAN 


Ziomek’s original RRA code was written in FORTRAN while the code for this thesis is written in 
Java. Since Ziomek's RRA code was validated against an accepted set of coupled ODE's, showing that 
the implementation in this thesis compares well to his data set provides validation of the Java 
implementation. Below are several tables that show the comparison of the RRA simulation in both 


FORTRAN and Java. Because differences between the solutions are negligible, no plots of these results 
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are provided. Thus the Java implementation produced by this thesis provides valid results for Ziomek’s 


recommended test cases. 


Method 
FORTRAN 
Java 
FORTRAN 
Java 
FORTRAN 
Java 


Method 
FORTRAN 
Java 
FORTRAN 
Java 
FORTRAN 
Java 


Method 
FORTRAN 
Java 
FORTRAN 
Java 
FORTRAN 
Java 


Method 
FORTRAN 
Java 


Table VHI.1 Sound Speed Profile 1: constant speed of sound 


yo (m) 


100 
100 
100 
100 
100 
100 


c(y) 7 1500 ms 0O<y<200 m 
Bo (deg) r (km) y (m) 
45 10.0 99.95 
45 10.0 100.00 
85 10.0 174.88 
85 10.0 174.87 
135 10.0 100.05 
135 10.0 100.00 


t (s) 
9.428091 
9.428090 
6.692132 
6.692132 
9.428091 
9.428090 


Table VIIL2 Sound Speed Profile 2: Linear SSP with a positive gradient 


yo (m) 
100 
100 
100 
100 
100 
100 


c(y) = 1500 m/s + (0.017/s)y 0<y<200 m 


Bo (deg) 
45 
45 
85 
85 
135 
135 


r (km) 
10.0 
10.0 
10.0 
10.0 
10.0 
10.0 


y (m) 

99.95 

100.36 
162.57 
162.74 
100.05 
100.03 


t (s) 
9.417422 
9.417592 
6.683496 
6.683505 
9.417422 
9.417403 


Table VIII.3 Sound Speed Profile 3: Linear SSP with a negative gradient 


yo (m) 
100 
100 
100 
100 
100 
100 


yo (m) 
500 
500 


c(y) = 1500 m/s + (-0.017/s)y ک0‎ > 200 


Bo (deg) 
45 
45 
85 
85 
135 
135 


Bo (deg) 
85 


85 


r (km) 
10.0 
10.0 
10.0 
10.0 
10.0 
10.0 


c(y) =1490 m/s + (4x10?/ms)(v-500m) 


r (km) 
10.0 
10.0 
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Table VIII.4 Sound Speed Profile 4: Parabolic SSP 


y (m) t (s) 
99.95 9.438795 
100.40 9.438980 
168.96 6.699409 
168.86 6.699402 
100.05 9.438794 
100.44 9.438583 
0 <y < 1000 m 

y (m) t (s) 
777.48 6.705893 
777.45 6.705896 


s (km) 
14.1421 
1.41421 
10.0382 
10.0382 
14.1421 
14.1421 


s (km) 
14.1421 
14.1424 
10.0374 
10.0374 
14.1421 
1.41421 


s (km) 
14.1421 
14.1424 
10.0380 
10.0380 
14.1421 
14.1418 


s (km) 
10.0150 
10.0150 


FORTRAN 500 95 10.0 222.52 6.705893 
Java 500 95 10.0 222 35 6.705896 
FORTRAN 500 135 10.0 545.00 9.448829 
Java 500 135 10.0 543.48 9.449543 
Table VIII.5 Sound Speed Profile 5: Classic SSP 
c(y) = 1500 m/s + (0.016/s)y O < y <100 m 
C(y) = 1501.6m/s + (-0.02956/sX(y-100m) 100 > ۲ > 1000 m 
c(y) = 1475m/s + (0.017/s)(y-1000m) y>1000 m 
Method Yo (m) Bo (deg) r (km) y (m) t (s) 
FORTRAN 50 45 10.0 1760.21 9.604047 
Java 50 45 10.0 1758.87 9.604941 
FORTRAN 50 85 10.0 1595.33 6.806873 
Java 50 85 10.0 1596.28 6.807168 
FORTRAN 50 125 10.0 840.08 8.295966 
Java 50 125 10.0 839.14 8.296582 


10.0150 
10.0150 
14.1104 
14.1115 


s (km) 

14.2771 
14.2781 
10.1264 
10.1266 
12.3284 
12.3290 


where y in the Ziomek coordinate system is identical to depth of water. Figure 8.1 shows a representative 


deep-water SSP. Each profile also reflects the effects of a pressure release surface and rigid bottom. 
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Figure 8.1 Representative deep-water sound speed profile 
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Minor differences between the solution implemented in Java and FORTRAN are noted in the 
above tables. These differences were made smaller by choosing a smaller step size of time in the 
simulation. The reason these results were presented and not the more accurate values was execution time 
related. The smaller step-size runs had run times on the order of several minutes while the larger step- 
size had run times on the order of tens of seconds. Since the main use of the RRA API is intended to be a 
real-time simulation, execution time does become a consideration as well as accuracy. The compromise 
was to choose the quickest time that had no more than 196 error, compared to the FORTRAN 
implementation of RRA. This resulted in a time step size of 6ms. 

These simulation results show that the RRA algorithm and the specific simulation of it presented 
in this paper are accurate from a time and positional sense. The next step is to consider the accuracy of 


energy transport in the RRA API. 


D. RRA ENERGY TRANSPORT ACCURACY 


In a complex environment, especially with ocean bottoms that are not smooth when compared to 
the wavelength of the sonar pulse being used, it is a fair assumption that the interference patterns that one 
expects will be effected by strong interaction with the rough ocean surface and bottom. This interaction 
with somewhat random surfaces will in essence result in the sound waves becoming incoherent. This 
break down in coherency allows the pressures of individual waves to be superimposed on one another 
without regard to phase. This is not necessarily true if the wavelength of the sonar is comparable to the 
acoustic length of environmental objects. Since the design of this thesis was based on finding mines as 
the smallest objects, the wavelength of the sonar must be less than the size of a mine. Assuming that the 
size of a mine is one meter in diameter and the speed of sound in sea water is 1500 m/s, the minimum 
frequency to be able to resolve the mine is 1500 Hz. This and higher frequencies are typically used for 
this type of identification work. Most ocean structures tend to have acoustic lengths that are much greater 
than about one meter, but these larger structures tend to have smaller structures on them that can be near 


the size of the acoustic wavelength of the sonar. Therefore, with extensive acoustic interaction with the 


74 


ocean bottom and its structure, it is a safe assumption to use incoherent summation of pressure waves. 
However, for the sake of comparison purposes, it is necessary to compare the energy transport of the RRA 
algorithm to that of another accepted accurate model. 

To validate the RRA energy transport results, a highly simplified environment is used. This 
environment has a pressure release surface and a rigid bottom at 100 meters of water depth. Both surface 
and bottom are infinite in extent and sound speed in the water is 1500 m/s. Since the sound speed is 
constant, frequencies much lower than the typical lower limit (500 Hz, [Etter 1996]) produce accurate 
output in the RRA algorithm. Thus a frequency of 100 Hz was used as the frequency of operation for 
energy transport calculations. This low frequency was also used in order to speed up processing in the 
normal mode model, since it performs most efficiently at low frequencies. From the set up of this 
environment, it is seen that coherency must be adhered to in order to obtain accurate results. Thus for the 
time being, the real world parameters that result in incoherency are neglected and the ideal RRA model 
energy transport is compared to the normal mode model energy transport. The normal mode model was 
chosen since it is viewed as being very accurate in predicting energy transport for water columns of very 
simple geometry. 

Figure 8.1 compares the energy transport of the RRA model with the energy transport of the 
normal modes model. Notice how closely the two data sets correlate to each other. Thus coherent energy 
transport of the RRA model is similar to energy transport in the normal modes model. Now if we assume 
incoherency and can agree that this assumption has been used for years with close correlation between 
simulation and reality, then the RRA output is accurate. That is to say as accurate as assuming 
incoherency can be. Figure 8.2 shows a logarithmic trendline of the RRA model using incoherent 
summation of the pressure waves. This trendline shows the expected behavior. It is now seen that not 
only is the RRA algorithm accurate in a temporal and spatial sense, it is also reasonably accurate in 
calculating the energy transport in a simple geometry and when incoherency can be assumed in a more 


complex geometry. 
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Figure 8.2 Transmission Loss of RRA model versus Normal Mode Model 


E. DETECTION OF VIRTUAL OBJECTS 


To complete the RRA API validation the simulation needs to be used to perform some good task 
other than showing itself to produce sound pressure fields comparable to other algorithms. It needs to be 
used in a task for which it has been designed. That task being, the searching of a complex minefield in 


order to ascertain whether or not the detection capabilities of the RRA model function as desired. 


1. Mine Detection Scenario 


In this scenario a single vehicle will traverse a shallow water area of 100 meters in depth. There 
are 50 mines placed in a 1 kilometer by 1 kilometer minefield. The mines are randomly placed in all 
three dimensions. The water column is characterized by a constant speed of sound, the surface is a 
smooth pressure release surface and the bottom is a smooth horizontal rigid surface. Additionally, any 


acoustic energy leaving the 1 kilometer square is considered lost. In this test scenario, the RRA-based 
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simulated sonar system detected 49 of the 50 mines in the area. One expects that all mines are detected 
since the conditions are ideal and the only energy losses are when the sound leaves the 1 kilometer square 
area. Interestingly the reason that one of the mines was missed was due to the handling of reflections. 
Since a beam is defined as the area enclosed by four rays, when some of the rays encounter the surface 
sooner than others, there is a period of time that the beam is not traveling as it does in the real world; the 
wave front is not perpendicular to the directions of propagation of the four rays. The magnitude of this 
effect can be minimized by keeping the beam width and ping interval small which keeps the area of the 
advancing beam front small. Ideally this effect can be removed as the beam widths become arbitrarily 
small. This however introduces a massive time overhead and thus another compromise is made, this time 
between beam width and computational time. Improved algorithms and faster computers may someday 
make this small error disappear. A promising topic for future work is adaptively combining higher 
resolution, more sophisticated collision detection and smaller step size when in the vicinity of reflections. 


The code for this scenario (Pinger.java) is found in Appendix B. 


F. VISUALIZATION RESULTS 


The last two pages of this chapter contain four example sonar renderings. These renderings 
begin to show the potential of high-dimensional plots. All of the figures show the area swept out by a 
sonar beam with four degrees beam width in both the horizontal and vertical directions. The beam is 
directed at twenty degrees west of north with an eighty degree elevation from straight down. It is 
propagated for eight seconds into the sloped ocean terrain. Figure 8.3 shows the beam with transmission 
loss mapped to the color field of the VRML beam object. Figure 8.4 shows the same sonar beam with 
time mapped to color and transmission loss mapped to intensity. This mapping has the effect of drawing 
attention to only the low transmission loss portions of the beam, while giving a good feel for time and 
hence distance of the ping as it traverses the environment. Figure 8.5 shows presumed detectability and 
counter-detectability limits mapped to color and time mapped to intensity. This mapping has the effect of 


showing that the reflections coming back from the slope have little chance of counter detection by an 
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platform. Figure 8.6, which shows the same detectability limits, has its transmission loss mapped‏ ا 
to intensity. This figure more clearly shows where the area of counter detection is by allowing only the‏ 
most intense parts of the beam be seen in the later stages of propagation. Taken together, these example‏ 
figures serve to show that with some ingenuity and several mapping parameters, much more information‏ 
can be displayed in the same volume.‏ 

In conclusion, these simulation results were exactly those desired when the problem started out. 
One unexpected result was the large number of rays that need to be calculated to achieve accurate phase- 
coherent results. If incoherent summation can be assumed then the number of rays necessary drops 
considerably. However, with the continued increase in speed of general purpose computers and the 


advances in massively parallel networks, even coherent summation will be rapid. 


G. CONCLUSIONS 


The accuracy of the RRA API is shown to be nearly identical to the original RRA program as 
presented in Ziomek [1993]. In addition, the accuracy of energy propagation is compared to the normal 
mode model. This comparison proves that the RRA API is reasonably accurate in energy propagatıon. 
Detection of virtual world objects provided expected results with a two percent error rate for randomly 
placed mines. Visualization results were very good, but they barely scratch the surface of the abilities of 
real-time 3D sonar visualization. For information on obtaining a video of important simulation results 
from this thesis see Appendix G and for information on obtaining a CD ROM of the thesis contents and 


accompanying code see Appendix H. 
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Figure 8. 3 Sonar beam showing transmission loss mapped to color 
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Figure 8. 4 Sonar beam showing time mapped to color and transmission loss to intensity 
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Figure 8. 6 Sonar beam showing detection thresholds mapped to color and transmission loss to intensity 
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IX. CONCLUSIONS AND RECOMMENDATIONS 


A. PRINCIPAL THESIS CONCLUSIONS 


This chapter gathers together the many conclusions that have been reached throughout the thesıs. 
The principal conclusion of the entire thesis is that realistic real-time 3D sonar simulations and 
visualizations can be constructed with current computer hardware. Through the combination of strong 
network programming languages, Java and VRML, and the Internet, almost unlimited processing power 
can be supplied to the RRA sonar server. The sonar server coupled with a proxy server is highly 
scaleable. Adding more processors to a LAN can provide more sonar pulses per second or a higher 
degree of accuracy as needed. Widely available and inexpensive 3D graphics rendering software allow a 
deeper understanding of sonar data, as well as the ability to create immersive virtual world battle space. 
The advent of real-time 3D sonar visualization is a major breakthrough in sonar system design and 


employment. 


B. SPECIFIC CONCLUSIONS 


1. Tactical Simulation 


Tactical simulations improve significantly when using more accurate sonar information. The 
goal of tactical simulations is that of accurate prediction of real-world tactics performance, for use in 
preliminary tactics evaluation or operator training. These performance evaluations are only as good as the 
simulation that they are running on. Likewise, the simulation is only as good as the components from 
which it is composed. Evaluating tactics with only geographic constraints in mind is limited in value 
when compared to a simulation that takes in to account geographic constraints and water column 
properties that affect sonar, the primary underwater sensor. Tactical simulation results are further 


enhanced by advanced information presentation paradigms. Advanced presentations allow for more 
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intuitive and interactive, high-density information representations that lends more insight into the tactical 


performance. 


2 Sonar Simulation 


Ray theory solutions such as RRA can produce fast real-time sonar data information. These 
solutions can be characterized as three-dimensionally accurate in position, time and energy transport. The 
energy transport accuracy is verified for short-range shallow-water conditions to be nearly the same as is 
produced by normal mode theory. The same characteristics that make it ideal for accurate sonar 
information also create the right conditions for accurate detection information on targets in the virtual 
world. An especially important conclusion drawn about the RRA model is that it is highly scaleable. The 
number of processors on a given LAN can be increased without bound (as long as the LAN capacity is 
large enough to handle inter-process communications) and the sonar lobe can be divided into smaller and 
smaller pieces (limited only by the accuracy of double-precision floating-point numbers) providing for any 
desired accuracy of sonar information. Astonishingly, when proper care is taken, this can be done in real 
time (or faster than real time). Currently such capability is generally considered unfeasible on any 


computer, so this is an important result. 


3. Visualization 


In addition to tactical advantages 3D visualization of sonar information gives enhanced feel for 
space, time and intensity interactions, especially in relation to geographic features and variable sound 
speed profile (SSP). Given the high-dimension data space associated with acoustic signal processing, new 
and different forms of sonar visualization must be developed. This thesis lays the foundation for what 
advanced sonar visualization techniques might provide to tactician and researcher alike. With sonar users 
encountering larger and larger data sets, meaningful and highly condensed information need to be 
presented in a highly condensed and intuitive way, in order for the human brain to be able to process all of 


the data. These information sets are even more important to the tactician than the researcher in that 
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properly designed information rendering can allow for quicker identification and correlation, leading to 


enhanced response time and thereby improving safety of ship. 


€: RECOMMENDATIONS FOR FURTHER WORK 


1. Tactical Simulation 


With the enhanced simulation capabilities of the RRA sonar server API, the Manta UUV and 
Phoenix AUV tactical minefield searches ought to be upgraded to get better, more physically based 
information on tactics evaluations for autonomous minefield searches. Modification of the search area to 
eo porte realistic boundary conditions, geographic constraints and environmental sound speed profiles 
will lead to even more realistic tactical evaluations. Such development can lead to an optimal search 


pattern development in situ. 


2. Sonar Simulation 


The RRA sonar simulation API might benefit from a number of additions that can enhance 
realism of the model compared to real world processes. The following short paragraphs summarize the 
needed modifications, some explicitly stated in the body of the thesis and some only implied. 

A more complex bottom needs to be modeled that allows for arbitrary complexity in the physical 
constraints imposed by the bottom. Along with this increased complexity several choices for bottom 
composition need to be included that model more complex structures such as shale, clay, sand, and gravel. 
Ideally these structures can be built into a VRML 97 ElevationGrid node with additional fields added to 
indicate bottom composition at the various grid points. The effect of surface reverberation needs to be 
added to induce the all-too-present clutter in the sonar information that is generated as the acoustic energy 


interacts with the bottom. 
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A more complex surface model needs to be created that includes the disruption in the surface 
induced by wave action and micro-bubble formation. These effects have great influence on reverberation 
from the surface of the ocean. 

The current SSP model incorporates a spatially varying model that must be specified at 
compilation time. This is a severe limitation. Ideally, a model that can access a historical environment 
database of values from real-world sampling and merge it in a meaningful way with a stream of recent 
sampling updates will provide a more realistic world, one that reflects the most up to date SSP 
information. This in and of itself forms the basis of an entire thesis. To this enhanced model other 
effects need to be added, volume reverberation and noise, which further the realism of the API. The 
extreme dependence of sonar on SSP in every respect makes such work a high priority. 

Improvement to the target model offers enough of a challenge so as to warrant additional thesis 
work. The current model assumes the targets are spherical and that more complex targets (e.g. 
submarines) can be constructed of multiple adjacent spheres. This model has severe limitations when 
objects much larger or much smaller than the wavelength of the acoustic signal are encountered. 
Advancements in the processing speed of collision detection of the ping-to-target interaction needs to be 
made before even larger numbers of objects might co-exist in the virtual world. 

A final improvement might be to offer different propagation models in the RRA sonar API. This 
can make it a more general sonar API for researchers and extend the range to which the API target 
information is valid. Probable candidates include the finite element and split-step Fourier parabolic 
equation models. Both models require stepping out in range which can be converted to the time-space 
domain and can fit into the basic structure of the current RRA sonar API. Nevertheless, it remains 
important to note that RRA is very general and accurate, so integration of additional sonar models is not a 
prerequisite for widespread use. The current RRA API appears ready to serve as a computational engine 


to model most (and perhaps all) sonars currently employed by naval ships. 
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3. Visualization 


The biggest, most important and most exciting task recommended for continued work ıs the 
exploration of new and more intuitive 3D sonar data representations. A detailed study of current sonar 
data representations needs to be joined together with high-dimensional space plotting concepts that utilize 
all the capabilities of current 3D rendering hardware and software. A study of this type may prove to be of 
vital interest to the Navy and military in general. Another possibility 1s to develop a 3D target motion 
analysis system. This system takes the processed sonar data from a sonar system and through combined 
use of the RRA sonar API and knowledge of the sonar environment, statistically calculates the position of 
targets of interest. A system of this sort will have tremendous operational value. 

The final improvement to the RRA sonar API will be to separate the visualization routines from 
the RRA sonar API and to move them to a Visualization API. Many techniques which might be applied 
to a Visualization API are found in a VRML 97 generation program, [Roehl et al., 1997] Chapter 21. It is 
expected that sonar operators and researchers alike will take great advantage of the power, extensibility 


and ubiquity of VRML-based 3D graphics to better visualize and understand 3D sonar environments. 
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APPENDIX A. RRA CODE 


A. BEAM.JAVA 


/* 
Eile: Beam.java 
Compiler: Jdk1.1.6 
E 


package mil.navy.nps.rra; 


import mil.navy.nps.rra.*; 
import java.lang. *; 


per 

*(g version 1.0 

*(Qauthor LT Timothy M. Holliday (<A HREF="http://www.stl.nps.navy.mil/—auv/holliday"> 
http:/Awww.stl.nps.navy.mil/~auv/holliday</A>) | 

” 


*<dt><b>Location:</b> 

*<dd><a href= "http://www. stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/ra/Beam.java"> 
* http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/npshra/Beam.java</a> 

* 


*<dt><b>Hierarchy Diagram:</b> 

*«dd»«a href-"images/RRAClassHierarchy.gif" IMG SRC-"images/RRAClassHierarchyButton.]pg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/C oordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Summary:</b> 

*<dd>Forms a beam from four rays, propagates the beam and can retum 

* a VRML representation of the beam. 

* 

*<dt><b>Explanation:</b> 

*<dd>Beams of energy can be considered to consist of bundles of rays. The 


* energy in a bundle can be shown not to diverge from the bundle. Thus the energy 

* in a bundle is constant and thus the product of the intensity and the 

* area of the bundle perpendicular to direction of propagation is a 

* constant as well. Therefore a beam is a fundemental building block for 

* a lobe of a sonar pattern.<P> . 

* 

*<dt><b>History:</b> 

*<dd> 15Nov97 /Timothy M. Holliday [New 

*<dd> 17Mar98 /Timothy M. Holliday /Added HTML comment convention 
*<dd> 12Apr98 /Timothy M. Holliday /Parameterless Constructors 

*<dd> 14Apr98 /Timothy M. Holliday /Sumplified VRML Routines 

*<dd> 21Apr98 /Timothy M. Holliday. /Fixed Problem with calculateSoundPressureLevel 
* 

*@see Ray 

*@see BeamExampleStatic 

*@see BeamExampleDynamic 

*@see Bottom 

*@see Surface 

*@see Vec3d 

si 


public class Beam ( 
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public static final int T_L = 0; 
public static final int TIME = 1; 
public static final int NONE = 2; 


per 


* Constructor for the Beam class. 
* 
* A beam is defined as the volume swept out by 
* four rays as they traverse the ocean environment 
i 
public Beam() { 
/* 
* Initialize member variables that require it. 
n 
position = new Vec3d(); 


rayl = new Ray), 
ray2 = new Ray(); 
ray3 = new Ray(); 
ray4 = new Ray(); 


segment! = new Vec3d(); 
segment2 = new Vec3d(),; 
segment3 = new Vec3d(); 
segment4 = new Vec3d(); 


TL= new double[MAX POINTS]; 
detectTime = new double[MAX_ POINTS]; 
detectEL = new double[MAX_POINTS], 


j 


jr 

* This method resets all of the beam parameters after 
* instanciation has occurred since reuse is more time 
* efficient than garbage collection and reallocation. 
E 

public void reset() ( 


tay 1.setPosition(position.get(0),position. get(1),position.get(2)); 
ray 1.setElevation(elevationthalfBeamW idth Y); 
ray1.setAzimuth(azimuth+halfBeamW idthxX); 

tay] .setDeltaTime(deltaTime), 
rayl.setDuration(pulseDuration); 

rayl.setBottom(bottom ); 

ray1.setSurface(surface), 

rayl.setSsp(ssp); 

ray] .reset(); 
ray2.setPosition(position.get(0),position.get(1),position.get(2 )); 
ray2.setElevation(elevation-halfBeam WidthY ), 
ray2.setAzimuth(azimuththalfBeamWidthx); 
ray2.setDeltaTime(deltaTime), 
ray2.setDuration(pulseDuration); 

ray2.setBottom(bottom); 

ray2.setSurface( surface); 

ray2.setSsp(ssp); 

ray2.reset(); 

ray3.setPosition( position.get(O),position.get( 1 ),position.get(2)); 
ray3.setElevation(elevation-halfBeamWidthY ), 
ray3.setAzimuth(azimuth-halfBeamWidthX); 
ray3.setDeltaTime(deltaTime), 
ray3.setDuration(pulseDuration); 

ray3.setBottom(bottom); 

ray3.setSurface(surface); 
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ray3.setSsp(ssp); 

ray3.reset(); 

ray4.setPosition(position. get(O),position.get(1 ),position.get(2)); 
ray4.setElevation(elevation+halfBeamWidthY ); 
ray4.setAzimuth(azimuth-halfBeamWidthX ); 
ray4.setDeltaTime(deltaTime), 
ray4.setDuration(pulseDuration), 
ray4.setBottom(bottom); 
ray4.setSurface(surface); 

ray4.setSsp(ssp); 

ray4.reset(); 


if initialArea = area 1 meter from source 

/f = ] meter * solid angle (in steradians) 

if = Imeter*(2*halfBeamWidthxX * pi/180)*(2*halfBeam Width yY * pi/180) 
if = 1.21846972E-003*halfBeamWidthX*halfBeamWidthY 

ınitialArea = 1.218469 72e-3*halfBeamWidthX *halfBeamWidthY,; 

count = -1; 


detectCount = 0; 
j 
| per 
* This method sets the azımuthal angle, which is the angle from 
* the x-axis to the z-axis rotating about the y-axis. 
7 
public void setAzimuth(double phi) { 
azimuth = phi; 


y 


per 
* This method returns the azimuthal angle. 
EN 
public double getAzimuth() { 
retum azimuth; 
j 


[** 

* This method sets the half beam width of the beam in the 

* azimuthal direction. 

E) 

public void setHalfBeamWidthX(double pHalfBeamWidthX) { 
halfBeamWidthX = pHalfBeamWidthX; 

j 


per 
* This method returns the half beam width of the beam in the 
* azimuthal direction. 
*/ 
public double getHalfBeamWidthX() { 
return halfBeamWidthX; 
j 


Ge 
* This method sets the elevation angle, which is the angle from the 
* y-axis to the x-axis rotating about the z-axis . 
ef 
public void setElevation(double beta) { 

elevation = beta; 
j 
[** 
* This method returns the elevation angle. 
*/ 
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public double getElevation() { 
return elevation; 


j 


I xX xk 

* This method sets the half beam wıdth of the beam ın the 

* elevatıon direction. 

2 

public void setHalfBeamWidthY(double pHalfBeamWidthY) { 
halfBeamWidthY = pHalfBeamWidthY; 

j 


/ ۴ 
* This method returns the half beam width of the beam in the 
* azımuthal direction. 
x/ 
public double getHalfBeamWidthY() { 
return halfBeamWidthY; 
} 


/ ۷ 

* This method sets the position of the beam. 

el 

public void setPosition(double x, double y, double z) { 
position.set(x,y,Z); 

j 


JE 
* This method returns the position of the beam. 
= 
public Vec3d getPosition() { 
return position; 
j 


per 

* This method sets the time step in the simulation. 

2 

public void setDeltaTime(double pDeltaTime) { 
deltaTime = pDeltaTime; 

j 


[** 
* This method returns the simulation time step. 
*/ 
public double getDeltaTime() { 
return deltaTime; 
j 


JW 

* This method sets the simulation end time. 

x 

* This value is reletive to the start time which 1s 0.0. 

*/ 

public void setEndTime(double pEndTime) ( 
endTime = pEndTime; 

j 


/ Fx 
* This method returns the simulation end time. 
*/ 
public double getEndTime() { 
return endTime; 
j 
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[= 
* This method sets the duration of the sonar pulse. 
x 
* The integer refers to the number of deltaTime increments. 
* currently only 1 or 2 are allowed. 
El 
public void setDuration(int duration) ( 
pulseDuration = duration; 


j 


per 
* This method sets the duration of the sonar pulse. 
* 
* The integer refers to the number of deltaTime increments. 
21 
public int getDuration() { 
return pulseDuration; 
j 


پا 
This method sets the handle to the bottom object.‏ * 
4 
public void setBottom(Bottom pBottom) (‏ 
bottom = pBottom;‏ 
j‏ 
os‏ 
This method returns the handle to the bottom object.‏ * 
7 
publıc Bottom getBottom() {‏ 


return bottom; 
j 


[** 
* This method sets the handle to the surface object. 
*/ 
public void setSurface(Surface pSurface) { 

surface = pSurface; 
j 


per 
* This method returns the handle to the surface object. 
E 
public Surface getSurface() { 

retum surface; 
j 


/۴۳ 
* This method sets the handle to the sound speed profile object. 
t 
public void setSsp(SSP pSsp) { 
ssp = pSsp; 
$ 
per 
* This method returns the handle to the sound speed profile object. 
e 
public SSP getSsp() ( 
return ssp; 


j 
[** 
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* This method sets BeamNumber. It 1s used to uniquely 
* identify each beam for ROUTEing in VRML. 
e 
public void setBeamNumber(int number) ( 
beamNum = number, 
} 
[** 
* This returns BeamNumber. Which 1s used to uniquely 
* identify each beam for ROUTEing in VRML. 
u 
public int getBeamNumber() { 
return beamNum; 
j 


es 

* This method writes to the console a VRML shape that 1s 

* the three dimensional representation of the beam that ıs 

* propagated. 

er 

public String static VRML(int colorChoice, int intensityChoice) { 
String temp = ""; 


// Prepend the header information 

temp += "Transform {" + appendage ; 

temp +=" rotation 1 0 0 3.14" + appendage ; 
temp +=" children” + appendage ; 


temp += "Shape {" + appendage ; 

temp += " geometry IndexedFaceSet {" + appendage; 
temp+=" coord Coordinate {" + appendage; 
temp+=" point [ " + appendage; 

3-20, 


// Print the positions stored in each ray 
for (j=0; j<rayl.getCount(); j += 1) { 
count += 4; 
temp += rayl .position(3) + appendage; 
temp += ray2.position(j) + appendage; 
temp += ray3.position(j) + appendage; 
temp += ray4.position(3) + appendage; 
j 


temp *-" |" + appendage; 
temp += " }" + appendage; 
temp +="  coordindex (" + appendage; 


// Define the faces of the beam 

for(j=3; j<count; } += 4) { 
temp += (j-3) + m sa + " "+ (+4) + "o TED " = + appendage ; 
temp += (1-3) + ٩۱ + (+1)+ "" tre) ٩۱ + (j-2) + H ei; + appendage ; 
temp += (j-2) + ٩ ۲ + 042) + "o. 4 (j+3) + ۱۱ + (-1) * " 25 + appendage ; 
temp += (3-1) +" POTS) TOE "4+3+" -1" + appendage ; 

} 

temp += " ]" + appendage; 


temp +=" color Color{" + appendage; 
temp+="  color[" * appendage; 


// Specify a color for each group of ray positions for each time step 
// that corresponds to the Transmission Loss in dB. 
for (3=0; ¡<ray1.getCount() } += 1) { 

switch (colorChoice) { 
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case T L: Print VRML setColor Value( TL [)]); break; 
case TIME: Print VRML.setColor Value(ray].getTime(])); break; 
case NONE: Print VRML.setColor Value(1.0); break; 

j 

switch (intensityChoice) { 
case T_L: PrintVRML.setIntensityValue(TL{j]), break; 
case TIME: Print VRML.setIntensity Value(ray .getTime(j)); break; 
case NONE: Print VRML.setIntensity Value( 1.0); break; 


j 

temp *- Print VRML.getColor( ; 
j 
temp +=" ]" + appendage; 


temp +=" }" + appendage; 
temp +=" colorIndex[" + appendage; 


// Specify the vertices and their colors for each face specifed above. 
for(j=1; j<rayl.getCount(); j++) { 
temp += G-1)+" "+(-1)+""+(3)+""+(5)+" -1" + appendage ; 
temp += (j-1) +" "TOT "HUEI (+ -1" + appendage ; 
temp += (j-1) +" EE ES TOF " “FFO” -1" + appendage ; 
temp += G-1)+" SE u a ón "apu -|" + appendage ; 
j 


temp +=" ]" + appendage; 

temp +=" solid FALSE" + appendage; 
tempt+="  colorPerVertex TRUE" + appendage; 
temp +=" }" + appendage; 

temp += "$" + appendage; 

temp += "}" + appendage; 

return temp; 


qe 

* This method creates a dynamic VRML string shape that is 

* the three dimensional representation of the beam pulse that ıs 
* propogated. 

er 

public String dynamic VRMI() ( 

String beamString; 


beamString = "Transform {" + appendage + 
" rotation 1 0 0 3.14" + appendage + 
" children" + appendage + 
"Shape {" + appendage + 
" appearance Appearance {" + appendage + 
"material DEF PingColor"+beamNum+" Material {" + appendage + 
" emissiveColor 1.0 0.0 0.0" + appendage + 
" A diffuseColor O O O" + appendage +. 
"  shininess 0" + appendage + 
"  ambientiIntensity 0.2" + appendage + 
" transparency .5" + appendage + 
" }" + appendage + ۱ 
" }" + appendage + 
" geometry IndexedFaceSet {" + appendage + 
"coord DEF Ping"+beamNum+" Coordinate {" + appendage + 
" point [ " + appendage; 
JU. 


// Print the first position stored in each ray and make a box corresponding 
// to the pulse duration by printing the trailing edge of the pulse 
beamString = beamString + rayl.position(j) + " "; 
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beamString = beamString + ray2.position(j) + " "; 

beamString = beamString + ray3.position(j) +" "; 

beamString = beamString + ray4.position(j) + appendage; 
beamString = beamString + rayl.trailingPosition(j) +" "; 
beamString = beamString + ray2.trailingPosition(3) +" "; 
beamString = beamString + ray3.trailingPosition(j) +" "; 
beamString = beamString + ray4.trailingPosition(j) + appendage; 


beamString = beamStnng+" |" + appendage + 


" }" + appendage + 
۲ coordIndex [" + appendage + 


// Define the faces of the pulse 

"0374-1" +appendage + 
"045 1-1" + appendage + 
"1562 -1" + appendage + 
"2 5 7 3 -1" + appendage + 
"0 1 2 3 -1" + appendage + 
"4 5 6 7 -1" + appendage + 
" J" + appendage + 


" sold FALSE" + appendage + 
" ۱۳ + appendage + 

۲۱۳ + appendage + 

"y" + appendage + 


"DEF Propagation"+beamNum+" Coordinatelnterpolator {" + appendage + 
" key[" + appendage; 


for(j=0; j<ray].getCount(); j++) { 
beamString = beamString + rayl normalizedTime(j,endTime*2) + appendage, 

j 

beamString = beamString + " 0.51 1.0 ]" + appendage + 

" keyValue[" + appendage; 

for(j=0; j<ray1.getCount(); j++) { 
beamString = beamStnng + rayl.position(j) + " "; 
beamString = beamString + ray2.position(j) +" "; 
beamString = beamString + ray3.position(}) + " "; 
beamString = beamString + ray4.position(j); 
beamString = beamString + ray1.trailingPosition(j) + " "; 
beamString = beamStning + ray2.trailingPosition(j) + " "; 
beamString = beamString + ray3.trailingPosition(j) + * ", 
beamString = beamString + ray4.trailingPosition(]); 


) 
lz 

beamString = beamString + rayl.position(j) + " "; 
beamString = beamString + rayl.position(j) * " "; 
beamString = beamString + rayl.position(j) +" ""; 
beamString = beamString + rayl.position(j) + " "; 
beamString = beamString + rayl.position(j) + ""; 
beamString = beamString + ray1.position(j) +" "; 
beamString = beamString + rayl.position(j) +" ”; 
beamString = beamString + rayl.position(j) + " "; 
beamString = beamString + ray1.position(3) +" "; 
beamString = beamString + rayl.position(j) * " "; 
beamString = beamString + rayl.position(j) + " "; 
beamString = beamString + rayl.position(j) + " "; 
beamString = beamString + rayl.position(j) + " "; 
beamString = beamString + rayl.position) + " "; 
beamString = beamString * rayl.position(j) + " "; 
beamString = beamString + rayl.posıtion(j) + " ",; 
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beamString = beamString + " |" + appendage + 
NN + appendage + 


"DEF TLColor"+beamNum+" ColorInterpolator (" + appendage + 
" key[" + appendage; 


for(j=0; j<rayl.getCount(); j++) { 
beamString = beamString + rayl.normalizedTime(j,endTime*2), 
} 
beamString = beamString + " 0.51 1.0 ]" + appendage + 
" keyValue[" + appendage; 


// Specify a color for each group of ray positions for each time step 
// that corresponds to the Transmission Loss in dB. 
for (j=0; j<rayl .getCount(); } += 1) { 
if (TL{j] <0 ) { 
beamString = beamString + "1.0 0 0" + appendage; 
j 
else { 
beamString = beamString + (1.0-TL([5]*(.8/125)) + " 0 0" + appendage; 
j 
j 


beamString = beamString+" 000000 ]" + appendage + 
۲۱۲ + appendage + 
"ROUTE Pinglnterval.fraction_changed TO Propagation”+beamNum+".set_fraction" + appendage 
+ 
"ROUTE Propagation"+beamNum+".value_changed TO Ping"+beamNum+".set_point" + 
appendage + 
"ROUTE Pinglnterval.fraction_changed TO TLColor"+beamNum+".set_fraction" + appendage + 
"ROUTE TLColor"+beamNum+".value_changed TO 
PingColor"+beamNum+".set_emissiveColor" + appendage; 
return beamString; 
j 


[=e 

* This method calculates the trajectory of the beam of energy 

* enclosed by the defining rays of the beam tube. It also determines 
* when there 1s a detection. 

*/ 

public void calculateBeam( Targets targets) { 


durationCount = 0; 
5 


// Record first point always 
ray] .recordPoint( ), 
ray2.recordPoint( ); 
ray3.recordPoint( ); 
ray4.recordPoint( ); 

TL{j} = calculateTL(); 


it 


while (ray1.getTime() « endTime) { 


// Propagate all rays through one time step 
rayl.Propagate(.006); 
ray2.Propagate(.006 ); 
ray3.Propagate(.006 ), 
ray4.Propagate(.006 y, 


// detect and record target information 
if (targets.isCollision(rayl ,ray2,ray3,ray4)) { 
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detect Time[detectCount] = rayl.getTime(), 
detectEL[detectCount] = -2 *calculateTL()+ 10*Math.log(targets. get CollisionArea()/4/3.1415); 


3 


durationCount++, 


// If the trailing edge has left the source record all points 
if (durationCount == pulseDuration) { 
rayl.recordPoint(); 
ray2.recordPoint(); 
ray3.recordPoint(); 
ray4.recordPoint(); 
TL[j] » calculateTL(); 
y; 
j 


// If any ray has reflected then record all of the 
// points in the beam 
if ( rayl.reflected() || 
ray2 reflected() || 
ray3.reflected() || 
ray4.reflected()) { 
ray] .recordPoint(); 
ray2.recordPoint(); 
ray3.recordPoint(); 
ray4.recordPoint( ); 
TL[5] = calculateTL(); 


py 
j 


// Ifthe curvature sum reaches the limit, record 

// all of the points in the beam 

if ( rayl.totalCurvature() > CURVATURE _ LIMIT || 
ray2.totalCurvature() » CURVATURE LIMIT || 
ray3.totalCurvature() » CURVATURE LIMIT || 
ray4.totalCurvature() » CURVATURE LIMIT) ( 
ray l.recordPoint(); 
ray2.recordPoint(); 
ray3.recordPoint(); 
ray4.recordPoint(); 

TL] = calculateTLQ); 


> 


کر 


// Record last point always 
rayl.recordPoint(); 
ray2.recordPoint(); 
ray3.recordPoint(); 
ray4.recordPoint( ); 

TL{jj = calculateTLQ); 


> 


targets. resetl argets(); 
j 


/ ** 
* This method calculates the trajectory of the beam of energy 
* enclosed by the defining rays of the beam tube. 


= 

public void calculateSoundPressureLevel(double field{}{][], 
double deltaRange, 
double deltaDepth, 


double frequency) { 
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double averageX, averageY, averageReflectionPhase, deltaTime, deltaPath; 
double nearestGridX,nearestGridY; 

Vec3d averageNormal = new Vec3d(); 

int count; 


if (deltaRange>deltaDepth) ) 
deltaTime = deltaDepth/2/1500; 
deltaPath = deltaDepth/2; 

j 

else 1 
deltaTime = deltaRange/2/1500; 
deltaPath = deltaRange/2; 


j 


averageNormal.set(ray1.getNormal()); 
averageNormal.add(ray2.getNormal()); 
averageNormal.add(ray3.getNormal()); 
averageNormal.add(ray4.getNormal()), 
averageNormal.scale(.25); 


if (averageNormal.get(1) == 0.0) { 
count = (intXdeltaRange/deltaPath); 


j 
else if (Math.abs(averageNormal.get(0YaverageNormal.get(1))»deltaRange/deltaDepth) { 
count - Math.abs((intY deltaRange/deltaPath/averageNormal. get(O))); 


else { 
count — Math.abs((int)(deltaDepth/deltaPath/averageNormal.get(1))); 


j 


// Propagate a constant x distance for all rays 
while (rayl.getPositionX() « endTime*1500) 1 


deltaTime = 0.002; 


// Propagate all rays through the proper number of time steps 


for (J-0«1:]*) ( 
rayl.Propagate(deltaTime); 
ray2.Propagate(deltaTime); 
ray3.Propagate(delta Time); 
ray4.Propagate(delta Time), 
j 


averageX = (rayl.getPositionX()+ray2.getPositionX()+ray3.getPositionX()+ray4. getPositionX())/4; 
averageY = (rayl.getPosition Y()+ray2.getPosition Y( + 1ay3.getPosition Y()+ray4.getPosition Y())/4 


segment] .set(averageX average Y ,0); 

nearestGridX = (double)Math.round(averageX/deltaRange)*deltaRange; 
nearestGridY - (double)Math.round(averageY/deltaDepth)*deltaDepth; 
segment?2.set(nearestGridX ‚nearestGridY ‚0); 


if (inTheBeam(nearestGridX ‚nearestGridY ,‚0)) { 
averageNormal.set(ray1.getNormal()); 
averageNormal.add(ray2.getNormal()); 
averageNormal.add(ray3.getNormal()); 
averageNormal.add(ray4.getNormal )); 
averageNormal.scale(.25); 


// calculate distance from wavefront to storage point for phase correction 


segment? .sub(segment 1); 
a = segment2 dot(averageNormal); 
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// calculate phase correction 
averageReflectionPhase = 

(rayl.getReflectionPhase()*ray2.getReflectionPhase( )*ray3.getReflectionPhase( )*ray4. getReflectionPhase( ))/4; 
averageReflectionPhase -- 2*Math.PI*frequency*(ray1.getTime( )-a/ssp.C(rayl.getPosition())); 


// calculate the cos and sin of the phase 
b = Math.cos(averageReflectionPhase); 
c = Math.sin(averageReflectionPhase),; 


// calculate the magnitude of the pressure 
d = Math. sqrt(initialArea/calculateArea()), 


// calculate the x component of the pressure field 
field[(int)Math.round(averageX/deltaRange)][(int)Math.round(averageY/deltaDepth)][0] += d*b; 
field[(int)Math.round(averageX/deltaRange)][(int)Math.round(averageY/deltaDepth)][] 1] += d*c; 


j 
j 
j 


per 
*. This method returns the total number of detects by the beam 
n 
public int getDetectCount() ( 
return detectCount; 
j 


jes 
* This method returns the detect time for a given detect in 
* the beam 
f 
public double getDetectTime (int N) { 
return detectTime[N]; 
j 


per 
* This method returns the echo level of the detected target 
ki 
public double getDetectEchoLevel (int N) { 
return detectEL[N]; 
j 


2. ۱ 

* This method returns a VRML timer string that contains 

* the appropriate timing information for the beam 

2. 

public String pingTimerVRML() { 

return "DEF PingInterval TimeSensor {" + appendage + 

" cyclelnterval "+(endTime*2) + appendage + 
" loop TRUE" + appendage + 
"}" + appendage; 

j 


x 


This is a static method used to indicate whether a line feed 1s desired 
at the end of every line. 


— 


'true' indicates a linefeed is desired and 'false' indicates that a space 
is desired" 


+ + + بے * 


u 

public static void setAppendLineFeed(boolean pAppendLineFeed) { 
appendLineFeed = pAppendLineFeed; 
Ray.setAppendLineFeed(appendLineFged ); 
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if (appendLineFeed) { 
appendage = LINE FEED; 
} 
else { 
appendage = SPACE; 
} 
} 


ee 
* This is a static method that returns the current line appendage. 
t 
public static boolean getAppendLineFeed()( 
return appendLineFeed; 
j 


جرد دع دع بج عو بد بد بد بد ید kk oec eek e e xoc e‏ ود ید بد بعد ود ہے ہے تہ ہیر 


nt Private Section 
EEE CECI ECE CER ELE دہ‎ oco eoo oot / 
private static final String LINE. FEED - "n"; 
private static final String SPACE =" "; 
private static boolean appendLineFeed = true; 
private static String appendage = LINE FEED; 


[** 
* This method determines 1f a grid point is in the beam in the 
* current step. This is used ın calculateSoundPressureLevel. 
* Ifthe given point is in the beam then true is returned. 
i 
private boolean inTheBeam(double x, double y, double z) { 
Vec3d gridPoint = new Vec3d(x,y,z); 
return inBetween(ray1.getPosition(),gridPoint,ray3.getTrailingPosition())&& 
inBetween(ray2.getPosition( ),gridPoint,ray4.getTrailingPosition()); 
j 


سی 
This method returns a boolean indicating whether or not the grid‏ * 
point in question is in the sonar pulse.‏ * 
x‏ 
private boolean inBetween( Vec3d a, Vec3d b, Vec3d c) (‏ 
if ( ((a.get(0) 7 b.get(0) && b.get(0) 7 c.get(0)) |l(c.get(0) » b.get(0) && b.get(0) > a.get(0)))& &‏ 
((a.get(1) > b.get(1) && b.get(1) > c.get(1)) ||(c.get(1) > b.get(1) && b.get(1)» a.get(1)))&&‏ 
((a.get(2) > b.get(2) && b.get(2) > c.get(2)) ||(c.get(2) > b.get(2) && b.get(2) > a.get(2)))) {‏ 
return true;‏ 
j‏ 
else {‏ 
retum false;‏ 
j‏ 
j‏ 


۰م 
This method calculates the transmission loss at the current point.‏ * 
* 
The algorithm uses the fact that‏ * 
initialArea x initialIntensity = finalArea x ea‏ — * 
js true for a beam tube.‏ * 
/* 
private double calculateTL()(‏ 
double Attenuation 7 0;‏ 
double Area;‏ 


Attenuation = (rayl.getAbsorption() + ray2.getAbsorption() + 
ray3.getAbsorption() + ray4.getAbsorption())/4; 
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// Calculate area of half of quadralateral 

segment]! .sub(ray1.getPosition(),ray2.getPosition()), 
segment2.sub(ray2.getPosition(),ray3.getPosition()), 
segment] .cross(segment2); 

Area = .5*segment] .length(); 


// Calculate area of other half of quadralateral 
segment] .sub(ray4.getPosition(),rayl .getPosition()); 
segment sub(ray3.getPosition( ),ray4.getPosition()); 
segment] .cross(segment2); 

Area += .5*segment].length(), 


if (Area < .01){ 
Area = .01; 
j 


Attenuation += 10*Math.log(Area/initialArea)/Math.log( 10), 


return Attenuation; 
j 


JE * 
* This method calculates the area of the beam front at the current time. 
* [tis used exclusively by calculateSoundPressureLevel. 
*/ 
private double calculateArea()1 
double Area; 


// Calculate area of half of quadralateral 

segment] .sub(ray 1.getPosition(),ray2.getPosition()); 
segment2.sub(ray2.getPosition(),ray3.getPosition()); 
segment] .cross(segmentz2 ); 

Area = .5*segmentl .length(); 


// Calculate area of other half of quadralateral 
segment] .sub(ray4.getPosition(),ray1.getPosition()); 
segment2 .sub(ray3.getPosition(),ray4.getPosition()), 
segment] .cross(segment2 ); 

Area += .5*segmentl.length(); 


if (Area < .000001)4 


Area = .000001; 
} 
return Area; 
} 
Je 


* Curvature of beam that is allowed to elapse before a 
* point is recorded. This is for rendering purposes. 

ii 

private double CURVATURE LIMIT -.004; 


/* 

* The position of the ray 

x 

private Vec3d position = null; 


* 
: The azimuth of the beam 
* 
nt double azımuth = 0.0; 
JE 
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* The elevatıon ofthe beam 
7 
private double elevation = 90.0; 


/* 

* The azimuthal beam half width 

*/ 

private double halfBeamWidthX = 1.0; 


/* 

* The elevation beam half width 

El 

private double halfBeamWidthY = 1.0; 
/ ترش‎ 

* Number of iterations in length that the ray 
* segment is 

di 

private int pulseDuration = 1, 

/* 

* Handle to the bottom object 

i 

private Bottom bottom = nuil; 

/* 

* Handle to the surface object 

= 


private Surface surface = null; 


/* 

* Handle to the sound speed profile object 
=] 

private SSP ‘ssp = null; 

/۴ 

*The simulation time step 

e 

private double deltaTime = 0.006; 
per 

* Ray one for the current beam 
*/ 


private Ray rayl = null; 


|+ 

* Ray two for the current beam 
*/ 

private Ray ray2 = null; 

ER 

* Ray three for the current beam 
x 

private Ray ray3 = null; 

Dx 

* Ray four for the current beam 
*/ 


private Ray ray4 = null; 


/** 
* Segment one of the wavefront formed by 
* the four rays 
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*/ 
private Vec3d segment] - null; 


per 

* Segment two of the wavefront formed by 
* the four rays 

= 

private Vec3d segment? - null; 


TEX 

* Segment three of the wavefront formed by 
* the four rays 

= 

private Vec3d segment3 = null; 


Jas 

* Segment four of the wavefront formed by 
* the four rays 

E 

private Vec3d segment4 - null; 


pe * 

* Number of the beam 
el 

private int beamNum; 


XE 

* Area of the wavefront one meter from the source 

el 

private double initialArea = 0; // area Im from array in sq. meters 


per 

* Number of points of information that can be saved along 
* the beam path 

t 

private static final int MAX POINTS - 200; 


مب 

* Array for storing the transmission loss along the beam path 
= 

private double(] TL = null; 


JR 

* Array for storing time that a sonar detect occurs 
7 

private double[] detectTime = null; 


= 

* Array for storing the echo level of a sonar detect 
M 

private double[] detectEL = null; 

[** 

* Number of sonar detects that occur 

f 

private int detectCount = 0; 

jet 

* The pulse duration in quanta of deltaTime 
+ 

private int durationCount; 


/** 
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* End time of the simulation 
*/ 
private double endTime = 1.0; 


Ea 

* Local counters 

E 

private int count = -1; 
private int j = 0; 


per 


* Local doubles 
E 
private double a,b,c,d; 


B. BOTTOM.JAVA 


= 

File: Bottom.java 
Compiler: jdk1.1.6 

Ji 
package mil.navy.nps.rra; 
import mil.navy.nps.rra.*; 
/ Kk 

*@version 1.0 


*@author LT Timothy M. Holliday (<A HREF="http://www.stl.nps.navy.mil/-auv/holliday"> 
http://www. stl. nps.navy.mil/—auv/holliday</A>) 
* 


*<dt><b>Location:</b> 

*<dd><a href="http://www.sti.nps.navy.mildis-java-vrml/mil/navy/nps/rra/Bottom.java"> 
* http. //www.stl.nps.navy.miV/dis-java-vrml/mil/navy/nps/rra/Bottom.java</a> 

* 


*<dt><b>Hierarchy Diagram:</b> 

*<dd<a href="images/RRAClassHierarchy.gif"><IMG SRC="images/RRAClassHierarchyButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/CoordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<db<b>Summary:</b> 

*<dd>This class acoustically simulates characteristics of 

* the ocean bottom. Currently only a rigid bottom is modeled. 

1 

*<dt><b>Explanation:</b> 

*<dd>Since a ray 1s the normal vector to a plane. wave, rays interact with 

* surfaces the same way that plane waves do. Rays obey Snell's law and 

* the equations of reflection and transmission of wave energy. This model 
* uses vector algebra to determine when the end of a ray has passed through 
* the bottom and to reflect any rays that have penetrated.<P> 


* 


*<dt><b>History:</b> 

*<dd> 8Nov97 /Timothy M. Holliday /New 

*<dd> 18Mar98 /Timothy M. Holliday /Added HTML comment convention 
* 

*@see Surface 
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*(Qsee Vec3d 
27 
' public class Bottom ( 


[<= 
* Constructor for the bottom class. The first argument is the bottom type, 
* which has two choices, 'noslope' and 'slope The first choice is a flat 
* horizontal bottom and the second is a flat horizontal bottom with a steeply 
* nsing shelf and a gentle slope to shore. The second argument applies 
* only to the no slope case and 1s the depth of the bottom. 
e 
public Bottom(String bottomType, double depth) { 
if (bottomType.equalsIgnoreCase("noslope")) ( 
environment = 1; 
bottomDepth = depth; 


j 
else if (bottomType.equalsIgnoreCase("slope")) { 
environment = 2; 
shallowSlope = (10 - shelfDepth)/(5000 - shelfLine); 
shelfSlope = (shelfDepth - bottomDepth)/(shelfLine - bottomLine); 
planeVector = new Vec3d(), 
temporary Vector = new Vec3d(); 
reflectedRay = new Vec3d(), 
zone lNormal = new Vec3d(-1,1/shallowSlope,0); 
zone2Normal = new Vec3d(-1,1/shelfSlope,0), 
zone3Normal = new Vec3d(0, -1, 0); 
zone 1 Normal.normalize(); 
zone2Normal.normalize( ), 
j 
else { 
System.out.printin("User defined not implemented yet, defaulting to constant"); 
environment = 1; 
j 
j 


7 kx 

* This procedure checks to see 1f the ray path has crossed the bottom 
* during the current time step and returns true if it did. 

1 

public boolean intersected(Vec3d Pos) { 


switch (environment) { 
case |: 
if (Pos.get(1) > bottomDepth) { 
intersected = true; 
j 
else ( 
intersected = false; 
j 
break; 
case 2: 
if (Pos.get(0) > shelfLine) { 
if (Pos.get(1) > shelfDepth + (Pos.get(0)- shelfLine)*shallowSlope) { 
intersected = true; 


zone = 1; 
else { 

intersected = false; 
j 


j 

else if (Pos.get(0) > bottomLine) { 
if (Pos.get(1) > bottomDepth + (Pos.get(0)- bottomLine)*shelfSlope) { 
intersected = true; 
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zone = 2, 
} 
else { 
intersected = false; 
j 
j 


else { 
if (Pos.get(1) > bottomDepth) { 
intersected = true; 


zone = 3; 
} 
else { 
intersected = false; 
} 
} 
break; 
E intersected; 
j 
[** 


* This method causes a Snell's Law reflection to occur. 
er 
public double reflect( Vec3d Pos, Vec3d normal) { 
switch (environment) ( 
case 1: 
double diff, 
diff = Pos.get(1) - bottomDepth; 
Pos.set(1,bottomDepth-diff); 
normal.set(1,-normal.get(1)); 
break; 
case 2: 
switch (zone) { 
case 1: 
planeVector.set(shelfLine,shelfDepth,0); 
plane Vector.sub(Pos), 
intersectedDistance = -zonelNormal.dot(plane Vector)/zonelNormal.dot(normal); 
reflectedRay.set(normal), 
temporary Vector.scale(2*zone 1 Normal.dot(normal),zone1 Normal), 
normal.sub(normal,temporary Vector); 
normal normalize( ), 
reflectedRay.sub(normal ), 
reflectedRay.scale(intersectedDistance); 
Pos.sub(reflectedRay), 
break; 
case 2: 
plane Vector.set(bottomLine,bottomDepth,0); 
plane Vector.sub(Pos); 
intersectedDistance = -zone2Normal.dot(plane Vector)/zone2Normal.dot(normal); 
reflectedRay.set(normal ); 
temporary Vector.scale(2*zone2Normal.dot(normal),zone2Normal); 
normal.sub(normal,temporary Vector ); 
normal.normalize( ); 
reflectedRay.sub(normal ); 
reflectedRay.scale(intersectedDistance ), 
Pos.sub(reflectedRay), 
break; 
case 3: 
planeVector.set(bottomLine, bottomDepth,0); 
plane Vector.sub(Pos); 
intersectedDistance = -zone3 Normal.dot(plane Vector)/zone3Normal.dot(normal), 
reflectedRay.set(normal); 
temporary V ector.scale(2 *zone3 Normal.dot(normal),zone3Normal), 
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normal.sub(normal,temporary Vector); 
normal .normalize(); 
reflectedRay.sub(normal); 
reflectedRay.scale(intersectedDistance); 
Pos.sub(reflectedRay); 
break; 
j 
break; 
j 
return 0.0; 


j 


ex 
* This method returns the bottom depth at the given 
* (x,z) coordinate. (Ziomek Coordinate System) 
*/ 
public double depth(double x, double z) { 
switch (environment) { 
case 1: 
return bottomDepth; 
case 2: 
if (x < bottomLine ) { 
return bottomDepth; 


} 
else if ((x > bottomLine) && (x < shelfLine)) { 
return bottomDepth + shelfSlope*x; 
j 
else 1 
return shelfDepth + shallowSlope*(x-shelfLine); 
j 
j 
return -1; 
j 


]** 

* This method returns the VRML string representing the bottom of the ocean. 
a 

public String VRMLBottom() ( 


int minX = -5000; 
int minZ = -5000; 
int maxX = 5000; 
int maxZ = 5000; 


String bottomStnng = "Transform {" + appendage + 
" rotation 1 0 0 3.14" + appendage + 
" children[" + appendage + 
"Shape {" + appendage + 
" appearance Appearance {" + appendage + 
" material Material {" + appendage + 
"  diffuseColor .3 .15 .15" + appendage + 
" transparency 0" + appendage + 
" ۱۳ + appendage + 
" }" + appendage + 
" geometry IndexedFaceSet {" + appendage + 
" solid FALSE" + appendage + 
" coord Coordinate {" + appendage + 


" point [" + appendage; 


switch (environment) { 
case 1: 
bottomString = bottomString + " "+ maxX +" "+ bottomDepth +" "+ maxZ +"," + appendage + 
j "+ maxX +" "+ bottomDepth +" "+ minZ +"," + appendage + 
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"+ minX ۷ nds bottomDepth =." "+ minz. +" + appendage + 
"+ minX +" "+ bottomDepth +" "+ maxZ +"" + appendage + 


" . ]" + appendage + 
" }" + appendage + 
" coordIndex [ 0,1,2,3 ]" + appendage; 

break; 

case 2: 

bottomString = bottomString + " "+ maxX +" "+ "0" +" "+ maxZ +"," + appendage + 
: "+ maxX +" "+ "0" +" "+ minZ +"," + appendage + 
: "+ shelfLine +" "+ shelfDepth +” "+ minZ +"," + appendage + 
i "+ shelfLine +" "+ shelfDepth +" "+ maxZ +"," + appendage + 
5 "+ bottomLine +" "+ bottomDepth +" "+ maxZ +"," + appendage + 
i "+ bottomLine +" "+ bottomDepth +" "+ minZ +"," + appendage + 
"+ minx +" "+ bottomDepth +" "+ minZ +"," + appendage + 
i "+ minX +" "+ bottomDepth +" "+ maxZ+ "" + appendage + 
"J" + appendage + 
" ۱۳ + appendage + 
"  coordindex [ 0,1,2,3,-1,3,2,5,4,-1,4,5,6,7,-1 ]" + appendage; 

break, 


bottomString 7 bottomString +" 3" + appendage + 
"3" + appendage + 


"]' * appendage+ 
")" + appendage+ 
VRMLScales(); 
return bottomStnng; 
j 
/ xx 


* This method applies scales to the VRML bottom 
* for the x, y and z directions. 


y 
public String VRMLScales() ( 
String ts — ""; 
ts *- Print VRML.scale("Meters", 
5000,-5000, 
0, -bottomDepth, 5000, 
0,0, 1,0, 
10), 
ts += Print VRML.scale("Meters", 
bottomDepth, 0, 
-5000, 0, 5000, 
0, 0, 1, -1.57, 
2); 
ts += PnntVRML.scale("Meters", 
-5000, 5000, 
-5000, -bottomDepth, 0, 
0, 1, 0, -1.57, 
10); 
return ts; 
j 
]** 


This is a static method used to indicate whether a line feed is desired 
at the end of every line. 


‘true’ indicates a linefeed is desired and ‘false’ indicates that a space 
is desired" 


E Y A &‏ دہ 


*/ 
public static void setAppendLineFeed(boolean pAppendLineFeed) { 
appendLineFeed = pAppendLineFeed; 
if (appendLineFeed) { 
appendage = LINE FEED; 
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j 


else { 
appendage = SPACE; 
j 


[es 
* This is a static method that returns the current hne appendage. 
*/ 
public static boolean getAppendLineFeed()( 
return appendLineFeed; 
j 
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rz Private Section 
3x Voc ox ooo oec oe doe oe OE E E E E EL E N VL رع‎ 


private static final String LINE. FEED - "in"; 
private static final String SPACE - " 

private static boolean appendLineFeed = true; 
private static String appendage = LINE FEED; 


/** standard environment that was initialized */ 
private int environment = 1; 


/** zone within the environment where the ray currently is */ 
private int zone = 0; 


/** depth of the ocean floor in a simple sloped model */ 
private double bottomDepth = 2000; 


/** depth of the ocean shelf in a simple sloped model */ 
private double shelfDepth = 500; 


/** slope of the plane from the shelf to the ocean bottom */ 
private double shelfSlope = 0.0; 


/** slope of the plane from the shore to the ocean shelf */ 
private double shallowSlope = 0.0; 


/** location of the line in the x direction where the ocean shelf ends */ 
private double shelfLine = 2500.0; 


/** location of the line in the x direction where the ocean floor ends */ 
private double bottomLine = 0.0; 


/** distance from ray intersection with the bottom to the current ray position */ 
private double intersectedDistance = 0.0; 


/** vector from a point on the reflection plane to the current ray position */ 
private Vec3d planeVector = null; 


/** temporary vector used in calculating the reflection */ 
private Vec3d temporaryVector = null; 


/** the vector expressing the ray that is reflected from the bottom */ 
private Vec3d reflectedRay = null; 


/** the normal to the bottom in zone 1 in the simple sloped model*/ 
private Vec3d zonel Normal - null; 


/** the normal to the bottom in zone 2 in the simple sloped model*/ 
private Vec3d zone2Normal = null; | 
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/** the normal to the bottom ın zone 3 1n the simple sloped model*/ 
private Vec3d zone3Normal = null; 


/** temporary boolean variable indicating 1f a reflection has occured*/ 
private boolean intersected; 


E: LOBE.JAVA 


/ * 
File: Lobe.java 
Compiler: Jdk1.1.6 
= 
package mil.navy.nps.rra; 
import mil.navy.nps.rra.*; 
import java.util Date; 
er 


*(Q version 1.0 

*(Qauthor LT Timothy M. Holliday (<A HREF="http://www. stl.nps.navy.mil/—auv/holliday"> 
http://www. stl.nps.navy.mi/—auv/holliday</A>) 

* 


*<dt><b>Location:</b> 

*<dd><a href="http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/tra/Lobe.java"> 
* http //www.stl.nps.navy.miV/dis-java-vrml/miVnavy/nps/rra/Lobe.java</a> 

* 


*<dt><b>Hierarchy Diagram:</b> 

*<dd<a href="images/RRAClassHierarchy.gif"><IMG SRC="images/RRAClassHierarchyButton.jpg" 
ALIGN=ABSCENTER></a> 

x 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem gif"><IMG SRC="images/CoordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Summary:</b> 

*<dd>This class forms a lobe by defining multiple beams. 

* 

*<dt><b>Explanation:</b> 

*<dd>Lobes are emitted in a pattern from an source. Ina 
properly designed active sonar, there 1s only one main 

lobe. This class simulates that lobe. Lobes are characterized 
by a direction that they are pointed in and typically a 
horizonatal and vertical beamwidth. Within this lobe resides 
the acoustic energy used for detection of distant objects. 

This class subdivides the lobe into beams. This 1s done to 
keep the wave front of each individual bundle of rays from 
getting significantly larger than the wavelength of the acoustic 
energy. This constraint 1s so that accuracy is maintained and 
so that many special problems with large area bundle reflection 
can be ignored. Thus, the lobe class divides the lobe into a 
matrix of beams and divides the energy in the lobe among the 
beams. Each beam is then calculated individually using the beam 
class.<P> 
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*<dt><b>Future Work:</b> 
*<dd>To optimize the execution time, repeated calculation of duplicate rays 
* needs to be eliminated. 


* 
*<dt><b>History:</b> 
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*<dd> 29Nov97 /Timothy M. Holliday /New 


*«dd» 17Mar98 /Timothy M. Holliday /Added HTML comment convention 
*<dd> 12Apr98 /Timothv M. Holliday /Parameterless Constructors 
*<dd> 14Apr08 /Timothy M. Holliday /Simplified VRML Routines 
* 

*@see Beam 

*@see ExampleLobeStatic 

*@see ExampleLobeDynamic 

*@see Bottom 

*@see Surface 

*@see Vec3d 

t 


public class Lobe 1 


pus 
* Constructor for the lobe class. 
¥ 
* A lobe is defined as an mxn array of beams. Due to memory and 
* initialization constraints the current limit is 5 by 5. 
7 
public Lobe() { 
for (1=0;i<MAX _X PARTITIONS;i++) { 
for =0;)}<MAX_Y_PARTITIONS;j++) { 
mainLobe[i][j]] = new Beam(); 
j 
j 
j 


er 


* This method resets all of the lobe parameters after 

* instantiation has occurred since reuse is more time 

* efficient than garbage collection and reallocation. 

* 

a 

public void reset() { 
beamWidthX = lobeWidthX/(double)numberXPartition; 
beamWidthY = lobeWidthY/(double number YPartition; 
double halfLobeWidthX = lobeWidthX/2.0; 
double halfLobeWidthY = lobeWidthY/2.0; 


for (1=0; i<numberXPartition; 1++) f 
for (j=0; }<number YPartition; j++) { 
mainLobe[i][j].setPosition(position.get(0), position.get(1), position. get(2)); 
mainLobe[{i][}].setElevation(elevation-halfLobeWidth Y+(double)(2*1+1 )/2.0*beamWidthY), 
mainLobe[i][j].setAzimuth(azimuth-halfLobeW idthxXx+(double\2*}+1)/2.0*beamWidthx); 
mainLobe{i][}].setHalfBeam Width Y(beam Width Y/2.0); 
mainLobe[i][5].setHalfBeamWidthX(beamWidthX/2.0); 
mainLobe{i}[j].setDuration(pulseDuration), 
mainLobe[i][}].setEndTime(endTime); 
mainLobe[i][j].setBottom(bottom), 
mainLobe[1][j].setSurface(surface); 
mainLobe[i][}].setSsp(ssp); 
mainLobe[i][j].setBeamNumber(i*MAX_X PARTITIONS+)); 
mainLobe[1][j].reset(); 
j 
j 
j 


Js 


* This method sets the azimuthal angle, which is the angle from 
* the x-axis to the z-axis rotating about the y-axis. 
*/ 
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public void setAzimuth(double phi) { 
azımuth = phi; 
} 


/** 
* This method to gets the Azimuthatl Angle. 
*/ 
public double getAzimuth() ( 
return azimuth; 
j 


Er 

* This method to sets the total lobe wıdth 

* in the azimuthal direction reletive to the Ziomek 

* Coordinate System. Argument is in degrees. 

1 

public void setLobeWidthX(double pLobeWidthX) ( 
lobeWidthX = pLobeWidthX, 

j 


per 
* This is an accessor method to get the total lobe width 
* in the azimuthal direction reletive to the Ziomek Coordinate 
* System. Value is in degrees. 
dl 
public double getLobeWidthX() ( 
return lobeWidthX; 
j 


om 
* This method sets the elevation angle, which is the angle from the 
* y-axis to the x-axis rotating about the z-axis . 
1 
public void setElevation(double beta) { 
elevation = beta; 
j 


per 
* This method returns the elevation angle. 
i 
public double getElevation() { 
return elevation; 


j 


pur 

* This is an accessor method to set the total lobe width 

* in the elevation direction reletive to the Ziomek 

* Coordinate System. Argument is in degrees. 

$ 

public void setLobeWidthY(double pLobeWidthY) { 
lobeWidthY = pLobeWidthY; 

} 


per 
* This is an accessor method to get the total lobe width 
* inthe elevation direction reletive to the Ziomek Coordinate 
* System. Value is in degrees. 
E 
public double getLobeWidthY() { 
return lobeWidth Y; 
j 


/** 
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* This method sets the lobe positon. 

sf 

public void setPosition( Vec3d pPosition) ( 
position.set(pPosition ); 


/** 
* This method returns the lobe position. 
=] 
public Vec3d getPosition() { 
return position; 
j 


per 

* This method sets the time step in the simulation. 

7 

public void setDeltaTime(double pDeltaTime) { 
deltaTime = pDeltaTime; 

j 


per 
* This method returns the simulation step time. 
ay 
public double getDeltaTime() { 
return deltaTime; 
} 


تاد 3 / 

* This method sets the number of beams 1n the azimuthal 

* direction that there are in the lobe. 

*/ 

public void setNumberXPartition(int number) { 
numberXPartition = number, 

} 


/ * xk 
* This method returns the number of beams in the azimuthal 
* direction that there are in the Jobe. 
el 
public int getNumberXPartition() { 
return numberXPartition; 
j 


ر7 
This method sets the number of beams in the elevation‏ * 
direction that there are in the lobe.‏ * 
/* 
public void setNumberYPartition(int number) {‏ 
numberY Partition = number,‏ 
j‏ 


a 
* This method returns the number of beams ın the elevation 
* direction that there are in the lobe. 
*/ 
public int getNumberYPartition() { 
return number Y Partition; 
j 


Jer 


* This method sets the simulation end time. 
x 


* This value is reletive to the start time which ıs 0.0. 
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= 
public void setEndTime(double pEndTime) { 
` endTime = pEndTime; 


j 


je" 
* This method returns the simulation end time. 
*/ 
public double getEndTime() { 
retum endTime; 
j 


pes 
* This is an accessor method. Duration is currently | or 2. 
* The integer refers to the number of deltTime increments. 
"1 
public void setDuration(int duration) ( 

pulseDuration = duration; 
j 


میا 
This is an accessor method. Duration is currently 1 or 2.‏ * 
The integer refers to the number of deltTime increments.‏ * 
El‏ 
public int getDuration() {‏ 

return pulseDuration; 
j 


per 
* This method sets the handle to the bottom object. 
= 
public void setBottom(Bottom pBottom) { 
bottom = pBottom; 
} 5 


[** 
* This method returns the handle to the bottom object. 
di 
public Bottom getBottom() { 
retum bottom; 
j 


TEF 
* This method sets the handle to the surface object. 
1 
public void setSurface(Surface pSurface) { 
surface = pSurface; 
} 


jt 
* This method returns the handle to the surface object. 
x 
public Surface getSurface() ( 
return surface; 


yet 
* This method sets the handle to the sound speed profile object. 
x 
public void setSsp(SSP pSsp) { 
ssp = pSsp; 
j 
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وه 
This method returns the handle to the sound speed profile object.‏ * 
” 
public SSP getSsp() {‏ 
return ssp;‏ 
j‏ 


JE 
* This method calculates the Lobe by calling the 
* calculateBeam method for each beam in the lobe. 
el 
public void calculateLobe(Targets targets) { 


for (i=0; i<numberYPartition; 1++) { 
for (j=0; j<numberXPartition; j++) { 
mainLobe[i |[j].calculateBeam(targets); 


pee 


* This method returns VRML objects that signify detections that 
* were made in the virtual world. 

x 

* The objects are collections of red spheres indicating a detect. 
f 

public String detection VYRML() { 


double lAzimuth, lElevation, x, y, z, range, EL; 


String temporaryString = 
"Group {" + appendage+ 
" children [" + appendage; 


for (170; 1XnumberY Partition; 17) { 
for (j=0; j<numberXPartition; j++) { 
for(k=0;k<mainLobe[i ][}].getDetectCount();k++) { 
1Azimuth = azimuth + beamWidthX/2.0*(2*(double)j+1~double )numberXPartition); 
lAzimuth *= (3.1415/180), 
lElevation = elevation + beamWidth Y/2.0*(2*(double)i+1-(double)numberY Partition); 
lElevation *= (3.1415/180); 
range = mainLobe[1][}].getDetect Time(k)* 1 500; 


// Convert from Ziomek to VRML coordinate system 

x = position.get(0) + range*Math.sin(lElevation)*Math.cos(lAzımuth), 

y = -(position.get(1) + range*Math.cos(lElevation)), 

Zz = -(position.get(2) + range*Math.sin(lElevation)*Math.sin(lAzimuth)); 


EL = mainLobe[1][j].getDetectEchoLevel(ky-200; 
if(EL > 0.9) { 
EL = 0.9; 


j 

else if (EL <0) { 
EL = 0; 

j 


temporaryString += 
" Transform {" + appendaget 
" translation "+x+" "+y+" "+z+ appendage + 
" children [" + appendage+ 
" Shape {" + appendage+ 
" appearance Appearance ("  appendage* 
" materıal Material {" + appendage+ 
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" emissiveColor 1 0 0" + appendage+ 
" transparency "+EL+ appendage + 
" ۱۳ + appendage+ 
" }" + appendaget 
" geometry Sphere { radius 20 }" + appendage+ 
" }" + appendaget 
n" y + appendage+ 
"3" + appendage; 
j 
j 
j 
return temporaryString + "] } " + appendage; 
$ 


¥ 

* This method writes the dynamic VRML representation of the Lobe by calling the 
* VRMLBeam routine for each beam. 

* Also Writes out a Viewpoint node. 

=) 

public String dynamic VRML() { 


String lobeString = "Viewpoint {" + appendage + 
" position " + (position.get(0)-100) +" "+ -position.get(1) +" "+ (position. get(2)}-20}+ '"" + appendage + 
" orientation 0 -1 0 " + ((90-azimuth)*3.14/180}+"" + appendage + 
" description \"Beam\" " + appendage+ 
"}" + appendage; 


lobeString = lobeString + mainLobe[0][0].pingTimer VRML(); 


Date tumecheck = new Date(); 
for (1=0; ¡<number YPartition; 14+) ( 
for (j=0; j<numberXPartition; j++) { 
lobeString = lobeString + mainLobe[i][j].dynamicVRML(); 
j 


j 
Date timecheck2 = new Date(); 


lobeString = lobeString + "# " + timecheck + "" + appendage; 
lobeString = lobeString + "# " + timecheck2 + "" + appendage; 
return lobeString; 


$ 


[** 
* This method writes the static VRML representation of the Lobe by calling the 
* VRMLBeam routine for each beam. 
*/ 
public String staticVRML(int colorMap, int intensityMap) { 


String temp = ""; 


for (i=0; i<numberYPartition; 1++) { 
for (j=0; j<numberXPartition; j++) { 
temp += maimLobe[1]{j].static VRML( colorMap, intensityMap), 
j 
j 
return temp; 


j 


pe * 
* This 1s a static method used to indicate whether a line feed 1s desired 
at the end of every line. 


* * 


* ‘true’ indicates a linefeed is desired and 'false' indicates that a space 
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* is desired" 
74 
public static void setAppendLineFeed(boolean pAppendLineFeed) { 
appendLineFeed = pAppendLineFeed; 
Beam.setAppendLineFeed(appendLineFeed); 
if (appendLineFeed) { 
appendage = LINE FEED; 
$ 
else ( 
appendage = SPACE; 
j 
j 


per 
* This is a static method that returns the current line appendage. 
f 
public static boolean getAppendLineFeed() { 
return appendLineFeed; 
j 
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eS Private Section 
رع ہے ہے کہہے ہے ےکک‎ 
private static final String LINE_FEED = "\n", 
private static final String SPACE - " "; 
private static boolean appendLineFeed = true; 
private static String appendage = LINE FEED, 


private static final int MAX_X_PARTITIONS = 5; 
private static final int MAX_Y_PARTITIONS = 5; 


private Vec3d position = new Vec3d(); 
private double azimuth = 0.0; 
private double elevation = 90.0; 
private double lobeWidthX = 1.0; 
private double lobeWidthY = 1.0; 
private double beamWidthX = 1.0; 
private double beamWidthY = 1.0; 
private int numberXPartition = 1; 
private int number YPartition = 1; 
private int pulseDuration = 1; 
private Bottom bottom = null; 
private Surface surface = null; 
private SSP ssp = null; 

private double deltaTime = 0.006; 
private double endTime; 


private Beam[][] mainLobe = new Beam[MAX X PARTITIONS][MAX Y PARTITIONS], 
private int IK 


D. PRINTVRML.JAVA 


pe 
File: Print VRML. java 
Compiler: jdk1.1.6 
*/ 


package mil.navy.nps.rra; 
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import mil.navy.nps.rra.*; 


7" 

*@version 1.0 

*@author LT Timothy M. Holliday (<A HREF="http//www. stl. nps.navy.mil/-auv/holliday"> 
http://www. stl nps.navy.mil/~auv/holliday</A>) 

* 


*<dt><b>Location:</b> 

*«dd»a href-"http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/Print VRML java"> 

* http;//www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/Print VRML.java«/a» 

x 

*<dt><b>Hierarchy Diagram:</b> 

*<dd><a href="images/RRAClassHierarchy.gif"><IMG SRC="images/RRAClassHierarchyButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/CoordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

E 


*<dt><b>Summary:</b> 
*<dd>Contains static methods that return standard VRML strings. 
* 


*<dt><b>Explanatıon:</b> 
*<dd>Creates strings in VRML for header, directionalLight and 
* navigationinfo.<P> 


* 


*<dt><b>History:</b> 
*<dd> 15Apr98 /Timothy M. Holliday /New 
*<dd> 14May98 /Timothy M. Holliday /Added color and intensity schemes 


*<dd> 18May98 /Timothy M. Holliday /Added VRML colorBarView and textNode 
* 


EJ 
public class Print VRML ( 


public static final int RAINBOW = 1; 

public static final int RED = 2; 

public static final int GREEN = 3; 

public static final int BLUE = 4; 

public static final int TRICOLOR = 5; 

public static final int CONSTANT = 6; 

public static final int LINEAR = 7, 

public static final String HORIZONTAL ="0 0 1 0"; 
public static final String VERTICAL = "00 1 1.57"; 
public static final String RIGHT = "\"END\""; 
public static final Sting CENTER = "\"MIDDLE\""; 
public static final String LEFT = "BEGIN"; 


public static String header() { 
return "#VRML V2.0 utf8" + appendage, 
j 


public static String directionalLight() ( 
return "DirectionalLight {"+appendage+ 

" ambientIntensity 0"+appendage+ 
* color 11 1"+appendage+ 
" direction O -1 0"+appendage+ 
" intensity 1“+appendage+ 
" on TRUE"+appendage+ 
"}"+appendage; 

j 


public static String navigationInfo() ( 
return "Navigationinfo [ "+appendage+ 
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" type [\"FLY\" \"EXAMINE\"]"+appendage+ 
" speed 50"+appendage+ 
" headlight TRUE"+appendage+ 
" visibilityLimit 50000"+appendage+ 
"}"+appendage; 

j 


public static String printLegendView(String title, 

String infol, 

String info2, 

String info3, 

String info4, 

String info5, 

double x, 

double y, 

double z) { 
String temp = appendage; 
temp += "EXTERNPROTO protoLegendViewpoint[" +appendage; 
temp += "eventin  SFBool  set_bind" +appendage; 
temp += "exposedField SFFloat fieldOfView " +appendage; 
temp += "exposedField SFBool jump " +appendage; 
temp += "exposedField SFRotation orientation " +appendage; 
temp += "exposedField SFVec3f position "+appendage; 
temp += "field SFString description" *appendage; 
temp += "eventOut  SFTime bindTime" +appendage; 
temp += "eventOut SFBool ısBound" +appendage; 
temp += "field MFString Title " +appendage; 
temp += "field MFString info " +appendage; 
temp += "] \"protoLegendVıiewpoint.wri\"" +appendage; 


temp += "protoLegend Viewpoint {" +appendage; 
temp += "Title [\""+ title +"\"]" +appendage; 
temp += "position ی‎ NE^ "TzT"" +appendage; 
temp += "info [V* +1nfol+ mM ee +info2+ m پیا‎ +info3+ un +appendage, 
temp += " XO +1nfo4+ TOS Jun +info5+ my" +appendage; 
temp += "3" +appendage, 
return temp; 
j 


public static String protoColorBar(String vLabel, 
String hLabel, 
String vSubLabell, 
String vSubLabel2) { 
String temp = appendage; 
temp += "PROTO protoColorBar[" +appendage; 
temp += "] ("+appendage; 


temp += "Transform {"+appendage; 
temp += "translation 2 1 4.5"+appendage; 
temp += "children ["+appendage; 
temp t7 colorBar(vLabel,hLabel,vSubLabell ,vSubLabel2 ); 
temp += "]" +appendage; 
temp += "}" +appendage; 
temp += "}" +appendage; 
return temp; 
j 


public static String colorBar View(String title, 
double x, 
double y, 
double z, 
double xAxis, 
double yAxis, 
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double zAxis, 
double theta) { 


String temp = "Transform {" +appendage; 
temp += "translation "+x+" "+y+" "+z+"" +appendage; 
temp += "rotation "+xAxist" "+yAxist" "+zAxist" "+theta +appendage; 
temp += "children [" +appendage; 
temp += "Viewpoint {" +appendage; 
temp += "description \""+ title +"\"" +appendage; 
temp += "}" +appendage; 
temp += "protoColorBar {}"; 
temp += "]" +appendage; 
temp += "3" +appendage; 
return temp; 
j 


public static void setAppendLineFeed(boolean pAppendLineFeed) ( 
appendLineFeed = pAppendLineFeed; 
Ray.setA ppendLineFeed(appendLineFeed), 
if (appendLineFeed) { 
appendage = LINE FEED; 


else { 
appendage = SPACE; 
j 
j 


public static boolean getAppendLineFeed() { 
return appendLineFeed; 


public static void setColorScheme(int scheme) ( 
colorMaximum = 1; 
colorMinimum = 0; 
colorScheme = scheme; 
colorDelta = 1; 
colorSubDelta = colorDelta/3; 


j 


public static void setColorScheme(int scheme, double maximum, double minimum) { 
colorMaximum = maximum; 
colorMinimum = minimum; 
colorScheme = scheme; 
colorDelta = maximum - minimum; 
colorSubDelta = colorDelta/3; 
j 


public static void setColorScheme(int scheme, 
double maximum, 
double detection, 
double counterDetection, 
double minimum) { 
colorScheme = scheme; 
colorMaximum = maximum; 
colorDetection = detection; 
colorCounterDetection = counterDetection; 
colorMinimum = mınımum; 
colorDelta = maximum - minimum; 
j 
public static void setColorValue(double value) { 


if ((colorScheme — RAINBOW/)ll(colorScheme -- TRICOLOR)) í 
colorValue = (value-colorMinimum)/colorDelta; 
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if (colorValue<0.0) ( 
color Value = 0.0; 


} 
else if(colorValue»1.0) { 
colorValue = 1.0, 


j 


j 
/lelse if (colorScheme == TRICOLOR) { 
// colorValue = value; 
IR 
else { 
// _System.out.printin("#PrintVRML: setColorValue: not a valid scheme"); 
colorValue = 1; 
j 
j 


public static void setIntensityScheme(int scheme) { 
intensityScheme = scheme; 
intensityMaximum = 1; 
intensityMinimum = 0; 
intensityDelta = 1; 

j 

public static void setIntensityScheme(int scheme, double maximum, double minimum) ( 
intensityScheme = scheme; 
intensityMaximum = maximum; 
intensityMinimum - minimum; 
intensityDelta = maximum - minimum; 


j 


public static void setIntensityValue(double value) { 
if (intensityScheme == LINEAR ) { 
intensity Value = (value-intensitvMinimum)/intensityDelta; 
if (intensity Value < 0.0) { 
intensity Value = 0.0; 


j 

else if(intensity Value > 1.0) { 
intensityValue = 1.0; 

j 


j 

else if (intensityScheme == CONSTANT) { 
intensity Value = 1.0; 

} 

else { 
System.out.printin("#PrintVRML: setIntensityValue: not a valid scheme"); 
intensityValue = 1.0; 

j 

j 


public static String getColor() { 
double red = 0; 
double green = 0; 
double blue = 0; 


if (colorScheme == RAINBOW) { 
if ( color Value > .66667) { 
red = (colorValue - .66667)/.33333; 
green = 1-red; 
blue = 0.0; 


j 
else if (colorValue > .33333) { 
red = 0.0; 
green = (color Value - .33333)/.33333; 
blue = 1-green; 
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else { 


green = 0.0; 
blue = colorValue/.33333; 
red = 1-blue; 


j 


j 

else if (colorScheme == RED) { 
red = 1.0; 
green = 0.0; 
blue = 0.0; 


j 

else if (colorScheme == GREEN) { 
red = 0.0; 
green = 1.0; 
blue = 0.0; 


j 
else if (colorScheme == BLUE) ( 
red = 0.0; 
green = 0.0; 
blue = 1.0; 
2 
else if (colorScheme == TRICOLOR) { 
if (colorValue > (colorDetection-colorMinimum)/colorDelta) { 
red = 1.0; 
green = 0.0; 
blue = 0.0; 


j 
else if (color Value > (colorCounterDetectionolorMinimum)/colorDelta) { 
red = 0.0; 


green = 1.0; 
blue = 0.0; 
} 
else { 
red = 0.0; 
green = 0.0; 
blue = 1.0; 
} 
j 


red *- intensity Value; 
green *- intensity Value; 
blue *- intensity Value; 


return (StringXred+" "+green+" "+blue+appendage); 


j 
public static String colorBar(String vLabel, 
String hLabel, 
String detectionLabel, 
String counterDetectionLabel) { 
double i = 0; 
double j = 0; 
String ts= an 


switch (colorScheme) { 
case RAINBOW: 
ts += textNode( VERTICAL ,-.75, 0, 0, vLabel, .25,.25, CENTER); 
ts += textNode(HORIZONTAL..6, .4, 0, ""+colorMaxımum, .15,.15,LEFT); 
ts += textNode(HORIZONTAL..6, -.4, 0, ""+colorMinimum, .15,.15,LEFT); 
break; 
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case TRICOLOR: 

ts += textNode(HORIZONTAL -.6, .5+(colorDetection-colorMaxımum)/colorDelta, 
0, detectionLabel, .15,.15,RIGHT), 

ts += textNode(HORIZONTAL.,-.6, .5+(colorCounterDetection-colorMaximum )/colorDelta, 
0, counterDetectionLabel, .15,.15,RIGHT), 

ts += textNode(HORIZONTAL.,.6, .5, 0, ""+colorMaximum, .15,.15,LEFT),; 

ts += texiNode(HORIZONTAL ,.6, .5*(colorDetection-colorMaximum YcolorDelta, 
0, ""+colorDetection, .15,.15,LEFT),; 

ts += textNode(HORIZONTAL,.6, .5HcolorCounterDetection-colorMaximum YcolorDelta, 
0, ""+colorCounterDetection, .15,.15,LEFT), 

ts += textNode(HORIZONTAL.,.6, -.5, 0, ""+colorMinimum, .15,.15,LEFT); 

break: 


case RED: 
ts += textNode(VERTICAL.,-.75, 0, 0, vLabel, .25,.25, CENTER y; 
break; 


case BLUE: 
ts += textNode( VERTICAL.-.75, 0, 0, vLabel, .25,.25,CENTER), 
break; 


case GREEN: 
ts += textNode(VERTICAL,-.75, 0, 0, vLabel, .25,.25,CENTER); 
break; 


j 


switch (intensityScheme) { 
case LINEAR: 
ts += textNode(HORIZONTAL,O , .75, 0, hLabel, .25,.25, CENTER); 
ts += textNode( VERTICAL -.4, -.6, 0, ""+intensityMaximum, .15,.15,RIGHT), 
ts += textNode( VERTICAL ..4, -.6, 0, ""+tintensityMinimum, .15,.15,RIGHT); 
break; 


case CONSTANT: 
ts += textNode(HORIZONTAL,0 , .75, 0, hLabel, .25,.25, CENTER); 
break; 
j 


ts += "Shape {"+appendage; 
ts += "geometry IndexedFaceSet {"+appendage, 
ts += "coord Coordinate {"+appendage; 
ts += "point ["+appendage; 
for (1=0;1<20;i++) { 
for (j=0:3<20:j++) { 
ts += (-.5Hdouble)1/20) +" "; 
ts += (.5-(double)j/20) +" 0.0"+appendage; 
j 
j 
ts += "]"+appendage; 
ts += "}"+appendage; 


ts += "coordIndex ["+appendage; 
for (1=0;1<19,14>+) ( 
for (3=0;3<19:3++) ( 
ts += (imtG*20+1+" "HintxG+1)*20+)+" "HintxG+1)*20+1+1)+" "+(mtG*20+1+1)+" -1"+appendage; 
j 


j 
ts += "]"+appendage; 


ts += "color Color {"+appendage; 
ts += "color ["+appendage; 
for (1=0;1<20;14+) ( 
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]01۲ )120:1>20:1+7( ) 
setIntensity Value(intensityMaximum-(double /20*1ntensityDelta); 
setColor Value(colorMaximum-(double)j/20* colorDelta ), 
ts += getColor(); 
j 
j 
ts += "]"+appendage; 
ts += "}"+appendage; 


ts += "colorIndex ["+appendage; 
for (i=0;1<19;1++) { 
for (3=0:3<19:++) ( 
ts += (int)G*20+1)+" "+(intxG+1)*20+1)+" "+(nt(G+1)*20+1+1)+" "+(intG*20+1+1)y+" -1"+appendage, 
j 
j 
ts += "|"+appendage; 
ts += "solid FALSE"+appendage; 
ts += "colorPerVertex TRUE"+appendage; 


ts += "}"+appendage; 
ts += "}"+appendage; 
return ts; 


j 
public static String scale(String label, 
double maximum, 
double minimum, 
double x, 
double y, 
double z, 
double xAxis, 
double yAxis, 
double zAxıs, 
double theta, 
int divisions) { 
int 170; 
String value = ""; 
String ts = ""; 
ts +="Transform {"+appendage; 
ts +="translation "+x+" "+y+" "+z+appendage, 
ts +="rotation "+xAxis+" "+yAxis+" "+zAxis+" "+theta+appendage; 
ts +="children ["+appendage; 
ts +="Shape {"+appendage; 
ts +="appearance Appearance {"+appendage; 
ts +="material Material {"+appendage; 
ts +="emissiveColor 1 1 1"+appendage, 
ts +="y"+appendage, 
ts +="}"+appendage; 
ts +="geometry IndexedLineSet {"+appendage; 
ts +="coord Coordinate {"+appendage; 
ts +="point ["+appendage; 
ts +=""+minimum+"” O 0"+appendage, 
ts +=""+maximum>+" O 0"+appendage, 
ts +="]"+appendage, 
ts +="}"+appendage; 


ts +="coordIndex ["+appendage; 
ts +="0 1 -1"+appendage; 

ts +="]"+appendage; 

ts +="}"+appendage; 

ts +="3"+appendage; 


for (1=0; i<=divisions,itt) { 
if (minimum > maximum) { 
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value = ""+ int) maxımum+ı*((minimum-maxımum)/divisions)); 
j 
else { 

value = ""+(intXminimum+H*((maximum-minimum)/divisions)); 
j 


ts *-textNode( VERTICAL, 
minimum-(maximum-minimum)*(double)/(double)divisions, 
0, 
0, 
value, 
Math.abs(((maximum-minimum 4.0 Y divisions), 
Math.abs(((maximum-minimum)/4.0)/ divisions), 
RIGHT); 


j 
ts +=textNode( HORIZONTAL, 

(minimum+maximum)/2.0, 
-Math.abs((minimum-maximum)/divisions), 
0, 
label, 
Math.abs(((maximum-minimumy/2.0 y divisions), 
Math.abs(((maximum-minimum)/2.0)/divisions), 
CENTER); 


ts +="]"+appendage; 
ts +="}"+appendage; 
return ts; 

j 
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= Private Section 
a بر ابر بای اب ما ی ی ایو سب اب‎ ERE ECE E | 
private static String textNode(String orientation, 
double x, 
double y, 
double z, 
String text, 
double size, 
double space, 
String justify) ( 


String ts = "Transform ("-appendage; 
ts += "translation "+x+" "+y+" "+z+appendage, 
ts += "rotation "torientationtappendage; 
ts += "children Shape {"+appendage; 
ts += "appearance Appearance {material Maternal { emissiveColor 1 1 133"+appendage; 
ts += "geometry Text {"+appendage; 
ts += "string \" "+text+" \" "tappendage; 
ts += "fontStyle FontStyle {"+appendage; 
ts += "size "+size+appendage; 
ts += "spacing "+space+appendage; 
ts += "justify "+justify+appendage, 
ts += "3 "+appendage; 
ts += "3 "+appendage, 
ts += "}}"+appendage; 
return ts; 
j 


private static final String LINE FEED 7 "in^; 
private static final String SPACE =" "; 
private static boolean appendLineFeed = true; 
private static String appendage = LINE_FEED, 
private static int colorScheme = 0; 
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private static double colorMaximum = 0.0; 
private static double colorMinimum = 0.0; 
private static double colorDelta = 0.0; 

private static double colorSubDelta = 0.0; 
private static double colorDetection = 0.0; 
private static double colorCounterDetection = 0.0; 
private static double color Value = 0.0; 

private static int intensityScheme = 0; 

private static double intensityMaximum = 0.0; 
private static double intensityMinimum = 0.0; 
private static double intensityDelta = 0.0; 
private static double intensity Value = 0.0; 


E. RAY.JAVA 


/* 

‚File: Ray.java 
Compiler: jdk1.1.3 
b 


package mil.navy.nps.rra; 


import mil.navy.nps.rra.*; 
import java.lang.Math; 


q 


*(g version 1.0 

*@author LT Timothy M. Holliday (<A 
HREF="http://dubhe.cc.nps.navy.mil/~tmhollid">http://dubhe.cc.nps.navy.mil/~tmhollid</A>) 
* e 

*<dt><b>Location:</b> 

*<dd><a href="http://www.stl.nps.navy.mildis-java-vrml/mil/navy/nps/rra/Ray.java"> 

* http//www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/Ray.java«/a» 

x 


*<dt><b>Hierarchy Diagram:</b> 
*<dd><a href="images/RRAC lassHierarchy.gif"><IMG SRC="1mages/RRAClassHierarchyButton. jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Coordinate System Diagram:</b> 
*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/CoordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Summary:</b> 

*<dd>This class incorporates the RRA (recursive ray acoustics) 
* algorithm from Professor Ziomek. 

* 


*<dt><b>Explanation:</b> 

*<dd>Any wave front can be decomposed into a set of plane waves 
* with the plane waves travelling in the direction of the normal 

to the wave front. The normal for each plane wave is defined by 
the normal to the wave front at the point where the front 

and the plane wave touch. This class uses the ray theory solution 
to the linear acoustic wave equation to simulate one ray from 

the surface of a wavefront as it propagetes through the ocean.<P> 


$ $ Y + $ a 


*<dt><b>History:</b> 

*<dd> 1Nov97 /Timothy M. Holliday /New 

*<dd> 17Mar98 /Timothy M. Holliday /Added HTML comment convention 
*<dd> 12Apr98 /Timothy M. Holliday /Parameterless Constructors 
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*<dd> 14Apr08 /Timothy M. Holliday /Simplified VRML Routines 


% 


*@see Surface 
*@see Bottom 
*@see Beam 
*@see ExampleRay 
*@see Vec3d 
*/ 

public class Ray { 


/** 
* Constructor for the Ray class. Constructs a fixed size structure to trace a 
* ray of acoustic energy through the ocean. 
t/ 
public Ray () { 
/* 


* [nitialize the member variables that require it 
= 

position = new Vec3d(); 

normal = new Vec3d(); 

startPosition = new Vec3d(); 

prevPosition = new Vec3d(); 

prevNormal = new Vec3d(); 

deltaPosition = new Vec3d(); 

starPosition = new Vec3d(); 

trailingPosition = new Vec3d(), 


int 1; 

for (=0;1<MAX_ POINTS;i++) { 
rayPath[i] = new Vec3d(), 
trailingPath[i] = new Vec3d(); 


j 


fer 

* This method resets all of the ray parameters after 

* instanciation has occurred since reuse is more time 
* efficient than garbabage collection and reallocation. 
2 

public void reset() { 


per 

* [nitialize the member variables that require it 
2 

prevPosition.set(position); 
setNormal(); 
prevNormal.set(normal); 
deltaPosition.set(position); 
starPosition.set(position ), 
trailingPosition.set(position); 
time = 0.0; 

pathLength - 0.0; 
reflectionPhase = 0.0; 
pointCount = 0; 

reflection = false; 
absorbtion = 0.0; 
totalCurve = 0.0; 


j 


er 
* This method takes the direction cosines of the ray 
* and forms the vector normal to the wavefront. 
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E 
private void setNormal() 1 
double u,v,w; 


u = (double)(Math.sin(elevation*Math.PV/180) * Math.cos(azımuth*Math.PV/180)); 
v = (doubleXMath.cos(elevation*Math. PI/180)); 
w 7 (double)( Math.sin(elevation* Math.PI/180) * Math.sin(azimuth* Math.PI/180)); 


normal.set(u,v,w); 
} 


/ 2 ٭‎ 
* This method returns the normal vector to the wavefront. 
*/ 
public Vec3d getNormal() { 
return normal; 
j 


= 
* This method sets the azimuthal angle, which is the angle from the x-axis 
* to the z-axis rotating about the y-axis. 
El 
public void setAzimuth(double phi) { 
azımuth = phı; 
j 


7/۰ 
* This method returns the azimuthal angle. 
2 
public double getAzimuth() 1 
return azimuth; 


j 


per 
* This method sets the elevation angle, which is the angle from the y-axis 
* to the x-axis rotating about the z-axis. 
= 
public void setElevation(double beta) { 
elevation = beta; 
} 


per 
* This method returns the elevation angle. 
D 
public double getElevation() { 
return elevation; 
j 


ner 
* This method sets the position of the ray. 
el 
public void setPosition(double x, double y, double > { 
startPosition.set(X,y,Z), 
position.set(x,y,Z); 
prevPosition.set(x,y,Z); 
time = deltaTime; 
} 


/ Kk 
* This method gets the position of the ray. 
El 
public Vec3d getPosition() { 
return position; 
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j 


/ xx 
* This method is used to change the time step of the ray after 
* instanciation. 
* 
public void setDeltaTime(double deltaTime) { 
time = deltaTime; 
j 


/* * 
* This method is returns the time step of the ray. 
el 
public double getDeltaTımef) { 
Teturn time; 
} 


je * 
* This method sets the ping duration. Duration is currently 1 or 2. 
* The integer refers to the number of deltaTime increments. 
*/ 
public void setDuration(int duration) { 
pulseDuration = duration; 
j 


pex 
* This returns the ping duration. Duration is currently 1 or 2. 
* The integer refers to the number of deltTime increments. 
E 
public int getDuration() f 

return pulseDuration; 
j 


/* * 
* This method sets the handle to the bottom object. 
x 
public void setBottom(Bottom pBottom) { 
bottom = pBottom; 
} 


[** 
* This method returns the handle to the bottom object. 
El 
public Bottom getBottom() { 
return bottom; 
j 


er 
* This method sets the handle to the surface object. 
*/ 
public void setSurface(Surface pSurface) { 
surface = pSurface; 
j 


/* * 
* This method returns the handle to the surface object. 
*/ 
public Surface getSurface() { 
retum surface; 
j 


/ ax 
* This method sets the handle to the sound speed profile object. 
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xi 

public void setSsp(SSP pSsp) { 
ssp = pSsp; 

j 


q, 
* This method returns the handle to the sound speed profile object. 
t) 
public SSP getSsp() 1 
return sSp; 
j 


per 
* This method returns the phase change of the wavefront 
a 
public double getReflectionPhase() { 
return reflectionPhase; 
j 


pur 
* This method retums the x-component of the position 
* of the wavefront 
*/ 
public double getPositionX() { 
return position.get(0); 
j 


pe 
* This method retums the y-component of the position 
* ofthe wavefront 
*/ 
public double getPositionY() ( 
return position. get(1); 
j 


e 
* This method returns the z-component of the position 
* ofthe wavefront 
«f 
public double getPositionZ() { 
return position.get(2), 
j 


per 
* This method returns the position of the trailing edge of 
* the wavefront 
d 
public Vec3d getTrailingPosition() { 
return trailingPosition; 
j 


/** 

* This method returns the current simulation time. 

*/ 

public double getTime() { 

return time; 

} 
per 

* This method returns a saved simulation time. 

e] 
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public double getTime(int index) { 
return ray l'ime[index]; 
j 


q 
* This method causes the ray to propagate one timestep into the future. 
= 
public void Propagate(double timeStep) { 
deltaTime = timeStep; 
nextPosition( ); 
nextNormal(); 
next TimeandpathLength(); 


} 


193 
* "This method retums the current number of points stored for the ray. 
4 
public int getCount() { 
return pointCount, 
۱ 
f 


[** 
* This method returns the current total relaxation absorbtion in dB. 
*/ 
public double getAbsorption() { 

return absorbtion; 
j 


per 
* This method returns whether or not the ray has been reflected in 
* the most recent time step. 
= 
public boolean reflected() { 
return reflection; 


j 


/** 
* This method forces a recording of the current ray position. It also forces the 
* total curvature of the ray since the last recorded point to 0.0. 
public void recordPoint() { 
rayPath{pointCount].set(position), 


if (pulseDuration == 2) { 
trailingPath(pointCount].set(trailingPosition); 
j 


else ( 
trailingPath{pointCount].set(prevPosition); 


j 


rayTime[pointCount] = time; 


pointCount++; 
totalCurve = 0.0; 
} 
[== 


* This returns the total curvature of the ray path since the last recorded point. 


* 


* The curvature of a small segment of curve can be approximated by the curvature 
* of a circle which is the change in angular position around the circle divided 

* by the pathlength change. K = deltaTheta/deltaPathLength. Thus total curvature 
* is the sum of the curvatures of the small segments along a ray trajectory. 
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*/ 
public double totalCurvature() 1 
return totalCurve; 


j 


per 

* This method returns the requested position as a String. 

4 

public String position(int N) { 

return (int)ravPath[N].get(0)  " " + 

(intJrayPath[N].get(1)+""+ 
(int)rayPath[N].get(2) + appendage; 

} 


ye. 
* This method returns the requested trailing edge of the ray as a String. 
8 
public String trailingPosition(int N) { 
return (int)trailingPath[N].get(0) 4 " " + 
(int)trailingPath[N].get(1) + + 
(int)trailingPath[N].get(2) + appendage; 
j 


yum 
* This method returns the requested normalized time as a string. 
E 
public String normalizedTime(int N, double endtime) { 

return Math.floor( 100*rayTime[N l/endtime)/100 + appendage; 
j 


jus 
* This is a static method used to indicate whether a line feed 1s desired 
* at the end of every line. 
* 
* 'true' indicates a linefeed 1s desired and ‘false’ indicates that a space 
* is desired" 
= 
public static void setAppendLineFeed(boolean pAppendLineFeed) { 
appendLineFeed = pAppendLineFeed; 
if (appendLineFeed) { 
appendage = LINE FEED; 


else { 
appendage = SPACE; 
} 
} 


[** 
* This is a static method that returns the current line appendage. 
x 
public static boolean getAppendLineFeed(){ 
retum appendLineFeed; 
j 
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7 Beginning of private section 
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private static final String LINE FEED = "\n"; 
private static final String SPACE - " "; 
private static boolean appendLineFeed = true; 
private static String appendage - LINE FEED; 
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/* 

* This method calculates the position of the ray after the time step. 
2/ 

private void nextPosition() { 


// save trailing position depending on duration 
if (pulseDuration == 2) { 
trailingPosition.set(prevPosition), 
i 
else { 
trailingPosition.set(position); 


} 


prevPosition.set(position); 
prevdeltaArc = deltaArc; 
prevNormal.set(normal); 
position.add(deltaPosition()); 
reflection = false; 


if (bottom.intersected(position)) { 
reflectionPhase += bottom.reflect(position, prevNormal); 
reflection = true; 


j 


if (surface.intersected(position)) { 
reflectionPhase += surface.reflect(position,prevNormal), 
reflection = true; 
j 
j 


/* 
* This method calculates the normal to the wavefront at the position of the ray 
* after the time step using the ziomek formula. 
2 
private void nextNormal() { 
updatestarPosition(); 
deltaArc = deltaTime*ssp.C(position); 
termi ..scale(ssp.C(position)/ssp.C(prevPosition),prevNormal); 
term2.scale(deltaTime,ssp.gradC(starPosition)); 
totalCurve += term2.length(); 
normal.sub(prevNormal,term2); 


j 


/* 
* This method the calcultes the total elapsed time and length of path travelled. 
*/ 
private void nextTimeandpathLength() ( 
time += deltaTime; 
pathLength += deltaTime*ssp.C(position); 
absorbtion = pathLength * 1e-3; 
j 


[* 
* This method calcultes the vector that expresses the change from current position 
* to the next. 
21 
private Vec3d deltaPosition() { 
temp.scale(prevdeltaArc, prevNormal); 
retum temp; 
1 
f 


/* 
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* This method calculates the midpoint between the two most recent points ın 
* the ray path. 
*/ 
private void updatestarPosition() { 
temp.scale(prevdeltaArc/2.0.prevNormal); 
starPosition.add(prevPosition,temp); 


j 


/* 

* This method sets the x position of the ray 

*/ 

private void setPositionX(double x) ( 
position.set(0,x); 


/* 

* This method sets the y position of the ray 

er 

private void setPositionY(double y) { 
position.set(1,y); 

j 


/* 

* This method sets the z position of the ray 

E 

private void setPositionZ(double z) ( 
position.set(2,z); 

j 


/* 

* Maximum number of points of of information along the ray 
* path. 

a 

private static final int MAX POINTS = 100; 


/* 

* Handle to the sound speed profile object 
t 

private SSP ssp= null; 

/* 

* Handle to the bottom object 

© 

private Bottom bottom = null; 
/* 

* Handle to the surface object 

= 

private Surface surface = null; 
7* 

* Starting position of the ray 

i 


private Vec3d startPosition= null; 


/* 

* Position of the ray after the last iteration 
E 

private Vec3d  prevPosition - null; 


y 
* Position ofthe ray in this iteration , 
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7, 


private Vec3d position = null; 
/۴ 
* Change in position in this iteration 
2 


private Vec3d deltaPosition= null; 


[* 

* Position halfway between position and 
* position+deltaPosition 

Ej 

private Vec3d starPosition= null; 


/* 

* Position of the trailing edge of the ray segment, not 
* necessarily the prevPosition 

a 

private Vec3d trailingPosition = null; 


Je 

* Direction of ray propagation in the last iteration 
3) 

private Vec3d prevNormal = null; 


/۴ 

* Direction of ray propagation in the current iteration 
*, 

private Vec3d normal = null; 


/* 

* Change in path length in the current iteration 
*/ 

private double deltaArc = 2.0; 


[Z 

* Change in path length in the previous iteration 
1) 

private double prevdeltaArc - 0.0; 


/* 

* Simulation time for the ray 

t 

private double time = 0.0; 


j* 

* Number of points saved for later use 
Ej 

private int pointCount = 0; 


/* 

* Total curavature of the ray path, used to 

* decide when a point should automatically be saved 
2 

private double totalCurve = 0.0; 


/* 

* Indicates if a surface or bottom reflection has occurred 
€] 

private boolean reflection = false; 


/* 
* Angular elevation of the ray from downward (Y-direction) 
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* toward the (X-Z plane). In the RRA coordinate system 
el 
private double elevation = 85.0; 


/ * 

* Angular displacement of the ray in the (X-Z plane) 
* about the (^ Y-axis) 

f 

private double azimuth = 0.0; 


ys 

* Change in time in one iteration 

7 

private double deltaTime = 0.006; 


/ M 

* Number of deltaTime segments from the beginning 
* ofthe ray segment to the trailing edge 

*/ 

private int  pulseDuration= 0; 


/* 
* Total relaxation absorbtion along the ray path 
E 


private double absorbtion = 0.0; 

pe 

* Total path length that the ray traverses 
E] 

private double pathLength = 0.0; 

/* 

* Total phase changes due to reflection 
= 


private double reflectionPhase = 0.0; 


f * 

* Storage of the positions of the ray at the desired points 

* along the ray path 

21 

private Vec3d[] rayPath= new Vec3d[MAX_ POINTS]; 


/ * 

* Storage of the trailing edge positions at the desired 
* points along the ray path 

ا 

private Vec3d[] trailingPath ? new Vec3d[MAX, POINTS]; 


pe 

* Storage of simulation time at the desired points along 

* the ray path 

e] 

private double[] rayTime = new double[MAX: POINTS]; 


۰ر7 

* Variable allocated at instanciation to speed execution 
* so that we do not have to wait for dynamic allocation 
* at every iteration 

*/ 

private Vec3d gradc = new Vec3d(); 


/* 
* Temporary variable allocated at instanciation so that 
135 


* time is not wasted allocating 1t dynamically. 
*/ 
private Vec3d termi = new Vec3d(); 


pe 

* Temporary variable allocated at instanciation so that 
* time is not wasted allocating it dynamically. 

*/ 

private Vec3d term2 = new Vec3d(); 


/* 

* Temporary variable allocated at instanciation so that 
* time is not wasted allocating it dynamically. 

*/ 

private Vec3d temp = new Vec3d(); 


/۴ 

* Temporary variable allocated at instanciation so that 
* time is not wasted allocating it dynamically. 

27 

private double temp3; 


F. SSP.JAVA 


pe 
File: SSP.java 
Compiler: jdk1.1.6 
E 
package mil.navy.nps.rra; 


import mil.navy.nps.rra.*; 


per 

*@version 1.0 

*@author LT Timothy M. Holliday (<A HREF="http://www.stl.nps.navy.mil/~auv/holliday"> 
http://www. sti. nps.navy.mil/~auv/holliday</A>) 

x 


*<dt><b>Location:</b> 

*<dd><a href="http://www.stl_nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/SSP.java"> 

* http//www.stl.nps.navy.mildis-java-vrml/mil/navy/nps/rra/SSP.java«/a» 

de 

*<dt><b>Hierarchy Diagram:</b> 

*<dd><a href="images/RRAClassHierarchy.gif"><IMG SRC="images/RRAClassHierarchyButton.jpg" 
ALIGN=ABSCENTER></a> 

* 

*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem. gif"><IMG SRC="images/CoordinateSystemButton. jpg" 
ALIGN=ABSCENTER></a> 

ae 


*<dt><b>Summary:</b> 

*«dd»* This class simulates the characteristics of the sound speed 
* ofthe volume of ocean under simulation. 

+ 

*<dt><b>Explanation:</b> 

*<dd> This class simulates five standard profiles: 

* 1) no gradient - constant sound velocity of 1500 m/s 
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2) positive gradient - constantly increasing sound velocity 
starting with 1500 m/s and having a slope of .017 m/s/m 
3) negative gradient - constantly decreasing sound velocity 
starting with 1500 m/s and having a slope of .017 m/s/m 
4) parabolic - velocity profile obeying 
1490 + 4e-5 *(500 - depth)”2 
5) traditional - velocity profile obeying 
at Om 1500 m/s 
Om- 100m + slope of .016 m/s/m 
1000m - 2000m - slope of .02956 m/s/m 
> 2000m + slope of .017 m/s/m<P> 
this profile froms the typical deep sound channel 
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*<dt><b>History:</b> 
*<dd> 250ct97 /Timothy M. Holliday [New 
*<dd> 17Mar98 /Timothy M. Holliday /Added HTML comment convention 
* 
*/ 
public class SSP { 


۰ر 


* Constructor for the SSP class. Its argument is one the the following 
* five sound speed profile types: 'constant', '+gradient’, '-gradient', 
* 'parabolic' and 'traditional' 
*/ 
public SSP (String sspType) { 
if (sspType.equalsignoreCase("constant”)) { 
environment 7 1; 


j 
else if (sspType.equalslgnoreCase("+gradient")) { 
environment = 2; 


} 
else if (sspT ype.equalsIgnoreCase("-gradient")) { 
environment = 3; 


} 

else if (sspType.equalsIgnoreCase("parabolic")) ( 
environment = 4; 

} 

else if (sspType.equalsIgnoreCase( "traditional")) f 
environment = 5; 

j 

else { 
System.out.printin("User defined not implemented yet, defaulting to constant"); 
environment = 1; 

} 

} 


[** 

* Method used to calculate the gradient of sound speed for each of the 
* five standard profiles 

71 

public Vec3d gradC(Vec3d starPosition) { 


switch (environment) { 
case 1: gradient - gradC l(starPosition); break; 
case 2: gradient = gradC2(starPosition); break; 
case 3: gradient = gradC3(starPosition); break; 
case 4: gradient = gradC4(starPosition);, break; 
case 5: gradient = gradC5(starPosition);, break; 

} 

return gradient; 

} 
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pex 
* Method used to calculate the gradient of sound speed from any sound 
* speed profile. 

*/ 

public Vec3d generalgradC(Vec3d starPosition) { 


position2.add(starPosition deltax); 
position1.sub(starPosition,deltax); 
xgradient = (C(position2)-C(position! ))/delta; 
position2.add(starPosition,deltay); 
position1.sub(starPosition,deltay); 
ygradient = (C(position2 }-C(position] ))/delta; 
position2.add(starPosition,deltaz); 
position] .sub(starPosition,deltaz); 
zgradient 7 (C(position2 -C(positionl ))/delta; 


gradient.set(xgradient,ygradient,zgradıent); 


return gradient, 
j 


= 

* Method used to calculate the speed of sound for each of the five 
* standard profiles. 

t 

public double C( Vec3d position) { 


switch (environment) { 
case 1: soundSpeed = SSP (position); break; 
case 2: soundSpeed = SSP2(position); break; 
case 3: soundSpeed = SSP3(position), break, 
case 4: soundSpeed = SSP4(position); break; 
case 5: soundSpeed = SSP5(position);, break; 

} 

return soundSpeed; 


} 


ا کا ا ا ا ا کا ا ا ا ا ا ا A‏ ا ا ا A A‏ 
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** Private setion begins here 


xx 


29 ۳ ۳ ۳ ۳ EECEASE SE RESALE LESL SLES ES ES REESE SEES ESSE AEST EH / 
private double SSP1(Vec3d position) { 
return 1500.0; 
j 


private double SSP2(Vec3d position) ( 
return (1500.0+0.017*position.get(1)); 
j 


private double SSP3(Vec3d position) ( 
return (1500.0-0.017*position.get(1)); 
j 


private double SSP4(Vec3d position) { 
retum (1490.0 + 4.0e-5 * (500.0 - position.get(1))*(500.0 - position.get(1))); 
j 


private double SSPS(Vec3d position) { 


if ( position.get(1) <= 100) { 
soundSpeed = 1500.0+0.016*position. get(1 ); 
j 
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else if ( position.get(1) <= 1000) { 
soundSpeed = 1501.6 - 0.02956 * (position.get(1) - 100); 
j 


else 
soundSpeed = 1475.0 + 0.03 * (position. get(1) -1000), 


retum soundSpeed; 
j 


private Vec3d gradC1( Vec3d position) { 
gradient.set(0,0,0); 
return gradient; 

j 


private Vec3d gradC2(Vec3d position) { 
gradient. set(0,7-0.01 7,0); 
return gradient; 

} 


private Vec3d gradC3(Vec3d position) { 
gradient.set(0,-0.017,0), 
return gradient; 

j 


private Vec3d gradC4( Vec3d position) { 
gradient.set(0,4.0e-5 * -2 *(500.0 - position.get(1)),0); 
return gradient; 

j 


private Vec3d gradC5(Vec3d position) { 


if ( position.get(1) <= 100) { 
gradient.set(0,0.016,0); 


j 

else if ( position.get(1) <= 1000) { 
gradient. set(0,-0.02956,0); 

j 

else 

{ 
gradient.set(0,0.017,0); 

j 

return gradient; 

j 


/** standard environment that was initialized */ 
private int environment = 1; 


/** class constant defining incrementally small variation in x direction*/ 

private final Vec3d deltax = new Vec3d(0.01, 0.0, 0.0); 

/** class constant defining incrementally small variation in y direction*/ 

private final Vec3d deltay = new Vec3d(0.0, 0.01, 0.0); 

/** class constant defining incrementally small variation in z direction*/ 

private final Vec3d deltaz = new Vec3d(0.0, 0.0, 0.01); 

/** class constant defining incrementally small variation in magnitude in any direction*/ 
private final double delta = 0.02; 


/** class variable used in calculating speed of sound */ 
private double soundSpeed; 


/** class variable used in calculating the gradient of sound speed*/ 
private Vec3d position] = new Vec3d(); 
/** class variable used in calculating the gradient of sound speed*/ 
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private Vec3d position2 = new Vec3d(); 

/** class variable used in calculating the gradient of sound speed*/ 
private Vec3d gradient = new Vec3d(); 

/** class variable used in calculating the gradient of sound speed */ 
private double xgradient; 

/** class variable used in calculating the gradient of sound speed*/ 
private double ygradient, 

/** class variable used in calculating the gradient of sound speed*/ 
private double zgradient; 


G. SURFACE.JAVA 


/* 

File: Surface.java 
Compiler: jdk1.1.6 
e 

package mil.navy.nps.rra; 


import mil.navy.nps.rra.*; 


[ux 

*@version 1.0 

*@author LT Timothy M. Holliday (<A 
HREF="http://www.stl.nps.navy.mil/~auv/holliday">http://www.stl.nps.navy.mil/~auv/holliday</A>) 
x 


*<dt><b>Location:</b> 

*<dd><a href="http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/tra/Surface.java"> 

* http://www. stl.nps.navy.miV/dis-java-vrml/miVnavy/nps/rra/Surface.java</a> 

x 

*<dt><b>Hierarchy Diagram:</b> 

*<dd><a href="1mages/RRAClassHierarchy. gif"><IMG SRC="images/RRAClassHierarchyButton. jpg" 
ALIGN=ABSCENTER></a> 

ak 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/CoordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

E 


*<dt><b>Summary:</b> 

*<dd>* This class acoustically simulates the characteristics of 

* the ocean's surface. It is simulated as being a pressure release 

* surface with no roughness 

* 

*<dt><b>Explanatıon:</b> 

*<dd>Since a ray is the normal vector to a plane wave, rays interact with 
surfaces the same way that plane waves do. Rays obey Snell's law and 
the equations of reflection and transmission of wave energy. This model 
uses vector algebra to determine when the end of a ray has passed through 
the bottom and to reflect any rays that have penetrated. Since the 
surface is being approximated, to a good degree of accuracy, as a perfect 
pressure release surface any wave striking it will have a 180 degree 
phase shift and a reflection coefficient of one.<P> 
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*<dt><b>History:</b> 
*<dd> 8Nov97 /Timothy M. Holliday /New 
*<dd> 18Mar98 /Timothy M. Holliday /Added HTML comment convention 
x 
*@see Ray 
*@see Bottom 
*(@see Vec3d 
2 
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public class Surface { 


I HKH 
* Constructor for the surface class. Its argument has one choice, smooth. 
* Smooth ındıcates a perfect pressure release surface. 
El 
public Surface (String surfaceType) { 
if (surface Type.equalsIgnoreCase("smooth")) f 
environment = 1; 
j 
else { 
System.out.printin("User defined not implemented yet, defaulting to smooth"); 
environment = 1; 
} 
j 


[** 

* This procedure checks to see 1f the ray path has crossed the surface 
* during the current time step and returns true if it did. 

E, 

public boolean intersected( Vec3d Pos) ( 


intersected = false; 
switch (environment) { 
case 1: 
if (Pos.get(1) < surfaceDepth) { 
intersected = true; 
} 
break; 


return intersected; 


j 


[EZ 
* This method causes a Snell's Law reflection to occur at the surface of the ocean. 
D 
public double reflect( Vec3d Pos, Vec3d normal) { 
switch (environment) { 
case 1: 
double diff, 
diff = Pos.get(1) - surfaceDepth; 
Pos.set( 1 ,surfaceDepth-dift); 
normal.set(1 ,-normal.get( 1)); 


break; 


j 
return Math.PI; 
j 


per 
* This method returns a VRML string representing the Surface of the ocean. 
e 

public String VRMLSurface() ( 


int minX = -5000; 
int minZ = -5000; 
int maxX = 5000; 
int maxZ = 5000; 


return "Transform {" + appendage + 
" rotation 1 0 0 3.14" + appendage + 
" children" + appendage + 


"Shape {" + appendage + 
" appearance Appearance {" + appendage + 
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" material Material {" + appendage + 

emissiveColor 0 0 .4" + appendage + 

" diffuseColor 0 0 0" + appendage + 

"  transparency .9" + appendage + 

}" + appendage + 

}" + appendage + 
" geometry IndexedFaceSet {" + appendage + 
" solid FALSE" + appendage + 

coord Coordinate {" + appendage + 

" point [" + appendage + 
" "+ maxX +" "+ surfaceDepth +" "+ maxZ +"," + appendage + 
5 "+ maxX +" "+ surfaceDepth +" "+ minZ +"," + appendage + 
z "+ minX +" "+ surfaceDepth +" "+ minZ +"," + appendage + 
: "+ minX +" "+ surfaceDepth +" "+ maxZ +"" + appendage + 
" ]" + appendage + 
" }" + appendage + 
۲ coordindex [ 0,1,2,3 ]" + appendage + 
" 1 + appendage + 
"y" + appendage + 
"3" + appendage; 

j 


pur 
* This 1s a static method used to indicate whether a line feed 1s desired 
* atthe end of every line. 


* 
* ‘true’ indicates a linefeed is desired and 'false' indicates that a space 
* 1s desired" 
* 
public static void setAppendLineFeed(boolean pAppendLineFeed) { 

appendLineFeed = pAppendLineFeed; 

if (appendLineFeed) { 

appendage = LINE_FEED, 


else { 
appendage = SPACE; 


j 


[te 
* This 1s a static method that returns the current line appendage. 
e 
public static boolean getAppendLineFeed() { 
return appendLineFeed; 
j 
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ar Private Section 
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private static final String LINE FEED - "n^; 
private static final String SPACE =" "; 
private static boolean appendLineFeed = true; 
private static String appendage » LINE FEED; 


/** standard environment that was initialized */ 
private int environment = 1; 


/** seems useless, but allows for variations due to sea state*/ 
private double surfaceDepth = 0; 


/** temporary boolean variable indicating if a reflection has occured*/ 
private boolean intersected; 
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H.  TARGETSJAVA- 


/* 
File: Targets.java 
Compiler: jdk1.1.6 
= 


package mil.navy.nps.rra; 
import mil.navy.nps.rra.*; 


/** 


*@version 1.0 

*@author LT Timothy M. Holliday (<A HREF="www.stl.nps.navy.mil/—auv/holliday"> 
http://www.stl.nps.navy.mil/~auv/holliday</A>) 

* 


*<dt><b>Location:</b> 

*<dd><a href="http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/Targets.java"> 

* http://www.stl.nps.navy.mil/dis-¡ava-vrml/mil/navy/nps/rra/Targets.java</a> 

* 

*<dt><b>Hierarchy Diagram:</b> 

*<dd><a href="1mages/RRAClassHierarchy.gif"><IMG SRC="images/RRAClassHierarchyButton.jpe" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/CoordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Summary:</b> 

*<dd>This class tracks the positions of targets in the VRML scene. 
* [t provides for random target generation and a method to see if 

* asonar ping has collided with a target. 

ik 

*<dt><b>Explanation:</b> 

*<dd>As acoustic waves traverse the ocean environment they impinge 
on objects. The bottom and suface are constraints on the water 
colomn extent and as such they are handled in special classes 

by themselves. All other objects, those ın the water column, 

are considered targets. When the sonar hits the target, energy 

is scattered off the target in a mainly isotropic manner. This 
scattered energy 1s what returns to the sending sonar giving 

away the location of the target. In the simulation the plane 
forming the front of the wave front is compared to the distance 

of the target awar from the plane by means of a vector 

algebra clculation. When this distance changes sign, there 

is a possibility that the sonar pulse hit the object, and 

continued on its way. Further checks are then done to verify 
whether or not a detection event has occurred. The corresponding 
information is then made available to the calling method as 

the return value of a method. This class also provides a method 
for randomly generating sample targets and is the logical place 

to add code to receive position information from external sources, 
like the internet via DIS.<P> 
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*<dt><b>History:</b> 

*<dd> 30Jan98 /Timothy M. Holliday /New 

*<dd> 17Mar98 /Timothy M. Holliday /Added HTML comment convention 
*<dd> 13May98 /Timothy M. Holliday /New Constructor for placing a 

P minefield. New statistic method. 
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* 


*@see Beam 
=) 
public class Targets { 


×۴ 


* Constructor for the targets class. The first two arguments supply the surface 
* and bottom types. The second two indicate how many mines and submarines to add 
* tothe target data base. The minestsubmarines are limited to 100 and they are 
* randomly placed. 
El 
public Targets (Surface s,Bottom b, int numberMines, int numberSubs) { 
int counter, 
double x,y,z; 


skipCount = 0; 

surfaceNormal = new Vec3d(); 
position = new Vec3d(); 
segment] = new Vec3d(); 
segment2 = new Vec3d(); 


numberofMines = numberMines; 
numberofSubs = numberSubs; 


if ((numberofMinest+numberofSubs >=MAX TARGETS) { 
System.out.pnntln("Too many Targets: " + (numberofMinestnumberofSubs)); 


else { 

for (counter=0;counter<numberofMines/2;counter++) { 
x 7 Math.random()*2500; //in ziomek coordinate system 
z - Math.random()*48004-200; 
y = Math.random()*b.depth(x.z); 
target[counter] » new Vec3d(x,y,z); 
targetSize[counter][0] 7 1.0; // x dimension 
targetSize[counter][1] = 1.0; // y dimension 
targetSize[counter][2] = 1.0; // z dimension 
targetArea[counter] = 3.1415; 
targetDistance[counter] = 0.0; 

j 

for (,counter<numberofMines;counter++) { 
x = Math.random()*2500; //in ziomek coordinate system 
z = Math.random()*-4800-200; 
y = Math.random()*b.depth(x,z); 
target[counter] = new Vec3d(x,y,z); 
targetSize[counter][0] = 1.0; // x dimension 
targetSize[counter][1] = 1.0; // y dimension 
targetSize[counter][2] = 1.0; //z dimension 
targetArea[ counter] = 3.1415; 
targetDistance[counter] = 0.0; 


for (;counter<numberofSubst+numberofMines;counter++) { 
x = Math.random()*-1000; 
z = Math.random()*800-400; 
y = Math.random()*b.depth(x,z); 
target[counter] = new Vec3d(x,y,z); 
targetSize[counter][0] = 100.0; // x dimension 
targetSize[counter][1] = 10.0; // y dimension 
targetSize[counter][2] = 10.0; // z dimension 
targetArea[counter] = 3.1415*100; 
targetDistance[counter] = 0.0; 
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[** 
* Constructor for the tartgets class. The first two arguments are 
* the bottom and surface objects. The next two are the max and min 
* x-direction boundaries for a minefield and the next two are for 
* the z-direction. The last is the number of mines. The maaximum 
* number of mines is 100 and they are randomly placed. 
El 
public Targets (Surface s, 
Bottom b, 
double maximumX, 
double minimumX, 
double maximumZ, 
double minimumZ, 
int numberMines) ( 
int counter, 
double x,y,z; 


skipCount = 0; 

surfaceNormal = new Vec3d(); 
position = new Vec3d(); 
segment] = new Vec3d(); 
segment? = new Vec3d(), 


numberofMines = numberMines, 
numberofSubs = 0; 


if ((numberofMines+numberofSubs)>=MAX TARGETS) { 
System.out.println("Too many Targets: " + (numberofMines+numberofSubs)); 
j 
else { 
for (counter=0;counter<numberofMines;counter++) { 
x = Math.random()*(maximumX-minimumX)+minimum;, //in ziomek coordinate system 
z = Math.random()*(maximumZ-minimumZ }hminimumZ,; 
y = Math. random()*b.depth(x,z), 
target[counter] = new Vec3d(x,y,z); 
targetSize[counter][0] 7 1.0; // x dimension 
targetSıze[counter][1] = 1.0; // y dimension 
targetSize[counter][2] 7 1.0; // z dimension 
targetArea[counter] = 3.1415; 
targetDistance[counter] = 0.0; 
numberOfDetects[counter] = 0; 
estimatedPosition[counter] = new Vec3d(); 
j 
j 
j 


/ * x 
* This method resets the parameters that speed up target collision detection. 
* It is used whenever a new beam is calculated. 
*/ 
public void resetTargets() { 
skipCount = 0; ۱ 
for (counter=0;counter<numberofSubstnumberofMines;counter++) { 
targetDistance[counter] = 0.0; 
} 
} 


Jer 

* This method returns a VRML string representing the targets in the ocean. 

21 

public String VRML Targets() { 
int counter, 
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String temporaryString, 


temporaryString = "EXTERNPROTO ContactMine [\n" + 
" exposedField SF Vec3f translation ]\"ContactMine. wrl\"\n" + 


"EXTERNPROTO DieselSub [\n" + 
" exposedField SF Vec3f translation ]\"DieselSub.wrl\"\n" + 


“Transform {\n" + 
" children [\n"; 
for(counter=0;counter<numberofMines;counter++) { 
temporaryString +=" ContactMine (translation " + (int)target[counter].get(0) +" " 
+ (-(1nt)target[counter].get(1)) + " " 
t (-(1nt)target(counter].get(2))  " " 
m 
j 
for(;counter<numberofMinestnumberofSubs;counter++){f 
temporaryString+=" DieselSub {translation " + (int)target[counter].get(0) + " " 
+ (-(1nt)target(counter].get(1))  " " 
* (-(1nÜtarget[counter].get(2))  " " 


+ in? : 
j 
return temporaryString += " Jin" + 
pons 
j 
pee 


* This method returns a boolean indicating whether or not a detection as occurred. 
* True indicates a detection and false indicates no detection. 


El 

public boolean isCollision(Ray ray1, 
Ray ray2, 
Ray ray3, 
Ray ray4) ( 


collisionCount = 0; 
collision = false; 


// speed up detection by skipping distance between targets 
if (skipCount > 0) { 
skipCount—; 
j 
else 1 
//reset the minimum distance to a large value 
minDistance = 100000; 


//calculate the surface normal 

segment} .sub(ray2.getPosition(),ray1.getPosition()); 
segment2.sub(ray4.getPosition(),ray1 .getPosition()); 
surfaceNormal.cross(segment] ,segment2 ); : 
surfaceNormal.normalize( ); 

collisionArea = 0.0; 


for (counter=0;counter<(numberofMinest+numberofSubs);counter++) { 


// determine the distance of the pulse from the target 
position.sub(target[counter],ray1.getPosition()); 
distance = position.dot(surfaceNormal ); 


// if the product of the distance in the last step and the 

// current step is negative then the beam may have passed 
// through the object 

if (targetDistance[counter |*distance « 0) { 
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collision = possibleCollision(rav 1, 
ray2, 
ray3, 
ray4); 
if (collision) { 
numberOfDetects[counter]++; 
estimatedPosition[counter].add(rayl.getPosition(),ray3.getPosition()); 
estimatedPosition[counter].scale(.5); 
collisionArea += targetArea[counter]; 


else { 


j 


j 
targetDistance[counter] = distance; 


// save the distance if 1t 1s the minimum for this time step 
if (Math.abs(distance) < minDistance) { 
minDistance = Math.abs(distance), 
j 
j 


// calculate the number of time steps that looking for detections can 
// skipped from the minimum distance to the next targer,time step and 
//average sound speed 
skipCount = (intXminDistance/12.0); 
j 


return collision; 
j 


er 
* This method returns the cross sectional area of the object that was detected. 
e 
public double getCollisionArea() { 
return collisionArea; 


j 


> 
* This method prints the mine number, position, estimated position and 
* number of detects on each mine in the data base. The numbers are 
* printed to the standard console. 
2 
public void printMineStatistics() { 
int 1; 
for(1=0;1<numberofMines;1++) ( 
System. out. printin( 
"Mine "++ 
" Position "+target([1].get(0)+" "+ 
target(1].get(1)+" "+ 
target[1].get(2)+ 
" Estimate "+estimatedPosition[1].get(0)+" "+ 
estimatedPositionfi].get(1)+" "+ 
estimatedPositionfi].get(2)+ 
" Detects "+numberOfDetects[1]); 


good deed dokedojolkex doc doce ex oc ee ok oto ee ole doo eee xoc oec ooo ore 
xx 

e*t Private Section 

** 
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[* 
* This method returns a boolean indicating that a possible interaction between 
* the sonar pulse and another registered object has occurred. 


E 

private boolean possibleCollision(Ray rayl, 
Ray ray2, 
Ray ray3, 
Ray ray4) { 


if ( inBox(ray1.getTrailingPosition(),ray3.getPosition()) && 
inBox(ray2.getTrailingPosition(),ray4.getPosition()) && 
inBox(ray3.getTrailingPosition(),ray1.getPosition()) && 
inBox(ray4.getTrailingPosition(),ray2.getPosition())) { 
collisionCount++; 
return true; 

j 

else ( 

return false; 


j 
y 


/* 
* This method returns a boolean indicating whether or not the any part 
* ofthe object is within the parallelpiped formed by the sonar pulse 
M 
private boolean inBox(Vec3d corner1, Vec3d corner2) { 
if ( inInterval(corner].get(O),corner2.get(0), target[counter].get(0)) && 
inInterval(corner] .get(1),corner2.get(1), target[counter].get(1)) && 
inInterval(comerl.get(2),corner2.get(2), target[counter].get(2))) { 
return true; 
j 
else { 
return false; 
j 
j 


Je 
* This method returns a boolean indicating whether or not the any part 
* of the object is within the extent of the pulse in one of the three 
* cartesian directions. 
i ۱ 
private boolean inInterval(double pos1, double pos2, double target) { 
if ( ((pos1 « target) && (target < pos2)) || 
((pos2 < target) && (target < pos1))) ( 
return true; 
j 
else { 
return false; 
j 
j 


private static final int MAX TARGETS = 100; 

private Vec3d target{] = new Vec3d[MAX TARGETS]; 

private double targetSize[][] = new double[MAX_TARGETS][3]; 
private int numberOfDetects[] = new int{MAX_TARGETS], 

private Vec3d estimatedPosition[] = new Vec3d{MAX_TARGETS]; 


private double collisionArea = 0.0; 
private double targetArea[] = new double[MAX_TARGETS],; 
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private int numberofMines; 
private int numberofSubs; 
private boolean collision = true; 
private int counter, 


private int collisionCount = 0; 

private Vec3d surfaceNormal = null; 

private Vec3d position = null; 

private Vec3d segment] = null; 

private Vec3d segment2 = null; 

private double distance = 0.0; 

private double targetDistance[] - new doubl|MAX TARGETS]; 
private int skipCount; 

private double minDistance, 
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APPENDIX B. EXAMPLE APPLICATIONS 


A. EXAMPLEBEAMDYNAMIC.JAVA 


/* 
File: ExampleBeamDynamic.java 
Compiler: jdk1.1.6 
21 

package mil.navy.nps.rra; 

import mil.navy.nps.rra. *, 

import java.util.Date; 

Je E 


*(g version 1.0 


*(Qauthor LT Timothy M. Holliday (<A HREF="http://www.stl.nps.navy.mil/~auv/holliday"> 
http://www. stl.nps.navy.mil/—auv/holliday</A>) 
* 


*<dt><b>Location:</b> 
*<dd<a href="http://www.stl.nps.navy.mil/dis-java-vrml/miV/navy/nps/rra/ExampleBeamDynamic.java"> 
* http://www.stl.nps.navy.miVdis-java-vrml/mil/navy/nps/rra/ExampleBeamDynamic.java</a> 


* 
*<dt><b>Hierarchy Diagram:</b> 


*«dd»«a href-"images/RRAClassHierarchy.gif"-IMG SRC-"images/RRAClassHierarchyButton.jpg " 
ALIGN=ABSCENTER></a> 
d 


*<dt><b>Coordinate System Diagram:</b> 


*<dd><a href="images/CoordinateSystem. gif"><IMG SRC="images/CoordinateSystemButton. jpg" 
ALIGN=ABSCENTER></a> 
* 


*<dt><b>ExampleBeamStatic:</b> 


*<dd><a href="ExampleBeamDynamic. wrl"><IMG SRC="images/ExampleBeamDynamic. jpg" 
ALIGN=ABSCENTER></a> 
* 


*<dt><b>Summary:</b> 
*<dd>Produces an example VRML scene showing a beam. 
* 


*<dt><b>Explanation:</b> 
*<dd>Shows a VRML scene with a surface, bottom and beam. The beam 
* show the volume traced out by the pulse during its transmission.<P> 


s 


*<dt><b>History:</b> 

*<dd> 15Apr97 /Timothy M. Holliday /New 
* 

*@see Beam 

2 


public class ExampleBeamDynamic ( 


px 
* This method initializes the environment and the sonar Lobe and 
* has the VRML reprsentation of each written to the console. 
er 
public static void main(String[] args) { 


SSP ssp - new SSP("traditional"), 
Bottom bott = new Bottom("slope",2000); 
Surface surf = new Surface("smooth"); 
Beam beam = new Beam(); 


151 


beam.setElevation( 80); 
beam.setAzimuth(20); 

beam. setHalfBeam Width Y(2); 
beam.setHalfBeamWidthX(2), 
beam.setPosition(-1000, 50, 0); 
beam.setDuration(1); 
beam.setBottom(bott); 
beam.setSurface( surf); 
beam.setSsp(ssp); 
beam.setEndTime(8); 
beam.reset(); 


// Create object with no targets for methods requiring a target 
// as a parameter 
Targets targets = new Targets(surf,bott,0,0); 


// record start time 
Date timecheck = new Date(); 


// print out a VRML header 
System.out.printIn("* VRML V2.0 utf8"); 
System.out. printin("\niniine {url{\"Header. wrl\"]}\n"); 


// call each class’ VRML print routine 
System.out.printin(bott. VRMLBottom()); 
System. out.println(surf. VRMLSurface()); 


System.out.printIn(beam.pıngTımerVRML()), 
beam calculateBeam(targets); 
beam.dynamic VRML(); 


// record stop time 
Date timecheck2 = new Date); 


System.out.printin( "# " + timecheck ); 
System.out.printin( "# " + tumecheck2 ); 


B. EXAMPLEBEAMSTATIC.JAVA 


/ شود‎ 
Pile: ExampleBeamStatic.java 
Compiler: jdk1.1.6 
E 


package mil.navy.nps.rra; 


import mil.navv.nps.rra.*; 
import java.util. Date; 


jae 

*@version 1.0 

*@author LT Timothy M. Holliday (<A HREF="http://www.stl. nps.navy.mil/~auv/holliday"> 
http://www.stl.nps.navy.mil/~auv/holliday</A>) 

* 


*<dt><b>Location:</b> 


*<dd><a href="http://www.stl.nps.navy.miVdis-java-vrml/mil/navy/nps/rra/ExampleBeamStatıc.Java"> 
* http://www.sti.nps.navy.mil/dis-java-vrml/miVnavy/nps/rra/ExampleBeamStatic.Java</a> 
* 
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*<dt><b>Hierarchy Diagram:</b> 

*<dd><a href="images/RR AClassHierarchy.gif"><IMG SRC="images/RRACblassHierarchyButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/CoordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>ExampleBeamStatic:</b> 

*<dd><a href="ExampleBeamStatic.wrl"><IMG SRC="images/ExampleBeamStatic. jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Summary:</b> 

*<dd>Produces an example VRML scene showing a beam. 

* 

*<dt><b>Explanation:</b> 

*<dd>Shows a VRML scene with a surface, bottom and beam. The beam 
* show the volume traced out by the pulse during its transmission.<P> 


* 


*<dt><b>History:</b> 

*<dd> 15Apr97 /Timothy M. Holliday /New 
* 

*@see Beam 

*/ 


public class ExampleBeamStatic { 


per 
* This method initializes the environment and the sonar Beam and 
* has the VRML representation of each written to the console. 
© 
public static void main(String[] args) { 


Vec3d pos = new Vec3d( -1000, 50, 0), 
SSP ssp = new SSP("traditional"); 
Bottom bott = new Bottom("slope",2000); 
Surface surf = new Surface("smooth"); 


Beam beam - new Bean), 
beam.setElevation(80); 
beam.setAzimuth(20), 
beam.setHalfBeamWidthY(2), 
beam.setHalfBeam WidthX(2), 
beam.setPosition(-1000, 50, 0), 
beam.setDuration( 1 ), 
beam.setBottom(bott); 
beam.setSurface( surf), 
beam.setSsp(ssp); 
beam.setEndTime(8); 
beam.reset(); 


// record start time 
Date timecheck 7 new Date(); 


PrintVRML.setColorScheme(Print VRML.RAINBOW 40,100); 
Print VRML. setIntensityScheme(Print VRML.CONSTANT ); 


// print out à VRML header 
System.out.printin(Print VRML .header()); 
System.out.printin(PrintVRML.navigationinfo()); 

System out.printin(Print VRML. directionalLight()); 
System.out.printin(Print VRML .protoColorBar("TL[db]","","","")), 
System.out.printin(Print VRML .colorBarView("From the East", 
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-500, -500, 14500, 

0, 0, 1, 0); 
System.out.printin(Print VRML.colorBar View("From the South", 

-14500, -500, 0, 

0, 1,0, -1.57)), 
System.out.printin(PrintVRML.colorBarView("From the West", 

-500, -500, -14500, 

0, 1, 0, 3.14); 
System.out.printin(Print VRML.colorBar View("From the Air", 

-500, 14500, 0, 

1707 0991-52): 


// call each class'' VRML print routine 
System.out.printin(bott. VRMLBottom()); 
System.out.println(surf. VRMLSurface( )); 
PrintVRML.setColorScheme(PrintVRML.RAINBOW 8,0); 
PrintVRML.setIntensityScheme(PrintVRML.CONSTANT); 
System.out.printin(Print VRML.protoColorBar("TL [dB}","","","")); 
System.out.println(beam. static VRML(Beam. T_L,Beam.NONE)); 
System.out.printin(Print VRML. printLegendView("ExampleBeamStatic", 

" This is an example beam.", 

"It has an initial elevation”, 

"of 80 degrees from", 

"vertical and 20 degrees", 

"west of north.", 

-500, -500, 15000) ); 


// record stop time 
Date timecheck2 = new Date( ); 


System.out.println( "# " + tumecheck ); 
System.out.println( "# "+ timecheck2 ); 


C EXAMPLELOBEDYNAMIC.JAVA 


7 * 
File: ExampleLobeDynamic.java 
Compiler: jdk1.1.6 
E" 


package mil.navy.nps.rra; 
import mil.navy.nps.rra.*; 


import java.util.Date; 


[=> 

*@version 1.0 

*@author LT Timothy M. Holliday (<A HREF="http://www.stl.nps.navy.mil/~auv/holliday"> 
http://www. sti.nps.navy.mil/-auv/holliday</A>) 

* 


*<dt><b>Location:</b> 

*<dd><a href="http://www.sti.nps.navy.mil/dis-Java-vrml/mil/navy/nps/rra/ExampleLobeDynamic.java"> 
* http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/ExampleLobeDynamic.java</a> 

* 


*<dt><b>Hierarchy Diagram:</b> 
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*<dd><a href="images/RRAClassHierarchy.gif'"><IMG SRC="images/RRAClassHierarchyButton.jpg" 
ALIGN=ABSCENTER></a> 
* 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/CoordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

= 


*<dt><b>ExampleLobeDynamic:</b> 

*<dd><a href="ExampleLobeDynamic. wrl"><IMG SRC="images/ExampleLobeDynamic.jpg" 
ALIGN=ABSCENTER></a> 

x 

*<dt><b>Summary:</b> 

*<dd>Produces an example VRML scene showing a lobe. 

* 

*<dt><b>Explanation:</b> 

*<dd>Shows a VRML scene with a surface, bottom and lobe. The scene 

* show the pulse during its transmission as it travels. 


* A lobe is composed of several individual beams.<P> 
* 


*<dt><b>History:</b> 

*<dd> 15Apr97 /Timothy M. Holliday /New 
* 

*@see Lobe 

E 


public class ExampleLobeDynamic { 


pen 
* This method initializes the environment and the sonar Lobe and 
* has the VRML reprsentation of each written to the console. 
”/ 
public static void main(String[] args) { 


Vec3d pos = new Vec3d( -1000, 50, 0); 
SSP ssp = new SSP("traditional"); 
Bottom bott = new Bottom("slope",2000); 
Surface surf = new Surface("smooth"); 


Lobe lobe 7 new Lobe(); 
lobe.setElevation(80); 
lobe.setAzimuth(20), 
lobe.setLobe Width Y (4); 

lobe. setLobe Widthx(4),; 
lobe.setNumber YPartition(2), 
lobe. setNumberXPartition(2); 
lobe. setPosition(pos); 

lobe. setDuration(1); 
lobe.setBottom(bott); 
lobe.setSurface(surf), 

lobe. setSsp(ssp); 
lobe.setEndTime($8), 
lobe.reset( ); 


// Create object with no targets for methods requiring a target 
// as a parameter 
Targets targets = new Targets(surf,bott,0,0); 


// record start time 
Date timecheck = new Date(); 


// print out a VRML header 

System.out. printin(PrintVRML.header()); 
System.out. printin(PrintVRML.navigationinfo()); 
System.out. printin(Print VRML. directionalLight()); 
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System.out.printin("\nInline {urif\"Header.wr!\"]}\n"); 


// call each class’ VRML print routine 

System.out.println(bott. VRMLBottom()); 

System.out.printin( surf. VRMLSurface()); 

lobe.calculateLobe(targets); 

System.out.printin(lobe.dynamicVRML()); 

System.out.println(PrintVRML.printLegendView("ExampleLobeDynamic", 
" This 1s an example lobe.", 
"It is made of 4 beams", 
"Elevation - 80 degrees", 
"from vertical and 20", 
"degrees west of north.", 
-500, -500, 15000) ); 


// record stop time 
Date timecheck2 = new Date(); 


System.out.println( "2 " + tumecheck ); 
System.out.println( "# " + timecheck2 ); 


D. EXAMPLELOBESTATIC.JAVA 


/۴ 
File: ExampleLobeStatic. java 
Compiler: jdk1.1.6 
4 


package mil.navy.nps.rra; 
import mil.navy.nps.rra.*; 
import java.util.Date; 


7 
*@version 1.0 


*@author LT Timothy M. Holliday (<A HREF="www.sti.nps.navy.mil/-auv/holliday"> 
http://www.sti.nps.navy.mil/-auv/holliday</A>) 
* 


*<dt><b>Location:</b> 

*<dd<a href-"http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/ExampleLobeStatic.java" 
* http//www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/ExampleLobeStatic.java«/a7 

* 


*<dt><b>Hierarchy Diagram:</b> 

*<dd><a href="1mages/RRAC lassHierarchy.gif"><IMG SRC="images/RRAClassHierarchyButton. jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="1mages/CoordinateSystem. gif"><IMG SRC="images/CoordinateSystemButton. jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>ExampleLobeStatic:</b> 
*<dd><a href="ExampleLobeStatic.wrl"><IMG SRC="1mages/ExampleLobeStatic. jpg" ALIGN=ABSCENTER></a> 


* 


*<dt><b>Summary:</b> 
*<dd>Produces an example VRML scene showing a lobe. 


x 
*<dt><b>Explanation:</b> 
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*<dd>Shows a VRML scene with a surface, bottom and lobe. The lobe 
* show the volume traced out by the pulse during its transmission. 
* Alobe is composed of several individual beams.«P» 


* 


*<dt><b>History:</b> 

*<dd> 15Apr97 /Timothy M. Holliday /New 
* 

*@see Lobe 

og | 


public class ExampleLobeStatic { 


Ez 
* This method initializes the environment and the sonar Beam and 
* has the VRML representation of each written to the console. 
4 
public static void main(String[] args) { 


Vec3d pos 7 new Vec3d( -1000, 50, 0); 
SSP ssp - new SSP("traditional"); 
Bottom bott = new Bottom("slope",2000); 
Surface surf = new Surface("smooth"), 
Targets targets - new Targets(surf,bott,0,0); 


Lobe lobe = new Lobe(); 
lobe.setElevation(80), 
lobe.setAzimuth(20), 
lobe.setLobeWidth Y(4); 
lobe.setLobeWidthX(4); 
lobe.setNumber YPartition(5); 
lobe.setNumberXPartition(5); 
lobe.setPosition(pos ), 
lobe.setDuration( 1 ); 
lobe.setBottom(bott), 

lobe. setSurface(surf); 
lobe.setSsp(ssp); 
lobe.setEndTime(3 ); 
lobe.reset( ); 
lobe.calculateLobe( targets); 


// record start time 
Date timecheck = new Date( ), 


// print out a VRML header 
System.out.printIn( Print VRML. header()); 
System.out.printin( Print VRML.navigationinfo( )); 
System.out.printin(Print VRML directionalLight( )), 
System. out. println(bott. VRMLBottom()); 
System.out.printin(surf. VRMLSurface()); 


/* System.out.printin("#Scheme 1"); 
Print VRML.setColorScheme(PrintVRML.RAINBOW 8,0), 
PrintVRML.setIntensityScheme(PrintVRML.CONSTANT), 
System.out.printIn(Print VRML.protoColorBar("Time [sec}","","","")), 
System.out.printin(lobe.staticVRML{Beam. TIME ‚Beam. NONE)); 


System.out.printin("#Scheme 2"), 
PrintVRML.setColorScheme(PrintVRML.RAINBOW 40,100); 
PrintVRML.setIntensityScheme(PrintVRML.CONSTANT), 
System.out.printin(PrintVRML.protoColorBar("TL[dB]","","","")); 
System.out.println(lobe.static VRML(Beam.T L,Beam.NONE)); 

2 
System.out.println("ZScheme 3"); 
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/ 


= 


PrıntVRML.setColorScheme(Print VRML.RAINBOW 40,100), 

Print VRML.setIntensityScheme(Print VRML.LINEAR,0,8); 
System.out.printin(PrintVRML.protoColorBar("TL[dB]","Time [sec]","","")); 
System.out.printin(lobe.static VRML(Beam.T L,Beam. TIME)), 


System.out.println("? Scheme 4"); 

Print VRML.setColorScheme(Print VRML.RAINBOW, 8,0); 

Print VRML. setintensityScheme(Pnnt VRML.LINEAR,40, 100); 
System.out.printIn(PrintVRML.protoColorBar("Time [sec]","TL[dB]","","")); 
System.out.printIn(lobe.static VRML(Beam. TIME,Beam.T_L)); 


System.out.printin("#Scheme 5"); 

Print VRML.setColorScheme(PrintVRML. TRICOLOR,40,60,80, 100); 

Print VRML.setIntensityScheme(PrintVRML.LINEAR,0,8), 

System.out.printin(Print VRML.protoColorBar("","Time [sec]","Detect","Counter Detect")); 
System.out.println(lobe.staic VRML(Beam.T L,Beam.TIME)), 


System.out.printin("#Scheme 6"), 

Print VRML.setColorScheme(Print VRML. TRICOLOR,40,60,80,100); 

Print VRML.setIntensityScheme(Print VRML.LINEAR,40,100); 
System.out.printin(Print VRML. protoColorBar("","TL[dB]","Detect","Counter Detect")); 
System.out.println(lobe. static VRML(Beam.T L,Beam.T_L)); 


System.out.printIn("#Scheme 7"), 
PrintVRML.setColorScheme(PrintVRML.RED,40,100), 

Print VRML. setIntensityScheme(PrintVRML.LINEAR,40,100); 
System.out.printin(Print VRML. protoColorBar("","TL[db]"."","")y 
System.out.printin(lobe.staticVRML(Beam.NONE,Beam.T_L)); 


System.out.printin(Print VRML. colorBarView("From the East", 
-500, -500, 14500, 
0, 0, 1, 0)); 
System.out.printIn(PrintVRML.colorBarView("From the South", 
-14500, -500, 0, 
0, 1,0, -1.57)); 
System. out.printin(Print VRML. colorBarView("From the West", 
-500, -500, -14500, 
0, 1, 0, 3.14), 
System.out.printin(Print VRML. colorBarView("From the Aur", 
-500, 14500, 0, 
1, 0, 0, -1.57((; 
System. out. printin(Print VRML. printLegend View("ExampleLobeStatic", 
" This is an example lobe.", 
"It is made of 4 beams", 
"Elevation - 80 degrees", 
"from vertical and 20", 
"degrees west of north.", 
-500, -500, 15000) ); 


// record stop time 
Date timecheck2 = new Date( ); 


System.out.printIn( "# " + timecheck ); 
System.out.printIn( "# " + timecheck2 ); 
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E. EXAMPLERAY.JAVA 


Je 
File: ExampleRay.java 
Compiler: jdk1.1.6 
d 


package mil.navy.nps.rra; 
import mil.navy.nps.rra.*; 
import java.util.Date; 


/* * 

*(Q)version 1.0 

*(Qauthor LT Timothy M. Holliday (<A HREF="http://www. stl.nps.navy.mil/—auv/holliday"> 
http://www. stl.nps.navy.mil/—auv/holliday</A>) 

* 

*<dt><b>Location:</b> 

*<dd><a href="http://www.stl.nps.navy.mildis-java-vrml/mil/navy/nps/sra/ExampleRay.java"> 

*. http//www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/ExampleRay.java«/a» 

* 

*<dt><b>Hierarchy Diagram:</b> 

*<dd><a href="images/RRAClassHierarchy.gif'><IMG SRC="images/RRAClassHierarchyButton.jpg" 
ALIGN=ABSCENTER></a> 

x 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/CoordinateSystemButton jpg” 
ALIGN=ABSCENTER></a> 

x 


*<dt><b>ExampleRay:</b> 

*<dd><a href="ExampleRay.wrl"><IMG SRC="images/ExampleRay.jpg" ALIGN=ABSCENTER></a> 
x 

*<dt><b>Summary:</b> 

*<dd>Produces an example VRML scene showing a ray. 

* 

*<dt><b>Explanation:</b> 

*<dd>Shows a VRML scene with a surface, bottom and ray.<P> 
* 

*<dt><b>History:</b> 

*<dd> 15Apr97 /Timothy M. Holliday /New 


* 


*@see Ray 
i 
public class ExampleRay (1 


/* * 
* This method initializes the environment and the sonar Ray and 
* has the VRML reprsentation of each written to the console. 
20 
public static void main(String[] args) { 


SSP ssp - new SSP("traditional"), 
Bottom bott = new Bottom("slope",2000); 
Surface surf - new Surface(" smooth", 


Ray ray = new Ray(); 

ray.setPosition( -1000, 50, 0); 
ray.setElevation(80); 
ray.setAzimuth(20); 
ray.setDeltaTime(.006); : 
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ray.setDuration( 1); 
ray.setBottom(bott); 
ray.setSurface(surf), 
ray.setSsp(ssp); 
ray.reset(); 


mt j = 0; 


// record start time 
Date timecheck 7 new Date(); 


// print out a VRML header 
System.out.printin(Print VRML. header()); 

System. out. printin(Print VRML. navigationinfo()); 
System.out.printin(Print VRML. directionalLight( )); 
System.out.printin("\nInline {url[\"Header.wri\"]}\n"); 


// call each class’ VRML print routine 
System.out.printin(bott. VRMLBottom()); 
System.out.printin(surf. VRMLSurface()); 


// propagate the ray and print it 
System.out.printin("Transform (" ); 
System.out.printIn(" rotation 1 0 0 3.14" ), 
System.out.printIn(" children" ); 


System.out.printin("Shape (" ); 
System.out.println(" appearance Appearance ("); 
System.out.printin(" material Material ("); 
System.out.printin" ^ emissiveColor 1.0 0.0 0.0"); 
System.out.printin(" — diffuseColor 0.0 0.0 0.0"); 
System.out.println(" }"); 

System.out.printin(" 3"); 

System.out.printin(" geometry IndexedLineSet {"), 
System.out.printin(" coord Coordinate {"); 
System.outprintln(" point ["); 

I 


// Record first point always 
ray.recordPoint(); 
dues 


while (ray.getTime() « 8.0) 1 


/| Propagate all rays through one time step 
ray.Propagate(.006); 


// If any ray has reflected then record all of the 
// points in the beam 
if ( ray.reflected()) 1 

ray.recordPoint(); 

pet 


j 


// Jf the curvature sum reaches the limit, record 
// all of the points in the beam 
if ( ray.totalCurvature() > .022 ) { 
ray.recordPoint(); 
m 
} 
j 
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// Record last point always 
ray.recordPoint(); 


// Print the positions stored in each ray 

for (j=0; j<ray.getCount() J] += 1) { 
System.out.print(ray.position(])), 

j 


System.out.printin("  ]"); 
System.out.printin(" 3") 
System.out.printIn("  coordIndex ["); 


// Define the segments of the ray 

for(j=0; j<ray.getCount(); j++) { 
System.out.println(j); 

j 

System.out.printin(" ]"); 

System.out.printin(" 3"); 

System.out.printin("}"); 

System.out.printin("}"),; 


System.out.printiIn(PrintVRML.printLegendView("ExampleRay", 
" This is an example ray.", 
"It has an initial elevation", 
"of 80 degrees from", 
"vertical and 20 degrees”, 
"west of north.", 
-500, -500, 15000) ); 


// record stop time 
Date timecheck2 = new Date( ); 


System.out.println( "£ " * timecheck ), 
System.out.printIn( "# " + timecheck2 ); 


F. PINGER.JAVA 


File: Pinger.java 
Compiler: jdk1.1.6 


package mil.navy.nps.rra; 


import java.util Date; 


import mil.navy.nps.rra.*; 


Ja 


*@version 1.0 


*@author LT Timothy M. Holliday (<A HREF="http://www.stl.nps.navy.mil—auv/holliday"> 


http://www.stl_nps.navy.mil/~auv/holliday</A>) 
* 


*<dt><b>Location:</b> 


*«dd»«a hrefz"http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/Pinger.java"7 
* http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/Pinger.java</a> 
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* 


*<dt><b>Hierarchy Diagram:</b> 

*<dd><a href="1mages/RRAClassHierarchy. gif"><IMG SRC="images/RRAClassHierarchyButton jpg" 
ALIGN=ABSCENTER></a> 

x 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/CoordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Summary:</b> 

*<dd>This stand alone program initializes the environment (bottom, surface and 
* SSP) and the sonar array radiation lobe. 

x 


*<dt><b>Explanation:</b> 

*<dd>This program provides a template for running a stand-alone simulation 

* that can run without user interaction to generate sonar pings and returns for a 
* virtual world.<P> 


* 


*<dt><b>History:</b> 

*<dd> 15Jan98 /Timothy M. Holliday /New 

*<dd> 17Mar98 /Timothy M. Holliday /Added HTML comment convention 
* 

*@see Lobe 

*/ 


public class Pinger { 


public statıc final double MINIMUMX = -5000; 
public static final double MINIMUMZ = 4000; 
public statıc final double MAXIMUMZ = 5000; 
public statıc final double MAXIMUMX = -4000; 
public static final double NORTH = 000.0; 
public static final double SOUTH = 180.0; 
public static final double EAST = 090.0; 

public static double deltaTime = .8; // in seconds 
public static double lastZLane = 4950; 

public static Vec3d position = null; 

public static Vec3d deltaPosition = null; 

public static Vec3d velocity = null; 

public static double orientation = 0.0; 


per 


* This method initializes the environment and the sonar Lobe and 
* has the VRML reprsentation of each written to the console. 
n 

public static void main(String[] args) { 


// Position velocity and orientation in Ziomek Coordinate System 
position = new Vec3d( -5000, 50, lastZLane); //in meters 
deltaPosition = new Vec3d( 0, 0, 0); //in meters 

velocity = new Vec3d( 2.8, 0, 0), // in m/s which is 5 kts 
orientation = NORTH; // 1n degrees 


SSP ssp = new SSP("constant"); 

Bottom bott = new Bottom("noslope",100); 

Surface surf = new Surface("smooth”), 

Lobe lobe = new Lobe); 

Targets targets = new Targets(surf,bott, MAXIMUMX,MINIMUMX,MAXIMUMZ,MINIMUMZ, 50), 


lobe.setElevation(90 y, 
lobe.setAzimuth(0); 
lobe.setLobe Width Y(40), 
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lobe.setLobeWidthX(52); 
lobe.setNumberYPartıtion(4); 
lobe.setNumberXPartition(4); 
lobe.setDuration( 1); 
lobe.setBottom(bott); 

lobe.setSurface( surf); 

lobe.setSsp(ssp); 

lobe.setEndTime(.133); // 200m/1500m/s 


// record start time 
Date timecheck 7 new Date(); 


while (notDoneSearching()) { 
lobe.setPosition(position); 
lobe.setAzimuth(orientation); 
lobe.reset(); 
//System.out.println("calculating"); 
lobe.calculateLobe(targets); 
update VehiclePosition(); 


j 


targets.printMineStatistics(); 
// record stop time 
Date timecheck2 = new Date(); 


System. out.printin( "4 " + timecheck ), 
System.out.println( "2 " * timecheck2 ); 


j 


public static boolean notDoneSearching() { 
if (position.get(2) <= MINIMUMZ) { 
return false; 


eise { 
return true; 
j 
j 


public static void update VehiclePosition () { 
if (orientation == NORTH) ( 
if (position.get(0) >= MAXIMUMX) { 
orientation = EAST; 
velocity.set(0,0,-2.8); 
position. print(); 
j 


j 
else if (orientation == SOUTH) { 
if (position. get(0) <= MINIMUMX) { 
orientation = EAST; 
velocity.set(0,0,-2.8); 
position. print(); 
} 


else if (orientation == EAST) { 
if (position. get(2) <= (lastZLane-100.0)) { 

if (position. get(0) >= MAXIMUMX) 1 
orientation = SOUTH; 
position. print(); 
velocity.set(-2.8,0,0); 

j 

else { 
orientation = NORTH; 
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velocity.set(2.8,0,0); 
position.print(); 


j 
lastZLane = position.get(2); 
j 


j 
deltaPosition.set( velocity); 


deltaPosition.scale(delta Time); 
position.add(deltaPosition); 
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APPENDIX C. CLIENT SERVER CODE 


A. BATTLESCENE.JAVA 


ye 

File: BattleScene.java 

Compiler: jdk1.1.6, Netscape 3.0, World View 2.0 
er 


import java.net.*; 
import java.10. *; 
import java.util. *; 
import vrml. *; 
import vrml.node. *; 
import vrml. field. *; 


import mil.navy.nps.rra. *; 


Jer 
*@version 1.0 


*@author LT Timothy M. Holliday (<A HREF="http:/Awww.stl.nps.navy.mil/~auv/holliday"> 
http://www.stl.nps.navy.mil/-auv/holliday</A>) 
* 


*<dt><b>Location:</b> 

*<dd><a href="http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/BattleScene.java"> 
* http//www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/BattleScene.java«/a» 

* 


*<dt><b>Summary:</b> 

*<dd>Interfaces between the browser and pinserver. This script also provides for 

* platform mobility and controls all timing of the interactions between ping 

* requests, VRML returns and platform movement. 

* 

*<dt><b>Explanation:</b> 

*<dd>A VRML world is in essence a static world, with some basic movements. To be able 
to emulate more real world processes the VRML world needs to have a script 
connection to a program running outside of the VRML rendering engine, that is 
running in the browsers virtual machine. This program is that extension of a 
particular VRML world. It also provides network connections to a sonar ping 

server runing on a separate computer. This network connection is provided since 
both rendering VRML and calculating real-time sonar are computationally intensive. 
Thus, the two functions are distributed two two machines over the net.<P> 


*X t * XX * A X* X 


*<dt><b>History:</b> 
*<dd> 15Jan98 /Timothy M. Holliday /New 
*<dd> 17Mar98 /Timothy M. Holliday /Added HTML comment convention 
x 
* @see BridgeServer 
*@see PingServer 
* / 
public class BattleScene extends Script implements Timed( 


// Browser Specific Stuff 

private Browser browser, 

ControlPanel panel; 

MFNode addChildren, removeChildren; 

MFNode addBeamChildren, removeBeamChildren; 
SF Vec3f beamTranslation; 

SFFloat yaw; 

SF Vec3f translation; 
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SFFloat beamYaw; 

BaseNode nodes|[] = null; 

private BaseNode oldBeamNodes[]{} = {null,null,null null null}; 
private BaseNode oldTargetsNodes[][] = {null,null,null,null,null}; 
private int activeBeamNode = 0; 

private int activeTargetsNode = 0; 


//Network specific stuff 

public final static int PORT = 4129; 

public final static String HOST = "localhost"; // server host name 
public final static float SET SONAR_PARAMETERS = (float)0.0; 
public final static float DETECT_TARGETS = (float)1.0; 

public final static float DETECT_BEAM = (float)2.0; 

public final static float DETECT BOTH - (float)3.0; 

Socket socket 7 null; 

DatalnputStream ın = null; // input stream from server to client 
DataOutputStream out = null; // output sream from client to server 


// Gereral Stuff specific to the vehicle in question 

private final static boolean DEBUG = true; 

private final static boolean LOCAL - false; 

private final static float START. DEPTH - 500; 

private final static float START X POS = -5000; 

private final static float START Y POS - 0; 

private float orderedSpeed = 0; 

private float orderedDepth = 0; 

private float orderedCourse = 0; 

private float position{] = new float[3]; // naval coordinate system 
private float drPosition[] = new float[3]; // naval coordinate system 
private float heading = 0; 

private float drHeading 7 0; 

private float speed = 0; 

private float drSpeed = 0; 

private float depth - START DEPTH; 

private float drDepth = 0; 


//Timing specific stuff 

private final static int interval- 100; 

private int tickCount = 0; 

private int pingOrderTick = 0; 

public boolean createVRMLThreadRunning = false; 
private float pingInterval = (float)4.0; 

private CreateVRMLBeamThread beamThread = null; 
private Create VRML TargetsThread targetsThread = null; 


pe 

* This method initializes the Java script. It establishes connections 
* to the browser, the control panel, the timer and to BridgeServer. 
E 

public void initialize()( 


if (DEBUG) 
System.out.println("BattleScene: initialize: starting"); 


// Connects to the browser and various VRML nodes in the browser 
browser = getBrowser(); 

// get the reference to the target node. 

Node root = (Node)((SFNode)getField("root”)).get Value(); 

Node manta = (Node X(SFNode)getField("manta")).getValue(); 
Node beam = (Node X(SFNode)getField("beam”)).getValue(); 


// Initiates the control panel 
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panel = new ControlPanel(thıs), // start the panel. 


// get the references to exposed fields, event-in/outs in the VRML nodes. 
addChildren = (MFNode)root.getEventIn("addChildren"); 
removeChildren = (MFNode)root.getEventIn("removeChildren"); 
addBeamChildren = (MFNode)beam.getEventIn("addChildren"); 
removeBeamChildren = (MFNode)beam.getEventIn("removeChildren"); 
beamTranslation = (SF Vec3 f)beam.getExposedField("translation"); 
beam Yaw - (SFFloat)beam.getEventIn( " yaw"); 

translation 7 (SFVec3f)manta.getExposedField( "translation" ); 
position[0]= translation.getX(); 

position[1]= translation.getZ(); // translate from VRML to NAVAL system 
position[2]= -translation.getY(); 

yaw = (SFFloat)manta. getEventin("yaw”); 


if (LOCAL) 1 
// open socket connections and establish data streams with BridgeServer 
try { 
// open network and input/output stream 
socket = new Socket("localhost", PORT); 
in = new DatalnputStream(socket.getInputStream()); 
out = new DataOutputStream(socket. getOutputStream()); 
} catch (UnknownHostException e) { 
browser.setDescription("Unknown host: " + HOST); 
} catch (Exception e) ( 
browser.setDescription( " Connection error"); 


j 


// initiate a callback timer 
new TimerThread(this, interval ).start(); 


depth 7 position[2 |; 
orderedDepth = depth; 


if (DEBUG) 
Systern.out. printin("BattleScene: initialize: ending"); 
j 


en 

* This method receives ticks from TimerThread and updates vehicle 
* position in the VRML scene. 

*/ 

public void tick( TimerThread t) | 


if (orderedSpeed-speed > .05) { 
speed += 0.1; 


} ۰ 

else if (orderedSpeed-speed < -.05) { 
speed -= 0.1; 

j 


if (orderedCourse-heading > .05) { 
heading += 0.1; 

j 

else if (orderedCourse-heading < -.05) 1 
heading -- 0.1; 

j 


position[0] += Math.cos(heading*3.14/180.0)*.1*speed; 
position[1] += Math.sin(heading*3.14/180.0)*. 1*speed; 
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j 


if (orderedDepth-depth > .05) { 
depth += 0.1; 

} 

else if (orderedDepth-depth < -.05) { 
depth = 0.1; 

j 


position[2] = depth; 


if (tickCountt+-%10 == 0) { 
panel.updateShipState(heading, speed, depth); 
browser.setDescription(" X=" + (int)position[0] + 
" Y=" + (int)position[1] + 
" Z=" + (int)positionf2]); 
update VehiclePosition(); 
j 


سی 


* 
* 
* 
* 
* 
* 
* 


This method calculates the position in the future that the vehicle 
will be at when the pulse 1s sent out. 


Since it takes a finite time to calculate sonar pings, all requests 
made to the server are future requests. This method calculates the 
vehicles position at that future position so that 1t matches the 
sonar pings starting position 


e, 
private void deadReckon(float interval) { 


float timer = 0; 


drPosition[0] = position[0]; 
drPosition[1] = position[1]; 
drPositionf2] = positionf2]; 
drSpeed = speed; 
drHeading = heading; 
drDepth = depth, 


while (timer « interval) ( 


if (orderedSpeed-drSpeed > .05) { 
drSpeed += 0.1; 


j 

else if (orderedSpeed-drSpeed « -.05) ( 
drSpeed = 0.1; 

j 


if (orderedCourse-drHeading > .05) { 
drHeading += 0.1; 


j 

else if (orderedCourse-drHeading < -.05) { 
drHeading = 0.1; 

j 


drPosition[0] += Math.cos(drHeading*3.14/180.0)*.1*drSpeed; 
drPosition[ 1] += Math.sin(drHeading*3.14/180.0)*.1*drSpeed; 


if (orderedDepth-drDepth > .05) { 
drDepth += 0.1; 


j 

else 1f (orderedDepth-drDepth « -.05) ( 
drDepth = 0.1; 

} 


drPosition[2] = drDepth; 
timer += 0.1; 
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j 


j 


je 


* This method processes events generated by the VRML Browser. 


* When the VRML world is entered, an 'entered' event is sent out 


* and this event handler processes ıt and starts the whole 


* simulation 1n motion. 
m 


public void processEvent(E vent ev) ( 


if (DEBUG) 
System.out.printin("BattleScene: processEvent: starting"); 


if(ev.getName().equals("entered"))( 
ConstSFBool v 7 (ConstSFBool)ev.get Value(); 
String receivedString = "#VRML V2.0 utf8\n"; 


// get VRML data from whichever source it is necessary 
if (LOCAL) 1 


try { 
boolean start = true; 


// send start Order 
out. writeBoolean( start); 


// receive new VRML from server 
receivedString += in.readLine(); 


browser.setDescription("Got a String!"); 
y catch (IOException e) { 
browser.setDescription("IOException: " 4 e); 


j 


// create VRML and add it to the appropriate VRML node 
try ( 
nodes = browser.create VrmlFromString(receivedString); 
if (null'- nodes) { 
addChildren.set Value(nodes), 
j 


j 

catch (Exception e) { 
browser.setDescription("can not create VRML node"), 
e.printStackTrace(); 


j 
} // end DEBUG 


// initialize the control panel and draw it 
if(v.get Value()){ panel.map(); } 
else { panel. hide(); } 


j 


if (DEBUG) 
System.out.println("BattleScene: processEvent: ending"); 


j 


yrs 


* This method cleans up memory and closes sockets at shutdown 
ui 


public void shutdown()( 


panel.dispose(); 
try { 
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out.close(); 
in.close(); 
socket.close( ); 
) catch (Exception e) ( 
browser.setDescription("Connection close error"); 
j 
j 


72 
* This method is called by the control panel to set the 
* requested maneuver in BattleScene. 
" 
public void setManeuver(float d, 
float s, 
float c) 1 
orderedDepth - d; 
orderedSpeed - s; 
orderedCourse - c; 


j 


per 


* This method is called by the control panel to set the 
* requested sonar ping inn BattleScene. 


2 
public void setSonarParameters(float elevation, 
float azimuth, 
float verticalBeamWidth, 
float horizontalBeamWidth, 
float verticalBeamConfiguration, 
float horizontalBeamConfiguration, 
float verticalSearchWidth, 
float horizontalSearchWidth, 
float pingl, 
float powerLevel)f 
if (DEBUG) 
System.out.printin("BattleScene: setSonarParameters: starting"); 
pingInterval — pingl, 
if (LOCAL) ( 
browser.setDescription( "Sending Ping"); 
try { 


out. writeFloat(SET SONAR. PARAMETERS); 
out. writeFloat(position[0]); 

out. writeFloat(position[2 ]); 

out. writeFloat(-position[1]); 

out. writeFloat(elevation ); 

out. writeFloat(-azimuth y, 

out. writeFloat( verticalBeamW idth y, 

out. wniteFloat(horizontalBeam Width); 

out. writeFloat(verticalBeamConfiguration y; 
out. writeFloat(horizontalBeamConfiguration ); 
out. writeFloat((pinglI/(float)2.0)); 

out. writeFloat(powerLevel); 


catch (Exception e) { 
browser.setDescription("IOException: " + e); 
j 
j 


if (DEBUG) 
System.out.printIn("BattleScene: setSonarParameters: ending"); 
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j 
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* This method ıs called by sonarDetection to decide ıf another 
* sonar request can be made. 
% 
* Another sonar request cannot be made if the threads for 
* adding VRML beams and targets are still alive. 
* When the threads die this indicates that the VRML has been 
* added to the scene graph and another request can be made. 
E 
public void waıtOnThreads( ){ 
if (targetsThread != null && targetsThread.isAlive()){ 
try { 
targetsThread.join(), 
y 
catch (InterruptedException e) { 
System.out.printin("BattleScene: waıtOnThreads: targetsThread "+e), 


j 


} 
_ if (beamThread != null && beamThread.isAlive()){ 
try { 
beam Thread.join(); 


j 
catch (InterruptedException e) ( 
System.out.printin("BattleScene: waitOnThreads: beamThread: "+e), 
j 
j 
j 


7۳۰ 

* This method 1s called from the control panel to initiate a 

* sonar ping. When a ping is recieved and the required threads 

* diea new sonar request is made and the corresponding VRML is 
* received. 

t 

public void sonarDetection(float order) ( 


if (DEBUG) 
System.out.printin("BattleScene: sonarDetection: starting"); 


waitOnThreads(); 
String receivedTargets = ""; 
String receivedBeam - ""; 
if (LOCAL) { 
try { 
// send Ping Order 
pingOrderTick = tickCount, 
System.out.printin("BattleScene: sonarDetection: order = "+order); 
out.writeFloat(order), 


// receive desired vrml strings; 
if (order == DETECT_TARGETS) { 
receivedTargets 7 in.readLine(); 


j 
else if (order == DETECT_BEAM) { 
receivedBeam = in.readLine(); 


- if (order -- DETECT BOTH) ( 
receivedTargets = in.readLine( ), 
receivedBeam 7 in.readLine(); 
j 
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j 
catch (IOException e) { 
browser.setDescription("IOException: " + e), 
j 
j 


else { 
// fill with fake VRML so that when the sonar server 1s being run 
// the local machine, the simulation won't bog down. Used for debugging. 
pingOrderTick = tickCount; 
receivedTargets = targets VRML({); 
receivedBeam = beamVRML({); 


targetsThread = new CreateVRMLTargetsThread(this,browser); 
beamThread = new CreateVRMLBeam Thread(this, browser,targets Thread); 


// perform a small wait if the sonar server is operating faster than 
// real time 
while (tickCount-pingOrderTick « (int)pingInterval* 10) ( 

int 1; 

for (1=0;1<10;1++)(3 


j 


// start the threads to produce the VRML, 

if (order -- DETECT TARGETS) | 
targetsThread.start(receivedTargets,old TargetsNodes|[(activeTargetsNode- 1 )%5]); 
updateBeamPosition( ); 


} 
else if (order == DETECT_BEAM) { 
beamThread. start(receivedBeam,oldBeamNodes{(activeBeamNode+1)%1}); 


} 

else if (order == DETECT_BOTH) { 
targetsThread.start(receivedTargets,oldTargetsNodes[(active TargetsNode+ 1 )%5]); 
beam Thread. start(receivedBeam,oldBeamNodes[(activeBeamNode+1)%1]); 
updateBeamPosition(); 

j 


if (DEBUG) 
System.out.printin("BattleScene: sonarDetection: ending"); 


j 


per 


* This method adds the beam VRML to the scene graph 
ah 
public void AddBeamChildren(BaseNode nodes[]){ 


if (DEBUG) 
System.out.println("BattleScene: AddBeamChildren: starting"); 


activeBeamNode++, 
oldBeamNodes[activeBeamNode%1] = nodes; 
addBeamChildren.set Value(nodes), 


if (DEBUG) 
System.out println("BattleScene: AddBeamChildren: ending"); 


j 

[** 

* This method decides whiche beam VRML to remove from the scene graph 
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el 
public void removeBeam() { 
waitOaThreads(); 
while (tickCount-pingOrderTick < (int)pingInterval*20) { 
int 1; 
for (1=0:1<10:147) (3 
j 
RemoveBeamChildren(oldBeamNodesl[activeBeamNode% 1 بر[‎ 


y 


per 

* This method removes the beam VRML from the scene graph 
7 

public void RemoveBeamChildren(BaseNode nodesl]) { 


if (DEBUG) 
System.out.println("BattleScene: RemoveBeamChildren: starting"); 


removeBeamChildren.set Value(nodes), 


if (DEBUG) ۱ 
System.out.printin("BattleScene: RemoveBeamChildren: ending”); 


} 


per 

* This method adds the target VRML to the scene graph 
E 

public void AddTargetsChildren(BaseNode nodes[])1 


if (DEBUG) 
System.out.printin("BattleScene: AddTargetsChildren: starting"); 


active TargetsNode++; 
oldTargetsNodes[activeTargetsNode965] = nodes; 
addChildren.set Value(nodes); 


if (DEBUG) 
System.out.printin("BattleScene: AddTargetsChildren ending"); 


j 


pun 

* This method removes the target VRML from the scene graph 
ai 

public void RemoveTargetsChildren(BaseNode nodes[]) ( 


if (DEBUG) 
System.out.printiIn("BattleScene: RemoveTargetsChildren: starting"); 


removeChildren. set Value(nodes); 


if (DEBUG) i 
System.out.println("BattleScene: RemoveTargetsChildren: ending"); 


j 


private void update VehiclePosition() { 
float transformPosition{] = new float[3}; 
transformPosition[0] = position[0]; 
transformPosition[1] = -position[2]; 
transformPosition[2] 7 position[1]; 
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translation. set Value(transformPosition), 
yaw.setValue(heading*(floatX-3.14/180)); 
return; 

j 


private void updateBeamPosition() { 
float transformPosition[] = new float[3]; 
transformPosition[0] = position[0]; 
transformPosition{1] 7 -position[2]; 
transformPosition[2] = position[ 1]; 


beamTranslation.set Value(transformPosition), 
beam Yaw. set Value(drHeading*(float)(-3.14/180)); 
return: 


private String targets VRML() { 
return “Transform {"+ 
" translation "+ (drPosition[0]+Math.cos(-drHeading)*50) +" "+ 
(-drPosition(2]) +" "+ 
(drPosition{ 1 }+Math.sin(-drHeading)*50) + 
" children ["+ 
" Shape {"+ 
" A appearance Appearance ("+ 
i material Material {"+ 
A emissiveColor 1.0 0 0"+ 
" 1T 
" } "d. 
" A geometry Sphere ( radius 20) "7 
" ) "+ 
" Pr 
D E 
j 


private String beamVRML() 1 
return "Transform {"+ 

" translation 50 50 0"+ 
" children ("+ 
" Shape {"+ 
" appearance Appearance {"+ 
material Material {"+ 
E emissiveColor 0 1.0 0"+ 
1 } "4+ 
" } "+ 
" A geometry Sphere ( radius 20) "7 
" } "4 
" + 


= ES 


B. BRIDGESERVER.JAVA 


Je 
File: BridgeServer.java 
Compiler: jdk1.1.6 
*/ 


package mil.navy.nps.rra; 
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import java.net.*; 
import java.io. *; 
import java. util. *; 


import mil.navy.nps.rra.*; 


¡ee 


*@version 1.0 

*@author LT Timothy M. Holliday (<A HREF="http://www.stl.nps.navy.mil/~auv/holliday"> 
http//www. stl.nps.navy.mil/—auv/holliday</A>) 

x 


*<dt><b>Location:</b> 

*<dd><a href="http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/BridgeServer.java"> 
* http:/Avww.stl._nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/BridgeServer.java</a> 

* 


*<dt><b>Hierarchy Diagram:</b> 

*<dd><a href="images/RRAClassHierarchy.gif"><IMG SRC="images/RRAClassHierarchyButton.jpg" 
ALIGN=ABSCENTER></a> 

Ye 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="1mages/CoordinateS ystem.gif"><IMG SRC="images/CoordinateS ystemButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Summary:</b> 

*<dd>Takes requests for sonar pings and passes them to PingServer and receives VRML 

* from the PingServer and returns it to BattleScene. 

* 

*<dt><b>Explanation:</b> 

*<dd>Since calculation of numerous ray trajectories is extremely time intensive, rendering and 
* ray tracing need to be on separate computers. This provides an interface from the rendering methods 
* to the ray tracing methods. The use of this intermediate bridge avoids the security 

* issues that are present in most browsers.<P> 

* 

*<dt><b>History:</b> 

*<dd> 15Jan98 /Timothy M. Holliday /New 

*<dd> 17Mar98 /Timothy M. Holliday /Added HTML comment convention 
se 

*@see PingServer 

*@see BattleScene 

i 


public class BridgeServer{ 
public final static int PORT = 4129; 
public final static int REMOTE_PORT = 4130; 
public final static float SET_SONAR_PARAMETERS = (float)0.0; 
public final static float DETECT TARGETS - (float)1.0; 
public final static float DETECT_BEAM = (float)2.0; 
public final static float DETECT BOTH - (float)3.0; 
public static DatalnputStream remoteln = null; 
public static DataOutputStream remoteOut = null; 
public static DataInputStream in = null; 
public static PrintStream out = null; 
public final static String HOST = "electric.stl.nps.navy.mil”; // remote server host name 


ger 


* The main loop ofthe bridge program. It accepts input and output until its 
* network connections are dropped or the program is kılled by the user. No arguments 
* are expected. 


pee 


public static void main(String[] args)í 


// Declare Network specific variables to browser 
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ServerSocket server_socket = null; 
Socket client_socket = null: 
boolean r = true; 

float serviceRequest = (float)-1.0; 


// Network specific variables to remote host 
Socket remoteSocket = null; 


// Sonar Ping specific variables 
int i= 0, 

int j = 0; 

Vec3d pos = new Vec3d(); 


System.out.println("Start server: " + PORT), 


// Declare the socket to connect to BattleScene 

try{ 
server, socket 7 new ServerSocket( PORT), 

} catch(IOException e){ 
System.out.println("Could not create socket on: " + PORT +", " + e), 
System.exit(1), 

j 


System.out.printin("Socket created: "+ PORT); 


// request to open socket and data streams to the PingServer 
try { 
// open network and input/output stream 
remoteSocket = new Socket(HOST, REMOTE, PORT), 
remoteln = new DatalnputStream(remoteSocket.getInputStream()); 
remoteOut - new DataOutputStream(remoteSocket.getOutputStream()); 
) catch (UnknownHostException e) ( 
System.out println( "Unknown host: " * HOST); 
} catch (Exception e) { 
System.out.printin("Connection error"), 
j 


System.out.printin("Waiting for client..."); 


// accept BattleScene request for establishing a socket connection 
uy( 
client socket = server socket.accept(); 
) catch(IOException e) ( 
System.out.printin("Accept failed: " + PORT +", " +e); 
System.exit( 1); 
j 


System.out.printin("Connection established: " 
* client socket.getinetAddress()), 
System.out.println("Open input/output stream..."); 


// establish data streams to and from the BattleScene 
try{ 
in = new DatalnputStream(chent_socket.getInputStream()); 
out = new PrintStream(client_socket.getOutputStream()); 
y catch (IOException e) { 
System.out.printin("Could not create input/output stream on: " 
+ PORTE Te, 
System.exit(1 ); 
j 


// wax for a signal from client to begin 
System.out.printin("Reading data from client..."), 
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try { 


T = in. readBoolean(); 
remoteOut. writeBoolean(r); 
System.out.printin("Relayed a Start Request"); 


j 

catch (IOException e) { 
System.out.println("Could not read data." ); 
System.exit(1); 


j 


// returning data Environment data from PingServer to BattleScene 
try i 
String tempString = remoteln.readLine(), 
System.out.printin(" Relaying new VRML environment”); 
out.printIn(tempString), 


j 

catch (IOException e) 1 
System.out.printin("Could not read data."); 
System.exit(1); 

j 


// Start an infinite loop of accepting ping requests and returning VRML objects 
while(true){ 


System.out.printin("Reading data from client..."); 
try { 
serviceRequest = in.readF loat(); 
remoteOut.writeFloat(serviceRequest); 
if (serviceRequest == SET_SONAR_PARAMETERS) { 
setParameters(); 


j 

else if (serviceRequest -- DETECT TARGETS) ( 
System.out.printin("Relaying a Ping Request..."); 
String tempString = remoteln.readLine(); 
System.out printin(” Relaying detection field VRML..."); 
out. println(tempS tring); 

} 

else if (serviceRequest == DETECT BEAM) { 
System.out.println("Relaying a Ping Request..."); 
String tempString = remoteln.readLine(); 
System.out.println(" Relaying VRML Beam..."); 
out.printin(tempString); 

} 

else if (serviceRequest == DETECT_BOTH) { 
System.out printin("Relaying a Ping Request..."); 
String tempString = remoteln.readLine( ); 
System.out.println(" Relaying detection field VRML..."); 
out.printin(tempString); 
tempString = remoteln.readLine(); 
System.out.println(" Relaying VRML Beam..."); 
out.printin(tempString), 

j 


} catch (IOException e) { 
System.out.printin( "Could not read data"); 
System.exit(1 ); 

j 

j 


/* 
* This method sends the pertainant data to PingServer 
* from BattleScene 
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= 
private statıc void setParameters () { 
try { 

remoteOut. writeF loat(in.readFloat()); 
remoteOut.writeF loat(in.readF loat()); 
remoteOut. writeFloat(1n.readFloat()); 
remoteOut.writeFloat(1n.readF1oat()); 
remoteOut.writeFloat(in.readFloat()); 
remoteOut.wnteF loat(in.readF loat()); 
remoteOut. writeFloat(in.readFloat()); 
remoteOut. writeF loat(in.readF loat()); 
remoteOut.wniteFloat(in.readFloat()); 
remoteOut. wnteF loat(in.readF loat()); 
remoteOut. writeFloat(in.readF loat()), 


j 

catch (Exception e) ( 
System.out.println( "Could not read data" y; 
System.exit( 1), 


j 


C: PINGSERVER.JAVA 


Je 
File: PingServer.java 
Compiler: jdk1.1.6 
1 


package mil.navy.nps.rra; 


import java.net. *; 
import java.io. *; 


import mil.navy.nps.rra.*; 


/ Xx 

*@version 1.0 

*@author LT Timothy M. Holliday (<A HREF="http://www.stl.nps.navy.mil/~auv/holliday"> 
http://www. stl.nps.navy.mil/—auv/holliday</A>) 

x 

*<dt><b>Location:</b> 

*«dd»«a href-"http://www.stl.nps.navy.mildis-java-vrml/mil/navy/nps/rra/PingServer.]ava"7 

* http://www.stl.nps.navy. mil/dis-java-vrml/mil/navy/nps/tra/PingServer.java</a> 

* 


*<dt><b>Hierarchy Diagram:</b> 

*<dd><a href="images/RRAClassHierarchy. gif"><IMG SRC="images/RR AClassHierarchyButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Coordinate System Diagram:</b> 

*<dd><a href="images/CoordinateSystem.gif"><IMG SRC="images/CoordinateSystemButton.jpg" 
ALIGN=ABSCENTER></a> 

* 


*<dt><b>Summary:</b> 

*<dd>Takes requests for sonar pings and retums VRML representations of the sonar ping and 
* any detections made. 

* 

*<dt><b>Explanation:</b> 

*<dd>Since calculation of numerous ray trajectories is extremely time intensive, rendering and 
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* ray tracing often need to be on separate computers to improve performance times. This provides an interface 
from the rendering methods 
* to the ray tracing methods.<P> 


* 


*<dt><b>History:</b> 

*<dd> 15Jan98 /T Eu. M. Holliday /New 

*<dd> 17Mar98 /Timothy M. Holliday /Added HTML comment convention 
xk 

*(Qsee Lobe 

E, 


class PingServer( 


public final static int PORT = 4130; 

public final static float SET_SONAR_PARAMETERS - (float)0.0; 
public final static float DETECT_TARGETS = (float)1.0; 

public final static float DETECT_BEAM = (float)2.0; 

public final static float DETECT_BOTH = (float)3.0; 


public static void main(String[] args){ 


// Network specific variables 
ServerSocket server socket — null; 
Socket client socket = null; 
DatalnputStream ın = null; 
PrintStream out = null; 
boolean r = true; 
float serviceRequest - (float)-1.0; 


// Sonar Ping specific variables 
inti- 5; 
int j = 110; 
Vec3d pos = new Vec3d( -5000, 500, 0); 
SSP ssp - new SSP("traditional"y, 
Bottom bott = new Bottom("slope" 2000), 
Surface surf = new Surface("smooth"), 
Lobe lobel = new Lobe(); 
lobel.setDuration(1), 
lobel.setBottom(bott); 
lobel .setSurface(surf), 
lobel.setSsp(ssp); 


//(doubley, 10, (double)ı, 1.0, 10, 10, 
// pos, 2, bott, surf, ssp, 
// 4), 
Targets targets = new Targets(surf,bott, 70,0), 


// call each class’ VRML print routine 
System.out.println( "Start server: " + PORT); 


// open socket on PORT 

)تا 
server_socket = new serverSocke(PORT) |‏ 

) catch(TOException e)( 
System.out.printin("Could not create socket on: " + PORT +", " + e); 
System.exit(1); 

j 


System.out println("Socket created: "+ PORT). 
System.out.printin("Waiting for client...”); 


// accept client's request 


try{ 


client socket — server, socket.accept(); 
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} catch(IOException e) { 
System.out.printin("Accept failed: "+ PORT +", "+e); 
System.exit(1 ), 

j 


System.out.printin("Connection established: " 
* chent socket.getInetAddress()); 
System.out.println("Open input/output stream..."): 


try{ 
in = new DatalnputStream(client_socket.getInputStream()); 


out = new PnntStream(client_socket.getOutputStream()); 
} catch (IOException e) { 
System.out.printin( “Could not create input/output stream on: " 
+ PORT +e), 
System.exit( 1); 


j 
// wait for a signal from client to begin 

System.out.println( "Reading data from chient..."); 

try 1 
r = in. readBoolean(); 
System.out.printin("Received a Start Request"), 

j 

catch (IOException e) ( 
System.out.printin("Could not read data. "); 
System.exit( 1); 


j 


bott. setAppendLineFeed( false); 

surf. setAppendLineFeed( false ); 

System.out.printIn(" Sending new VRML environment"): 
String tempString = bott. VRML Bottom()+surf. VRML Surface( y; 
out.printin(tempString ); 


lobe 1 .setAppendLineF eed( false), 


while(true){ 
System.out.printin("Reading data from client..."); 


try { 
serviceRequest = in.readFloat(); 
if (serviceRequest == SET SONAR PARAMETERS) { 
setParameters(lobel, in, pos); 


j 

else if (serviceRequest -- DETECT. TARGETS) ( 
System.out.printin("Received a Ping Request..."); 
lobel.calculateLobe(targets); 
System.out.printin(" Sending detection field VRML ..."); 
out.println(lobel detection VRML()); 


j 

else if (serviceRequest = DETECT_BEAM) { 
System.out.println("Received a Ping Request..."); 
lobel.calculateLobe(targets); 
System.out.printin(" Sending VRML Beam..."); 
out.printin(lobel .dynamic VRMI ()); 


j 

else if (serviceRequest -2 DETECT BOTH) 1 
System.out.printin("Received a Ping Request...” ); 
lobel.calculateLobe(targets); 
System.out.printiIn(" Sending detection field VRML ..."), 
out.println(lobel .detection VRML()). 
System.out.printin(" Sending VRML Beam..."); 
lobel.setNumberY Pertition( 1 ); 
lobel .setNumberXPartition( 1); 
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lobel .reset(), 

lobel.calculateLobe(targets); 

out.println(lobel .dynamic VRMI ()); 
} 


+ catch (IOException e) { 
System.out.printin("Could not read data"), 
System.exit(1); 

$ 

j 
$ 


private static void setParameters (Lobe lobel, 
DatalnputStream in, 
Vec3d pos) { 
try { 
pos.set((double)in.readF loat(), 
(double)in.readFloat(), 
(double)in.readFloat()); 
float elevation = in.readFloat(); 
float azimuth = in.readFloat(); 
float verticalBeam Width = in.readFloat(); 
float horizontalBeam Width = in.readFloatt ), 
float verticalBeamConfiguration = in.readFloat(); 
float horizontalBeamConfiguration 7 in.readFloat(); 
float pingInterval = in.readFloat(); 
float powerLevel = in.readFloat( ); 
lobel.setElevation(elevation), 
lobel .setAzimuth(azımuth); 
lobel.setLobeWidthY(verticalBeamWidth ); 
lobel.setLobeWidthX(horizontalBeamWidth), 
lobel .setNumberY Partition((1nt)verticalBeamConfiguration ); 
lobel.setNumberXPartition((int)horizontalBeamConfiguration), 
lobel .setPosition(pos); 
lobe] .setEndTime(pingInterval); 
lobel .reset(); 
System. out.printin("Beam azimuth angle "+azimuth), 
System.out.printin("Beam elevation angle "+elevation); 


j 
catch (Exception e) { 
System.out.println("Could not read data"), 
System.exit( 1); 
j 
j 
j 
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APPENDIX D. TACTICAL VISUALIZATION POWER POINT SLIDE SET 


A. INTRODUCTION 


This appendix presents the annotated slide set for a Manta tactical visualization scenario. This 
slide set gives an overview of the Manta design, a presumed tactical scenario and how the scenario 
develops. This slide set is integrated with the Hvper Text Mark-up Language (HTML) and Virtual Reality 


Modeling Language (VRML) web pages shown in Appendix E. 
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Tactical Visualization 
of the Environment: 
Manta Minefleld Search 


Don Brutzman 
Naval Postgraduate School 


brutzman(a)nps.navy.mil 


September 97 


The purpose of this project is to examine a 
Manta minefield search mission. We will show 
how tactical visualization of the environment 
improves mission planning and understanding. 


You are reading a PowerPoint presentation. 
It corresponds to the 2D/3D slide show written 
in HTML (the Hypertext Markup Language) 
and VRML (the Virtual Reality Modeling 
Language). | 
Annotations for each slide describe MPEG 
video dialog and VRML scene viewpoints. 
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Manta Tactical Visualization 


W Goal: study tactical encounter in 2D & 3D 
B Produced by NPS faculty and students for 
Naval Undersea Warfare Center (NUWC) 
98 2D Hypertext Markup Language (HTML) 
> Links to other information placed in context 


m 3D Virtual Reality Modeling Language 
(VRML) scenes 


> Links, multiple viewpoints and animation in 3D 





Left video: Principal investigator. 


"This is a realistic scenario, both for 
submarines and for multiple Mantas. ' 


2 


VRML scene: Simple Manta model plus 
moving sonar beam and spatialized sound. 
Green is detection, red is counterdetection. 


Viewpoints: Manta close-up, over-the-shoulder 
view, looking back towards Manta sonar beam, 
Manta plus sonar from 1000m off track. 
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Credits 






NUWC sponsor Erik Chaum 
Principal investigator Don Brutzman 
Theatre Commander RADM Marsha Evans USN 
Thesis, OOD LT Tim Holliday USN 
Commanding Officer CDR Mike Holden USN 
Weapons Officer LT Ben McNeal USN 
Sonar Officer LT Kevin Byrne USN 
Intelligence Officer CPT Russell Storms USA 


Enk Chaum works in NUWC Code 22. 
He supervised the design and construction of the 
Warfare Systems Presentation Facility (WSPF). 


Don Brutzman is a computer scientist and assistant 
professor at NPS. His research interests include 3D 
graphics, underwater robotics and internetworking. 


RADM Marsha Evans USN is NPS Superintendent. 


(This is an optional link from the web page. 


There is no corresponding VRML scene.) 
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Manta Mission Overview 


B Intelligence report, deployment decision by 
theater commander 


m SSN arrives on station, Mantas are deployed 
B Minefield searched & mapped by 4 Mantas 
B Diesel sub found torpedo fire/counterfire 


B Enemy torpedo versus NPS low-cost 
acoustic torpedo countermeasure 


w Mission completion report, Mantas & SSN 





Left video: Principal investigator. 


“All of the robot searches and software shown 
in this mission can be demonstrated today. 

We are using existing capabilities of the NPS 
Phoenix AUV. 


VRML scene: A 25-meter Manta searching 
for a much smaller l-meter mine. A white 
. vertical marker helps viewers locate the mine. 


Viewpoints: Over-the-shoulder view, 
Manta plus sonar beam 1000m off search track, 
Manta close-up, looking back towards Manta 
sonar beam. 
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Manta Design Assumptions 


m Length: 80 feet. Endurance: several days. 


m Optimum search speed versus mines, 
bottomed diesel: 5 knots 


m Single Manta search rate: 1 km?/hour 

m Payload: sensors, lıght-weight torpedoes, 
countermeasures 

m Acoustic communications: reports, 
occasional coordination only 





Left video: Principal investigator. 


"A variety of potential Manta designs are being 
evaluated by NUWC. engineers. " 


VRML scene: Simple Manta model plus 
moving sonar beam and spatialized sound. 


Viewpoints: Manta Close-up, 
Manta plus sonar beam, Manta Aft Port View, 
Manta Fore Port View. 
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Manta Operational Scenario 


m Hostile aggression by Orange Nation 

m Amphibious landing in Tangerine Harbor 

m USS MONTEREY needed: Manta vehicles 

B Covertly map enemy harbor for mines, subs 

m Tactical problem is well-suited to multiple 
Manta search 

m Realistic scenario, cannot be performed 
covertly by today’s fleet 





Left video: Principal investigator. 


“The slides, scenes & videos in the rest of this 
presentation walk you through the entire Manta 
tactical encounter. 


VRML scene: A simple VRML sphere 
wrapped with an image texture of the Earth. 
Geographic maps of Tangerine Harbor and 
SUBLANT Norfolk are hidden inside. 


Viewpoints: Click Earth to see maps appear. 
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Intelligence Report 


m One diesel sub plus 100 mines protect the 
hostile harbor 


m Fleet commander preparing for amphibious 
assault 


w= Need rapid covert Manta mapping of 
minefield 


m Call sent to Manta-capable attack 
submarine: USS MONTEREY 


ا س m‏ 


Left video: Intelligence officer, reporting 
to Admiral. “Admiral, intelligence shows a 
diesel sub and 100 mines have been laid in 
Tangerine Harbor.” 


Right video: Theater commander. 


“That could delay the amphibious landing. Get 
me the skipper of the Monterey on the satellite 
link." | 


VRML scene: Simple Manta scene. 
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Manta Mission Tasking 


m Admiral calls SSN commanding officer via 
satellite link 
> Get on station 
> Prepare for Manta minefield search 
> Prepare for possible hostilities 


B CO acknowledges, SSN dives & transits 





Left video: Theater commander. 
“Skipper, we need you to search Tangerine 
Harbor with your onboard Mantas. We can’t 
afford any delay in the amphibious invasion. ” 


Right video: Commanding officer. “We re 
on our way, Admiral... Officer of the Deck, take 
us down and proceed at flank speed to the 
landing zone at Tangerine Harbor.” 


VRML scene: Highly detailed submarine at 
periscope depth. Mountains in background. 


Viewpoints: Multiple viewpoints: use 
PgUp/PgDn to walk camera around the sub. 
Click on periscope to lower it. 
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USS MONTEREY Arrives 


mw Theatre mission planning center prepares 
Manta mission package 


= SSN location: safe standoff distance 


= Periscope depth ops following rapid transit 
> Nav fıx, package download, intelligence update 


m Get sound speed profile (SSP), visualize 
sonar effectiveness 


m Verify mission packages, prepare for launch 


Left video: Intelligence officer. 


“The MONTEREY will only have 24 hours to 
map Tangerine Harbor. We'll evaluate search 
tactics using the virtual world.” 


VRML scene: 
Submarine proceeds to periscope depth (PD). 
Viewpoints: 


Sub 500m away, scope up at periscope depth. 
Close-up: sub approach to periscope depth. 
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Mantas Away 





B Four Mantas launch from submarine 


= SSN stays at max communication range to 
remain undetected 


m Mantas complete transit to minefield and 
update SSP 


mg Mantas commence coordinated search 


m Tactical analysis developed by NPS 
Undersea Warfare officers 


Left video: Officer of the Deck (OOD). 
"Captain, we re on station, and the Manta 
mission package download is in progress. " 


Right video: Commanding Officer (CO). 
“Weps, you may launch Mantas when ready. ” 


Weapons Officer (Weps). 
“The boat is maintaining constant depth and 
speed... Mantas away. " 


VRML scene: Two of four Mantas are 
shown separating from the forward hull and 
proceeding to the search area. 


Viewpoints: Various camera angles. 
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Manta Minefield Search 


mw Visualize large-scale minefield search 


E Visualize small-scale Manta sensor 
interactions 


= Monitor progress of long-term search in 
three dimensions 


B All real-world data has corresponding place 
in virtual world 


m Expect better operator sense of presence 


Left video: Weapons Officer (Weps). 
"All Mantas are away, all communication links 
are up, and the Mantas are finding mines. ” 


VRML scene: This ıs an extended (24 hour) 
minefield search, reenacted by playing back the 
Manta mission scripts in the 3D VRML scene. 
After 100 mine boxes appear, Manta will start. 


Viewpoints: A variety of camera angles 
allow big-picture visualization ofalarge (10 
kilometer by 10 kilometer) harbor, as well as 
close-up views of Manta-mine interactions. Red 
bounding boxes & spheres help viewers | 
visualize the search and localize small mines. 
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Diesel Sub Found 


m Diesel sub on bottom in minefield, awaiting 
landing force 

m Manta reports diesel location to SSN and 
other Mantas 

B Manta requests permission to fire 

m Covert rules of engagement: self defense only 


m If permission to fire not given, the diesel sub 
engages Manta 





Left video: Officer of the Deck (OOD). 
"Captain, Manta Alpha has found a diesel 
submarine and requests permission to fire.” 
Commanding Officer (CO). 
"Permission denied. . Our current rules of 
engagement only allow firing in self defense. ’ 
VRML scene: The Manta passes over a 
hostile diesel submarine hiding on the bottom. 


7 


Viewpoints: Various aspects relative to 
submarine and Manta enable visualizing the 
encounter. 
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Fire / Counterfire / 
Torpedo vs. Countermeasure 


————————— ÓGÓ————————————————— ——————  ''OOOCoÓÓ—— À' 


m Diesel sub fires torpedo at Manta 
B Manta counterfires torpedo at diesel, evades 
m Manta launches low-cost countermeasure 


monitoring by SSN during engagement 
mg Low-cost DSP-based acoustic anti-torpedo 


m Manta communication relay enables covert | 
countermeasure: ongoing NPS work for ONR 


سس س ہے تسس — u‏ سس سح 


Left video: Sonar: “Torpedo in the water, 
launched from the diesel submarine. ” 


Weps: “Second torpedo counterfired by Manta.’ 
Sonar: “Manta is evading with countermeasures. 


OOD: “Manta’s torpedo is in terminal homing, 
on the diesel. ” 


Sonar: "Loud explosion down the bearing of the 
diesel submarine. " 


VRML scene: This one-quarter project shows 
embedded sound ın a torpedo-countermeasure 
interaction. Click on the torpedo to start the run. 


Viewpoints: Note NTDS symbols and 
equipment labels embedded 1n the environment. 
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Manta Mission Complete 





m Mantas destroy diesel, evade counterfire 
and map minefield 


m SSN reports mission complete to theatre 
commander 


Left video: Commanding Officer (CO). 
“Admiral, the minefield at Tangerine Harbor 
has been mapped, one diesel submarine has 
been destroyed, and all four Mantas are 
recovered & back on board." | 


Risht video: Theater commander. 


“Naval Postgraduate School student officers 
and faculty can make big contributions to real 
fleet challenges.” 


VRML scene: Manta search pings. 


Viewpoints: Over-the-shoulder viewpoint. 
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Conclusions 


m Mantas might perform essential missions 


m Visualization conclusion: add networked 
3D graphics, vehicles, physics and sensors 
just like any other Web-based content 


m "Building content" is better, more scalable 
than "programming" 


m NPS ready (and eager!) to continue 
work for NUWC 


Our conclusions match the expected results of 
the NUWC-NPS research proposal. 


These dramatic results show that construction of 
composable physics-based 3D scenes 1s possible 
using VRML. Tactical visualization of the 
environment improves understanding for forces 
afloat and the engineering community ashore. 


The star of this show is 3D visualization. 

We have demonstrated the viability of 2D/3D 
tactical visualization of the environment. A lot 
of promising work lies ahead. 
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Contact information 


Don Brutzman 


brutzman(@nps.navy.mil 
http://www.stl.nps.navy.mil/—-brutzman 


Code UW/Br Naval Postgraduate School 
Monterey California 93943-5000 USA 
408.656.2149 voice 
408.656.3679 fax 
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APPENDIX E. TACTICAL VISUALIZATION HTML DOCUMENTS 


A. INTRODUCTION 


This appendix presents Hyper Text Mark-up Language (HTML), MPEG-2 digitized video and 
Virtual Reality Modeling Language (VRML) web pages which show how HTML and VRML may be 
combined to form a tactical scenario presentation system. One advantage to using VRML in the 
presentation system is that behaviors can be given to the objects in the scene. This allows the scene to 


have a time dependent nature, instead of being just a static picture. 
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APPENDIX F. RRA API JAVADOC DOCUMENTATION 


A. INTRODUCTION 


This appendix gives a sample of RRA API documentation produced by javadoc, a Sun 
Microsystems product. This program comes with the Java distribution. This program can process a Java 
package (collection of Java classes) and produce detailed and standardized reference documents. These 
documents come in the Hyper Text Mark-up Language format and thus can be read by any standard web 
browser. This documentation generation program allows fast development of a cohesive API, since all 


parties can have rapid access to the most up-to-date API documentation. 
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All Packages Index - 





Class Hierarchy 


e class java.lang.Object 
o class mil.navy.nps.rra.Beam 
class mil.navy.nps.rra.Bottom 
class mil.navy.nps.rra.ExampleBeamDynamic 
class mil.navy.nps.rra.ExampleBeamStatic 


class mil.navy.nps.rra.ExampleLobeDynamic 
class mil.navy.nps.rra.ExampleLobeStatic 


class mil.navy.nps.rra.ExampleRay 
class mil.navy.nps.rra.Lobe 

class mil.navy.nps.rra.Print VRML 
class mil.navy.nps.rra.Ray 

class mil.navy.nps.rra. SVP 

class mil.navy.nps.rra. Surface 
class mil.navy.nps.rra. Targets 
class mil.navy.nps.rra. Vec3d 


000000000000 0 
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All Packages Class Hierarchy This Package Previous Next Index 








Class mil.navy.nps.rra.Beam 


3ava.lang.Object 
| 


+=-=mıl.navy.nps.rra.Beam 


public class Beam 
extends Object 


Author: 
LT Timothy M. Holliday ( http://www.stl.nps.navy.mil/-auv/holliday) 
Location: 


http://www.stl.nps.navy.mil/dis-java-vrml/mil/navy/nps/rra/Beam.java 


Hierarchy Diagram: 
LM ETE أ۱.‎ 
Coordinate System Diagram: 
ITIN 


Summary: 
Forms a beam from four rays, propagates the beam and can return a VRML 
representation of the beam. 

Explanation: 


Beams of energy can be considered to consist of bundles of rays. The energy in a 
bundle can be shown not to diverge from the bundle. Thus the energy in a bundle 
is constant and thus the product of the intensity and the area of the bundle 
perpendicular to direction of propagation is a constant as well. Therefore a beam 
is a fundemental building block for a lobe of a sonar pattern. 


History: 


15Nov97 /Timothy M. Holliday /New 

17Mar98 /Timothy M. Holliday /Added HTML comment convention 
12Apr98 /Timothy M. Holliday /Parameterless Constructors 
14Apr98 /Timothy M. Holliday /Simplified VRML Routines 
21Apr98 /Timothy M. Holliday /Fixed Problem with 
calculateSoundPressureLevel 


See Also: 


Ray, BeamExampleStatic, BeamExampleDynamic, Bottom, Surface, Vec3d 
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Variable_Index 


e NONE 

TL 

ə TIME 

Constructor /ndex 
» Beam() 


Constructor for the Beam class. 
Method. Index 


» calculateBeam(Targets) 
This method calculates the trajectory of the beam of energy enclosed by the 
defining rays of the beam tube. 

e calculateSoundPressureLevel(double[][ ][], double, double, double) 
This method calculates the trajectory of the beam of energy enclosed by the 
defining rays of the beam tube. 

e dynamic VRML() 
This method creates a dynamic VRML string shape that 1s the three dimensional 
representation of the beam pulse that is propogated. 


e getAppendLineFeed() 
This is a static method that returns the current line appendage. 


e getAzimuth() 
This method returns the azimuthal angle. 


e getBeamNumber() 
This returns BeamNumber. 


e getBottom() 
This method returns the handle to the bottom object. 


» getDeltaTime() 
This method returns the simulation time step. 


e getDetectCount() 
This method returns the total number of detects by the beam 


» getDetectEchoLevel(int) 
This method returns the echo level of the detected target 


e getDetectTime(int) 
This method returns the detect time for a given detect in the beam 


» getDuration() 
This method sets the duration of the sonar pulse. 


» getElevation() 
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This method returns the elevation angle. 


e getEndTime() 
This method returns the simulation end time. 


e getHalfBeamWidthX() 
This method returns the half beam width of the beam in the azimuthal direction. 


e getHalfBeamWidthY() 
This method returns the half beam width of the beam in the azimuthal direction. 


e getPosition() 
This method returns the position of the beam. 


* getSsp() 
This method returns the handle to the sound speed profile object. 


e getSurface() 
This method returns the handle to the surface object. 


e pingTimerVRML/) 
This method returns a VRML timer string that contains the appropriate timing 
information for the beam 
e reset() 
This method resets all ofthe beam parameters after instanciation has occurred 
since reuse 1s more time efficient than garbage collection and reallocation. 
» setAppendLineFeed(boolean) 
This is a static method used to indicate whether a line feed is desired at the end 
of every line. 
e setAzimuth(double) 
This method sets the azimuthal angle, which is the angle from the x-axis to the z- 
axis rotating about the y-axis. 
« setBeamNumberfint) 
This method sets BeamNumber. 
e SetBottom(Bottom) 
This method sets the handle to the bottom object. 
e setDelta Time(double) 
This method sets the time step in the simulation. 
e setDuration(int) 
This method sets the duration of the sonar pulse. 
e SetElevation(double) 
This method sets the elevation angle, which is the angle from the y-axis to the x- 
axis rotating about the z-axis . 
e set£ndTime(double) 
This method sets the simulation end time. 
e setHalfBeamWidthX(double) _ 
This method sets the half beam width of the beam in the azimuthal direction. 
e setHalfBeamWidthY(double) 
This method sets the half beam width of the beam in the elevation direction. 
e setPosition(double, double, double) 
This method sets the position of the beam. 
e setSsp(SSP) 
This method sets the handle to the sound speed profile object. 


e setSurface(Surface) 219 


This method sets the handle to the surface object. 

e Static VRML(nt, int) 
This method writes to the console a VRML shape that is the three dimensional 
representation of the beam that is propagated. 


variables 


eTL 
public static fine! کال‎ ۵ 


9 TIME 


public static fina" rnt TIME 


9 NONE 


public static final int NONE 


COMSÉLUCÉOLS 


9 Beam 


public Beam() 


Constructor for the Beam class. A beam is defined as the volume swept out by 
four rays as they traverse the ocean environment 


u -< 


9 reset 


public void reset() 


This method resets all of the beam parameters after instanciation has occurred 
since reuse is more time efficient than garbage collection and reallocation. 


9 setAzimuth 
public void setAzimuth(double phi) 


This method sets the azimuthal angle, which is the angle from the x-axis to the z- 
axis rotating about the y-axis. 


i 220 


9 getAzimuth 
public double getAzimuth () 
This method returns the azimuthal angle. 
9 ۷/٥11۸ 
public void setHalfBeamWidthX (double pHalfBeamWidthX) 
This method sets the half beam width ofthe beam in the azimuthal direction. 
9 zetHalfBeamWidthX 
public double getHalfBeamWidthX() 
This method returns the half beam width of the beam in the azimuthal direction. 
9 setElevation 
public void setElevation(double beta) 


This method sets the elevation angle, which is the angle from the y-axis to the x- 
axis rotating about the z-axis . 


9 zetElevation 
public double getElevation() 
This method returns the elevation angle. 
9 setHalfBeamWidthY 
public void setHalfBeamWidthY (double pHalfBeamWidthY) 
This method sets the half beam width of the beam in the elevation direction. 
9 gzetHalfBeamWidthY 
public double getHalfBeamWidthY() 
This method returns the half beam width of the beam in the azimuthal direction. 


@ setPosition 


public void setPosition(double x, 
double y, 
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double z) 
This method sets the position of the beam. 
@ getPosition 
public Vec3d getPosition() 
This method returns the position of the beam. 
9 setDeltaTime 
public void setDeltaTime (double pDeltaTime) 
This method sets the time step in the simulation. 
Ə setDeltaTime 
public double getDeltaTime () 
This method returns the simulation time step. 
9 setEndTime 
public void setEndTime (double pEndTime) 


This method sets the simulation end time. This value is reletive to the start time 
which is 0.0. 


9 getEndTime 
public double getEndTime () 
This method returns the simulation end time. 
9 setDuration 
pubie vord setburatjon(rnt duratron) 


This method sets the duration of the sonar pulse. The integer refers to the 
number of deltaTime increments. currently only 1 or 2 are allowed. 


9 getDuration 
public int getDuratčion() 


This method sets the duration of the sonar pulse. The integer refers to the 


number of deltaTime increments. 
222 


9 


public void setBottom(Bottom pBottom) 
This method sets the handle to the bottom object. 
@ getBottom 
public Bottom getBottom() 
This method returns the handle to the bottom object. 
@ setSurface 
public void setSurface (Surface pSurface) 
This method sets the handle to the surface object. 
9 getSurface 
public Surface getSurface() 
This method returns the handle to the surface object. 
9 setSsp 
public void setSsp(SSP pSsp) 
This method sets the handle to the sound speed profile object. 
9 getSsp 
public SSE cget>sspi) 
This method returns the handle to the sound speed profile object. 
9 setBeamNumber 
public void setBeamNumber (int number) 


This method sets BeamNumber. It is used to uniquely identify each beam for 
ROUTEing in VRML. 


@ getBeamNumber 
public int getBeamNumber () 
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This returns BeamNumber. Which is used to uniquely identify each beam for 
ROUTEing in VRML. 


9 static VRML 


public String staticVRML(1int colorcborce, 
int intensityChoice) 


This method writes to the console a VRML shape that is the three dimensional 
representation of the beam that is propagated. 


9 dynamic VRML 
public String dynamicVRML() 


This method creates a dynamic VRML string shape that is the three dimensional 
representation of the beam pulse that is propogated. 


9 calculateBeam 


public void calculateBeam(Targets targets) 


This method calculates the trajectory of the beam of energy enclosed by the 
defining rays of the beam tube. It also determines when there is a detection. 


@ calculateSound Pressure Level 


public void calculateSoundPressureLevel (double field[]i]Il, 
double deltaRange, 
double deltaDepth, 
double frequency) 


This method calculates the trajectory of the beam of energy enclosed by the 
defining rays of the beam tube. 


9 getDetectCount 


public int getDetectCount () 
This method returns the total number of detects by the beam 


9 getDetectTime 


public double getDetectTime(int N) 


This method returns the detect time for a given detect in the beam 


224 


9 setDetectEchoLevel 


public double getDetectEcholevel (int N) 
This method returns the echo level ofthe detected target 
9 pingTimerVRML 
public Stringi ping IMSrVRML () 


This method returns a VRML timer string that contains the appropriate timing 
information for the beam 


# setAppendLineFeed 


public static void setAppendLineFeed (boolean pAppendLineFeed) 
This is a static method used to indicate whether a line feed is desired at the end 
of every line. 'true' indicates a linefeed is desired and 'false' indicates that a space 


1s desired" 


2 getAppendLineFeed 


public static boolean getAppendLineFeed() 


This is a static method that returns the current line appendage. 


All Packages Class Hierarchy This Package Previous Next Index 
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APPENDIX G. RRA VIDEO INFORMATION 


A. INTRODUCTION 


This appendix provides information on how to obtain a copy of a video that shows the abilities of 
the RRA API and the simulations and visualizations obtained from it. The point of contact for obtaining a 


copy is Dr. Don Brutzman, brutzman@nps.navy.mil. 
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APPENDIX H. RRA CD ROM INFORMATION 


A. INTRODUCTION 


This appendix provides information on how to obtain a copy of a CD ROM that has this thesis, 
the RRA API source code, example programs using the RRA API, the simulations and visualizations 
obtained from using the RRA API and the contents of Appendix D, Appendix E and Appendix F. The 


point of contact for obtaining a copy is Dr. Don Brutzman, brutzman@nps.navy.mil. 
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